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IMPACT  OF  AIRCRAFT  EMISSIONS  ON  AIR  QUALITY 
IN  THE  VICINITY  OF  AIRPORTS 

VOLUME  III:  AIR  QUALITY  AND  EMISSION  MODELING  NEEDS 

by 

Donald  M.  Rota 


1  INTRODUCTION 


Estimates  of  emission  rates  and  atmospheric  concentrations  of  various 
airborne  substances  In  and  around  aviation  facilities  are  required  for 
determining  compliance  vlth  regulations,  testing  the  efficacy  of  potential 
control  strategies,  and  a  vide  variety  of  other  purposes. 

Precisely  hot*  these  estimates  should  be  made,  that  Is,  what 
computational  procedures  should  be  used,  has  been  the  subject  of  both 
theoretical  and  field  studies  for  over  a  decade*  When  the  Airport  Vicinity 
Air  Pollution  (AVAP)  and  the  Air  Quality  Assessment  Model  (AQAM)  were  first 
designed  some  10  years  ago,  for  commercial  and  military  facilities  respec¬ 
tively,  the  major  emphasis  was  on  providing  a  user-oriented,  state-of-the-art 
tool  that  would  allow  the  user  to  treat  virtually  all  sources  of  pollution 

with  as  much  detail  as  he  saw  fit.  The  computer  codes  were  developed  and 

tested  on  large,  relatively  fast  mainframe  computers*  It  was  assumed  that  the 
users'  machines  would  be  similar  and  that  the  users  would  have  staff  committed 
to  the  maintenance  and  operation  of  such  complex  computer  codes  on  a  long-term 
basis*  Such  staff  were  expected  to  be  familiar  with  the  aviation  facilities, 
air-pollution  regulations,  and  the  principles  of  emission  and  alr-quallty 

computation* 

Experience  has  shown  that  these  user  requirements  were,  in  fact,  too 
restrictive  to  permit  as  wide  use  as  might  otherwise  be  the  case*  Even  though 
computer  use  has  greatly  expanded  over  the  years,  such  of  this  growth  has  been 
In  the  area  of  the  new,  inexpensive  microcomputers  and  minicomputers*  Such 
machines,  while  convenient  and  eaay  to  use,  are  not  fast  enough  nor  do  they 

possess  the  necessary  memory  capacity  to  handle  the  AVAP  and  AQAM  computer 
codes  In  their  present  configurations*  Hence,  the  use  of  codes  like  these  has 
neither  grown  with  the  use  of  computers  nor  with  the  need  to  perform  emission 
or  alr-quallty  calculations*  Therefore,  along  with  the  need  to  maintain  these 
early  codes  and  to  periodically  perform  technical  updates  and  refinements  to 
keep  pace  with  developments  in  the  stats  of  the  art  of  modeling,  there  is  also 
the  clear  need  to  design  new  versions  that  meet  ths  needs  and  constraints  of  a 
greater  number  of  users  operating  in  the  new  environment  of  the  microcomputers 
and  minicomputers* 

The  purpose  of  the  present  report  is  to  address  the  subject  of  modeling 
needs  and  how  those  needs  can  be  systematically  satisfied*  It  is  noted  that 
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while  e owe  of  these  need#  have  already  been  wet  with  existing  models  and  their 
corresponding  computer  codes,  others  remain  to  be  satisfied.  The  approach  to 
this  subject  used  here  Is  essentially  appllcations-orlented.  It  begins  with  a 
listing  of  the  types  of  problems  Involving  aviation  facilities  that  are  likely 
to  require  the  use  of  computer  models  and  then  examines 9  In  some  detail,  the 
technical  characteristics  of  these  problems  that  determine  which  features  an 
applicable  model  oust  possess.  The  position  Is  taken  that  It  is  the  nature  of 
the  application,  together  with  various  user  and  machine  requirements,  that 
dictates  the  type  of  model  to  be  used.  Given  this  position,  the  ability  of 
the  currently  available  versions  of  the  AVAP  and  AQAM  computer  codes  to 
satisfy  these  technical  and  operational  requirements  Is  examined  In  detail. 
The  strengths  and  weaknesses  of  these  models  and  their  computer  codes  and 
supporting  documentation  are  noted.  Various  necessary  and  desirable  changes 
to  these  models  to  make  them  applicable  to  a  greater  variety  of  problems,  to 
update  them,  and  to  make  them  more  usable  are  then  discussed. 

As  a  guide  to  possible  future  technical  Improvements,  the  procedures 
and  algorithms  used  by  the  AVAP  and  AQAM  computer  codes  are  compared  and  the 
best  procedures  are  identified.  Where  necessary,  alternative  procedures  are 
recommended. 

A  "decision  tree"  Is  then  presented  as  a  systematic  means  of  laying  out 
some  of  the  alternative  courses  of  action  that  face  the  decision  makers  in 
trying  to  determine  how  best  to  resolve  the  discrepancies  that  exist  between 
the  present  modeling  needs  and  the  current  modeling  capabilities,  weaknesses, 
and  limitations.  Both  separate  and  joint  civilian  and  military  agency 
alternative  actions  are  presented. 

Finally,  the  modeling-needs  discussion  ends  with  a  presentation  of  a 
design  of  a  new  computational  system  that  la  proposed  to  satisfy  one  of  the 
more  important  new  user  requirements  not  met  by  any  current  version  of  AVAP  or 
AQAM.  The  new  system  is  designed  to  be  implemented  on  a  small  computer,  which 
should  greatly  increase  its  usability  by  in-house  agency  staff. 

The  report  is  divided  into  11  sections.  The  introductory  section  Is 
followed.  In  Sec.  2,  by  e  brief  summery  of  types  of  potent Isl  problems 
requiring  the  use  of  computer  models.  Section  3  examines  those  technical 
characteristics  of  the  various  model  applications  that  actually  determine  what 
featuraa  an  appllcabla  modal  should  poaaaaa.  The  operational  requirements  for 
models  are  dlecuseed  in  Sec.  A. 

The  AVAP  model  end  ite  various  currently  available  versions  ere 
described  in  See.  3,  along  with  e  summery  end  evaluation  of  its  intended  usee, 
strengths,  limit at lone,  end  weaknesses.  Updates  needed  to  Improve  the  model 
are  also  presented.  The  AQAH  and  ite  various  versions  ere  summarised  in  Sec. 
6  using  the  seme  format  ee  used  for  the  AVAP  model.  Section  7  contains  e 
detailed  f set u re -by-feature  comparison  of  the  AVAP  model  end  AQAM  in  which  the 
beet  approaches  or  alternatives  ere  identified. 


The  "decision  tree"  is  presented  in  Sec*  8,  end  Sec*  9  is  devoted  to  an 
outline  of  the  proposed  new  Joint  computational  system  to  be  Implemented  on 
small  computers*  Section  10  is  a  report  on  a  sensitivity  study  of  emissions 
calculations*  It  is  intended  to  provide  some  insight  into  the  importance  of 
including  various  source  categories  in  an  emission  Inventory  of  an  air  base  or 
airport*  Finally,  Sac*  11  contains  a  summary  and  recommendations  for  future 
tasks  needed  to  Improve  the  current  models  from  both  the  technical  and  user 
points  of  view* 


2  POTENT l AL-PROBLKM  TYPE8  AND  MODEL  USES 


The  two  mIq  problem  of  air-quality  analysis  ere  estimation  of 
emission  rates  and  computation  of  atmospheric  concentrations .  For  some 
purposes  It  la  sufficient  to  simply  estimate  the  emission  rates,  while  other 
problems  require  computation  of  the  atmoapherlc  concentrations  as  well* 
Examples  of  problems  requiring  emission  estimates,  but  not  necessarily 
atmospheric  concentrations,  are: 

e  Determining  compliance  with  state  Implementation -plan  (SIP) 
regulations  or  emission  limits* 

e  Using  rough-cut  or  screening  procedures  to  determine 
whether  a  potential  alr-quallty  problem  exists  or  could 
exist  and  whether  a  more  detailed  analysis  Is  required* 

e  Estimating  the  overall  effectiveness  of  emission-control 
strategies. 

Civilian  and  military  problems  requiring  computation  of  emission  rates 
and/or  atmospheric  concentrations  of  various  substances  are  generally  either 
of  a  regulatory  or  research  nature*  There  are  also  a  number  of  specialized 
problems  that  pertain  only  to  the  military  sector.  These  latter  problems  are 
essentially  beyond  the  scope  of  the  present  discussion  and  are  therefore  only 
briefly  referred  to  for  completeness*  Some  potential  model  uses  falling  under 
the  three  main  categories  are  listed  below. 


2*1  GENERAL  RESEARCH  PROBLEMS 

Investigation  of  potential  health  and  environmental  Impacts  of 
pollutant  emissions  from  various  source  types  (receot  Interest  has  focused  on 
photochemlcally  reactive,  odor  causing,  and  hazardous  airborne  gases  and 
aerosols)*  Some  possible  model  uses  are: 

e  An  aid  to  assessing  the  nature  and  spatial  extent  of  a 
particular  impact  or  hazard* 

a  An  aid  to  alr-quallty  and  emission  data  analysis  and 
diagnosis  of  emission-inventory  errors  and  deficiencies. 

a  Delineation  of  physical  and/or  chemical  effects  when 
complex  mixtures  of  phenomena  are  involved* 


e  Development  or  refinement  of  special  dispersion  algorithms 
or  submodels* 


5 


•  Development  and  Casting  of  simplified  or  "field"  models. 

a  Identification  of  eaieslon  or  air-quality  "hoc  spots." 

a  Aid  In  dsslgn  of  flald  taats  or  experiments  (a.g. , 
monitoring-instrument  daployaant). 


2.2  GENERAL  REGULATORY-RELATED  PROBLEMS 

a  Aid  In  perforalng  ragulatory-lapact  analyses. 

s  Detsralnatlon  of  need  for  and  efficacy  of  eaisslon  redac¬ 
tions. 

s  Evaluation  of  aalsslon-control  strataglaa. 

a  Preparation  of  preconstruction  anvlronaental  lapact  atate- 
aants  (ElSa)  for  new  aircraft  facilities  or  major 
aodlf lcatlons  to  existing  facilities. 

a  Screening  to  separate  cases  that  are  not  probleas  froa 
those  requiring  sore  detailed  treataent. 

e  Testing  for  coapllance  with  local,  state,  and  federal 
regulations  (e.g.,  national  aablent  air  quality  standards 
[NAAQSs] ,  SIPs,  and  prevention  of  significant  deterioration 
[PSD)  regulations). 

e  Calculations  related  to  legal  actions,  dtlxen  coaplalnts, 
etc. 

e  Analysis  of  lapllcatlons  of  proposed  new  regulations  for 
aircraft  operations,  fuel  use,  engine  design,  thrust 
settings,  etc. 


2.3  SPECIALIZED  MILITARY  APPLICATIONS 

e  Investigations  of  dispersion  of  gases  and  aerosols  of 
■liltary  laportance.  Including  obscurants,  sprays,  etc..  In 
both  terrestrial  and  aarltlne  environments. 

e  Investigations  of  alrcraft-pluas  dispersion  In  relation  to 
detection  using  elect ro-optleal  sensors. 
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•  Investigation  of  the  dispersion  of  effluents  froa  routine 
operations  at  rocket- launch ,  rocket-sled,  end  engine  test 
facilities* 

e  Investigation  of  accident  scenarios* 
e  Evacuation-corridor  or  hazard-zone  prediction* 
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3  TECHNICAL  REQUIREMENTS  FOR  MODELS 

In  view  of  Che  substantially  different  cechnlcel  requirements  for 
emission-rate  end  atmospheric-concentration  computations,  It  le  worthwhile  to 
exealne  theee  requirements  separately. 

3.1  SOURCE-EMISSION- INVENTORY  COMPUTATION  REQUIREMENTS 

The  eoaplletlon  of  e  source  emission  Inventory  tends  to  be  stralght- 
forwsrd  but  manpower  Intensive.  The  eaount  of  asnpower  required  Is 
essentially  directly  proportional  to  the  number  of  sources  to  be  considered. 
For  each  Individual  source,  the  emission  rate  is  generally  the  product  of  an 
appropriate  ealselon  factor  (expreased  in  terms  of  mass  of  emissions  per  unit 
tine  per  unit  of  source  activity)  end  a  measure  of  source  activity.  For 
automobiles,  the  unit  of  activity  is  vehicle  miles  (number  of  vehicles  times 
number  of  idles  traveled  by  esch  vehicle).  For  a  single  aircraft,  the  unit  of 
activity  Is  fuel  flow  rate  at  a  particular  thrust  setting.  In  an  airport  or 
air  baae,  aircraft  activity  generally  refers  to  the  number  of  aircraft  of  a 
given  type  In  a  given  mode  of  operation  (l.e.,  having  a  given  thruet 
setting).  In  practice,  because  the  emission  factors  for  each  pollutant  depend 
on  the  particular  type  of  source  as  well  as  on  its  mode  of  operation  and  fuel 
type,  the  enieslon  rate  calculations  can  become  far  too  time-consuming  and 
tedious  to  be  performed  by  hand,  eapedally  when  the  source  inventory  is 
large.  Furthermore,  if  a  large  number  of  tedious  hand  calculations  are 
required,  the  chances  of  making  errors  are  large  and  the  effort  required  to 
check  for  and  correct  such  errors  is  almost  as  laborious  as  the  original  hand 
calculations  themselves.  Bence,  except  in  cases  where  only  one  or  two 
isolated  sources  ere  being  considered  or  where  only  net  changes  in  total 
emissions  due  to  a  few  changes  in  selected  sources  are  being  considered.  It  Is 
highly  desirable  to  have  a  source -emission -Inventory  computer  model.  The 
principal  technical  features  required  of  such  a  model  are  listed  below. 

1 .  Up-to-date  emission  factors  for  all  the  criteria 
pollutants  should  be  Included  for  all  of  the  sources 
likely  to  be  encountered  at  an  aircraft  facility.  These 
emission  factors  should  preferably  be  based  on  actual 
source-emission  measurements  or,  lees  preferably,  on 
engineering  estimates,  and  should  be  adjustable  but 
protected  from  inadvertent  adjustment  by  general  users. 

2.  Emission  ratss  for  each  source  type  should  be  calculated 
and  stored  so  an  hourly  basis  for  subsequent  processing  or 
for  use  by  dispersion  models* 
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3*  It  aay  also  bt  desirable,  strictly  for  accounting  pur- 
poaaa ,  to  coapute  cal aa Iona  on  an  annual  or  other  long- 
tcra  basic.  For  such  purposes,  landing  and  takeoff  (LTO)- 
eyele-type  ealsslon  factors  should  be  eaployed. 

4.  For  purposes  of  providing  Input  data  for  dispersion  calcu¬ 
lations,  It  la  necessary  to  provide  hourly  ealsslona  for 
each  physical  source  according  to  the  source  geoaetry 
(point,  line,  area). 

5.  Since  civilian  and  allltary  facilities  use  different 
operetlng  procedures  and  service  different  aircraft  types. 
It  la  probably  aost  efficient  to  have  separete  but 
elallarly  structured  ealsslon  aodels  for  the  two  types  of 
facilities. 

6*  In  view  of  the  fact  that.  In  aoae  cases,  pollutant  concen¬ 
trations  aay  not  be  required  and  In  order  to  reduce  the 
overall  coaputer  requlreaents ,  It  would  be  best  to  keep 
the  ealsslon  aodels  separate  froa  the  dispersion  aodels. 

7.  Provisions  should  be  available  for  the  user  to  Include 
t rest went  of  aircraft,  aircraft  support  or  service 
equlpaent,  access  vehicles,  and  other  facility  subtle  and 
stationary  sources.  The  user  should  have  the  option  to 
Ignore  any  coablnatlon  of  these  source  types. 

8*  The  spatial  resolution  require aent  of  eource-ealsslon 
Inventories  depends  on  whether  dispersion  calculations  are 
required.  For  ealsslon  calculations  only,  the  spatial- 
resolution  requlreaents  are  gsnarally  aintaal  or 
nonexistent.  For  dispersion  calculations,  especially  of 
the  Gauss lan-aodel  type,  the  spatial  resolution  aust  be 
Increased  as  the  source-receptor  distances  of  Interest 
decrease.  Generally,  sines  public  access  Is  Halted  to 
certain  areas  of  the  airport,  the  spatial  resolution 
requireaent  is  not  the  seas  for  all  sources. 

9.  Airborne  aircraft  sources  are  generally  regarded  as 
Insignificant  coepared  to  ground-baaed  aircraft  sources 
for  ground-level  pollutant  concentrations.  However,  for 
purposes  of  souree-eaisalon-lnventory  coaplateneas,  or  to 
test  the  significance  of  such  sources,  the  user  should  be 
provided  with  the  option  of  including  or  oaittlng  thea. 
In  either  ease,  the  spatlal-teoolutloa  requireaent  for 
such  sources  Is  not  high,  (See  See.  10.3.) 


/ 


9 


10.  Inclusion  of  mporttln  ooooethaae  hydrocarbons  (MNC) 
sources  should  bo  optlonnl.  Ths  spstlsl- resolution 
requirement  for  such  sources  Is  afalmel.  Ibis  follows 
for  two  rsssone.  First,  those  sources  generally 
contribute  only  s  earn  11  fraction  of  the  total  MK 
snlsslons  fron  all  sources  (see  Sec*  10.$).  Second,  the 
NWiC  pollutant  category  Is  not  governed  by  e  health- 
related  anblent-alr  quality  standard  and  therefore  does 
not  require  detailed  source-receptor  dispersion-model 
calculations.  On  the  other  hand,  for  purposes  of 

enlsalon-lnventory  conpleteness,  or  to  assess  their 
potential  for  contributing  to  photochemical  snog,  their 
Inclusion  asy  be  warranted. 

In  addition  to  the  above  features,  there  era  s  nuaber  of  other  features 
that  are  required  by  an  emission  nodal  if  Its  outputs  are  to  serve  as  Inputs 
to  other  computer  codes.  Four  types  of  codes  that  are  of  special  Interest 
are: 


s  Display  codes  that  generate  special  foras  of  graphical  and 
tabulated  data. 

e  Statistical  analysis  packages. 

e  Source-oriented  atnospherlc-pollutant  concentration  nodels 
(Including  Ceusslan-pluas  nodels). 

e  Grid  or  cell-type  pollutant  concentration  nodels. 

Special  features  that  asy  be  required  by  these  codes  are: 

1.  Detailed  Individual  source  descriptions,  Including  source 
location  and  geometry.  Such  descriptions  ecu  required  for 
source -oriented  dispersion  models,  such  as  ths  Gausslan- 
plums-type  models,  if  source  effects  on  atmospheric 
dispersion  most  be  taken  Into  account  (see  Sec.  3.2.4). 
For  such  nodels  it  Is  convenient  to  separsto  the  source 
geometries  Into  points,  linns,  and  arses.  It  Is  also 
necessary  to  provide  the  Information  needed  to  compute 
plums  rise.  If  any,  and  Initial  plums  sisa. 

2.  Tims -dependence  of  source  activity  and  omission  rates. 
Hourly  averages  ace  adequate  for  most  purposes. 

3.  Orldded  omission  rates  (aggregated,  for  esanple,  to 

uniform  1  km  1  km  grid  cello).  Thaos  arm  required  foe 
grid  or  coll  concentration  nodels  used,  for  eaampls,  for 
some  photochemical-smog  sins  lotions.  Orldded  emission 
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rates  asy  also  ba  required  for  aoaa  types  of  display 
packages  (e.g.,  cootour-plottlng  packages)* 

Finally,  It  la  leportant  to  recognise  that  emission  rates  for  non- 
routine  pollutants  nay  be  required  for  special  purposes.  A  few  of  these 
purposes  and  their  requirements  are  described  below. 

1.  Calculation  of  MO  and  HOj  concentrations  requires  separate 
ealsslon  rates  for  NO  and  NO2,  In  contrast  to  the  routine 
reporting  of  NO^  ealsslon  factors  (NO  +  NOj). 

2.  Calculation  of  photocheaical-snog  products  requires 
separation  of  hydrocarbon  ealsslons  Into  several  separate 
classes  of  coapounde.  The  particular  classlf lest Ion 
ache as  to  be  used  has  not  been  uniquely  defined  as  yet. 

3.  Calculation  of  visibility  effects  and  odors  requires 
estlaatee  of  ealsslons  of  selected  gas  and  condensed-phase 
pollutants.  The  specific  effluents  of  greatest  concern 
have  not  been  fully  Identified  as  yet. 

4.  Calculations  of  concentrations  of  special  chemical 
aeroeols  of  allltary  or  agricultural  laportance  would  also 
require  appropriate  ealsslon  rates  but  are  best  handled  by 
special-purpose  hand  calculations  rather  than  by  general 
ealsslon  aodele. 


3.2  AIR -QUALITY  OR  POLLUTANT-CONCENTRATION  MODEL  REQUIREMENTS 

The  second  step  In  a  detailed  alr-quallty  analysis,  after  the  compila¬ 
tion  of  an  ealsslon  inventory,  la  coaputatlon  of  pollutant  concentrations  for 
comparison  with  NAAQSa  or  other  measures  of  significance.  Because  of  the 
coaplexlty  of  aost  probleas,  one  or  aore  computer  models  are  usually  required 
to  perform  these  calculations.  Although  the  term  "dispersion  models"  has 
often  been  used  to  refer  to  such  aodele,  it  Is  soaswhat  of  a  misnomer,  since 
other  phenomena  besides  atmospheric  dispersion  mast  often  be  considered  In  the 
calculation  of  ataospherlc-pollutaat  concentrations.  Hence,  a  aore 
appropriate  term  would  be  alr-quallty  or  pollutant-coocentratlon  models.  The 
specific  type  of  alr-quallty  model  required  by  civilian  and  military  users 
depends  largely  on  the  nature  of  the  application,  rather  than  00  the  user  per 
as.  Several  general  features  of  applications  that  determine  the  corresponding 
modeling  requirements  arm  described  below. 


3.2.1  level  of  Petell 

The  main  issue  regarding  the  level  of  do  tell  required  of  a  particular 
application  in  whether  a  screening  or  a  detailed  alr-quallty  modeling 
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calculation  la  required.  For  ragulatory  applications,  It  la  oftan  uaaful  to 
determine  whether  or  not  a  potantlal  probiaa  exists.  Thla  can  ba  dona  through 
application  of  a  conservative  acraanlng  procadura.  If  a  potantlal  probiaa 
doas  exist,  than  a  aora  datallad  analysis  nay  ba  raqulrad. 

Ttaa  aaln  raqulraaant  for  a  acraanlng  procadura,  aalda  froa  eaaa  of  uae, 
la  that  It  ba  conaarvatlva.  That  la.  It  auat  provlda  an  uppar  bound  to  tha 
aagnltuda  of  tha  aalaalon  rataa  and/or  ataoapherlc  concant rat Iona.  Tha  cloaar 
thla  conaarvatlva  uppar  bound  la  to  tha  expected  aaxlnua  valua,  tha  battar. 
However,  caution  la  always  raqulrad  In  applying  acraanlng  tools,  slnca  thay 
aay  ba  ovarly  conaervatlva.  In  which  case  a  naad  for  aora  datallad  analysis 
aay  ba  Indicated  aora  oftan  than  Is  raally  necessary.  On  tha  other  hand, 
because  of  the  Inherent  slapllclty  of  screening  procedures,  thay  aay  fall  to 
treat  thoaa  aspects  of  a  particular  probiaa  that  are  aost  laportant.  The 
applicability  of  screening  procedures  and  tha  Interpretation  of  tha  results 
should  always  be  subjected  to  expert  review. 

Another  Issue  of  laportanca  Is  whether  It  la  neceasary  to  treat  Indi¬ 
vidual  sources  In  detail  or  whether  It  la  satlafactory  to  aggregate  source 
ealsslona  up  to  soae  less-detailed  level.  Generally,  the  level  of  detail 
required  for  the  treataant  of  source  configurations  Is  determined  by  the 
relevant  source-receptor  distances.  If  these  dlstancee  are  relatively  large 
compared  with  the  lntereouree  separation  distances  of  coaparable  sources,  the 
sources  can  be  aggregated.  However,  if  one  la  interested  in  studying  the 
contributions  of  selected  source  types,  or  If  adjacent  sources  have  physically 
different  characteristics,  then  Individual  source  treataant  aay  still  be 
required. 


3.2.2  Teaporal  and  Spatial  Scales 

The  teaporal-  and  spatial-scale  combinations  pertinent  to  four  general 
types  of  appllcatlona  are  displayed  In  Table  1.  Mote  that  tha  NAAQS-related 
applications  Involve  hourly  to  annual  averaging  times.  In  spite  of  the  fact 
that  tha  shortest  averaging  tine  for  MAAQSa  la  one  hour,  It  la  important  to 
raalisa  that  for  sons  substances  atmospheric  concentrations  are  governed  by 
phenomena  having  much  shorter  characteristic  times.  For  example,  WO  interacts 
within  to  form  MOj  in  tha  presence  of  ambient  levels  of  Oy  Hence, 

simulation  time  scales  that  must  ba  treated  in  a  model  aay  differ  froa  the 
required  output  tine  scales  and  depend  largely  on  the  characteristics  of  tha 
air-quality  variables  of  concern. 

Spatial  scales  are  not  specified  in  the  HAAQSa,  but  it  Is  unlikely  that 
the  public  would  routinely  be  closer  them  a  fee  hundred  maters  to  military  or 
civilian  aircraft  activity.  Hence,  the  nser-f laid,  source -dominated  turbu¬ 
lence  sons  is  usually  net  subject  to  the  HAAQSa.  However,  this  tons  aay  be 
subject  to  occupational  health  end  safety  standards.  At  the  other  end  of  tha 
spatial  scales  are  tha  long-range  transport  processes,  including  global-scale 
processes .  It  is  unlikely  that  aircraft  operations  will  contribute  signifi¬ 
cantly  to  ground-level  concent rations  on  scales  greater  than  tens  of 
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Table  l  Temporal-  and  Spatial-Scale  Combinations 
Relevant  to  Air-Quality  Problems* 


Spatial  Scales 

Temporal  Scales 
(averaging  tiae) 

Up  to 

10s  of 
Meters^ 

10s  to 

100s  of 
Metersc 

100s  of 
Meters 
to  10  k^> 

Up  to 

100a  of 

K1 loaeters 

Global 

Scale 

Praction  of  an 
hour  or 
continuous 

1.3,4 

1,3,4 

1,4 

0 

0 

l  Hour 

1,3,4 

1,2, 3,4 

(AVAP?,  AQAM?) 

1,2,4 

(AVAP,  AQAM) 

1 

0 

3-24  Hours 

1.3,4 

1,2, 3,4 

(AVAP?,  AQAM?) 

1,2,4 

(AVAP,  AQAM) 

1 

0 

1  Month  to  1 

Year  or  More 

1,3 

1.2.3 

(AQAM?) 

1,2,4 

(AQAM) 

1 

1 

*Prob lea-type  code; 

0  Not  applicable 

1  Research 

2  Regulatory  —  NAAQS 

3  Regulatory  —  Occupational  Safety  and  Health  Adainlstratlon  (OSHA) 

4  Special  military. 

bThe  region  in  which  aource  characteristics  significantly  affect  or  dominate 
dispersion. 

cThe  transition  region  in  which  aablent  atmosphere  plays  an  increasingly 
laportant  role  in  dispersion  but  in  which  single  events  or  sources  remain 
distinguishable. 

^The  airport  or  air-base  vicinity  (urban  scale)  in  which  sultisource 
contributions  dominate. 


) 
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kilometers  except  in  epeclal  ceeee.  One  epcdel  ceee  of  coneldereble  concern 
over  the  pent  decode  hee  been  the  N0X  contrlbutlone  of  high-flying  aircraft  to 
upper  tropoepheric  end  lower  stratoepherlc  phenomena.  These  applications  are 
beyond  the  acope  of  the  present  report. 

Table  1  also  shows  which  scale  combinations  are  epplicable  to  the  AVAP 
model  and  AQAM.  The  entries  with  question  marks  mean  the  applicability  of  the 
models  la  not  certain. 


3.2.3  Source  Configuration 

Source  configuration  refers  both  to  the  spatlel  distribution  of  the 
sources  of  Interest  relative  to  the  receptors  and  to  the  geometry  of 
individual  sources.  The  spatial  distribution  and  number  of  sources  Is 

important  for  several  reasons.  First,  if  the  number  of  sources  la  large, 
machine-oriented  computational  procedures  are  almost  certainly  required  to 
reduce  the  tedium  Involved  and  to  reduce  the  likelihood  of  making  errors. 
Second,  the  distribution  of  source  types  determines,  to  a  large  extent, 
whether  source-by-source  or  aggregated  source  analysis  la  needed.  In 
addition,  in  cases  where  emiaslons  from  adjacent  sources  tend  to  overlap 
(e.g.,  overlapping  aircraft  plumes),  it  may  be  necessary  to  investigate  the 
consequences  of  that  overlapping  in  detail.  If  the  pollutant  of  interest  is 
relatively  inert,  the  overlapping  can  be  handled  by  a  simple  superposition  of 
contrlbutlone  from  individual  sources.  However,  if  the  emissions  are 

reactive,  then  the  simple  superposition  principle  may  not  apply  and  some 

alternative  to  the  simple  addition  of  single-source  contrlbutlone  may  be 
called  for. 

Sources  are  generally  divided  b if  geometry  into  line  sources,  point 
sources,  and  area  or  volume  sources.  Line  sources  ere  complicated  by  the  feet 
Chat  their  orientation  relative  to  wind  direction  and  ground  level  must  be 
taken  into  account.  In  addition,  if  a  line-source  geometry  la  chosen  to 
represent  a  runway,  then,  since  the  aircraft  operating  on  the  runway  move  at 
nonuniform  speeds,  the  emission  density  along  the  line  source  will  be 

correspondingly  nonuniform. 

In  view  of  the  inportance  of  line-source  geometries  for  simulating 
aircraft  operations,  the  greatest  attention  to  detail  and  accuracy  is  required 
of  line-source  pollution-concentration  models.  However,  in  the  case  of 
routine  air-quality  problame,  it  hee  been  demonstrated  that  for  ground-level 
eoncent rat lone,  airborne  aircraft  amissions  ere  not  very  important.  Hence,  it 
is  not  necessary  to  provide  detailed  model  treatments  of  elevated-aircraft 
lima  sour erne.  (See  See.  10.3.) 

Finally,  aa  noted  earlier,  spatial-scale  and  source  configuration  are 
closely  linked.  As  the  spatial  scale  is  Increased,  leas  attention  to  details 
of  the  source  configuration  is  required.  This  suggests  that,  since  NAAQSs 
generally  apply  only  to  locations  mall  removed  from  aircraft,  oas  can  avoid 
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detailed  treataant*  of  aircraft  pluae  dynaalca  for  such  regulatory  applica¬ 
tion*  t  Thla  la  a  danger oua  aaauaptlon,  alnce  for  ion  caaea  ouch  phenoaana  aa 
initial  pluaa  dispersion,  pluaa  rise,  and  enhanced  vertical  pluaa  diepereion 
can  have  algnlf leant  lapact  on  downwind  concent rat Iona.  It  la  therafore 
laportant  to  carefully  exanlna  the  level  of  treataant  of  pluaa  dynaalca 
required  before  proceeding  with  a  particular  application.  Theoretical 
Investigations,  Including  eeneltlvlty  teata  with  aore  aophlatlcated  node la, 
are  perhape  the  only  aound  way*  to  provide  guidance  on  thla  laaue. 


3.2.4  Source  Effect* 

Application*  requiring  data  on  ataoepherlc  concentration*  near  aourcea 
require  that  aodela  provide  reaaooable  treataant  of  aource  effect*  auch  aa: 

*  building  and  atack  downwaeh  and  wake  effect*. 

e  Aircraft  and  aircraft -engine -gen* rated  turbulent  wake 
effect*. 

*  Pluae  rle*  and  enhanced  vertical  dlaperalon. 

Such  effect*  Influence  both  the  pluaa  trajectory  and  It*  Initial 
dlaperalon,  and  therefore  elgniflcantly  affect  the  near -fie Id  concent ration*. 


3.2.5  Pollutant  Characterlatlca 


The  following  pollutant  characterlatlca  play  a  doalnant  role  In 
deteralnlng  the  aodeling  requlreaanta  for  a  particular  application: 

1.  Pollutant-reactivity  tiaa  acale  relative  to  the  tlaa  acale 
of  the  problea  (e.g.,  hourly  average  coace at rat loo). 

2.  Priaary  vs.  secondary  pollutants. 

3.  Phase  (gas,  liquid  or  solid  aerosol) • 

4.  luoyaacy. 

Pollutant*  that  arc  audited  directly  into  the  at ensphere  frea  aewree* 
aio  called  priaary  pollutants.  Sana  pollutant* ,  each  aa  Cl^  or  nan reactive 
hydrocarbons,  are  relatively  Inert  over  tins  scale*  of  hear a,  uhila  others, 
such  aa  00,  are  relatively  Inert  over  tins  scales  of  day*;  HO,  on  the  other 
hand,  reacts  on  the  tine  scale  of  second*  with  background  on  oua.  In  case*  In 
which  travel  tines  are  chart  coupe red  with  reaction  ties*  and  concent ration 
averaging  tlaa*  are  no  aore  than  one  hoar,  a  number  of  pollutant*,  including 
00,  0(4,  nonraactlve  hydrocarbons  (and  soon  node rate ly  reactive  hydrocarbons ) , 


SO  2 ,  and  NOx  (NO  +  NOj)  can  be  Created  aa  chough  they  were  Inert. 
Consequently,  they  can  be  handled  with  non reactive  pollutant  concentration 
uodels. 

Secondary  pollutanta,  ouch  aa  NO2  (aone  NO2  la  alao  prlnary),  f  lne- 
partlcle  aerosols,  some  photochemical-smog  compounds,  and  ozone  are  produced 
froe  Interact  Iona  between  prlaary  pollutanta  (gaaeo  and  vaporo)  and  background 
ataoepherlc  conatltuenta.  Since  their  foreatlon  time  acales  vary  froe 
seconds,  for  NO2  converalon  froe  NO  In  the  preaence  of  ozone,  to  an  hour  or 
eore,  for  the  foreatlon  of  photocheelcal-aeog  conatltuenta,  and  hours  to  daya 
for  foreatlon  of  sulfate  aeroaola,  the  node ling  requlreeents  depend  quite 
strongly  on  the  particular  pollutant  of  Interest*  Whereas  the  NO  to  NOj 
conversion  process  can  occur  within  an  Individual  Urcraft  pluee,  eany  of  the 
photochealcal-aeog-foreatlon  processes  and  aerosol -forming  processes  occur 
long  after  the  individual  plueas  have  earged  with  each  ocher  and  with 
background  pollutants.  Hence,  photochemical-smog  models  are  generally 
required  to  operate  over  urban  to  lnterurban  spatial  scales,  while  sulfate 
models  involve  long-range  transport  over  hundreds  to  thousands  of  kilomatera. 

The  pollutant's  phase,  that  is,  whether  It  la  a  gas  or  condensed  liquid 
droplet  or  solid  particle,  affects  the  pollutant  transport  In  the  atmosphere 
and  Its  rate  of  deposition  onto  surfaces.  Whereas  particles  with  diameters  In 
excess  of  10  pm  t end  to  be  removed  rather  rapidly  due  to  gravitational 
settling,  those  with  diameters  less  than  1  pm  are  deposited  out  only  very 
slowly.  In  addition,  many  pollutanta  can  undergo  phase  changes  over  periods 
ranging  from  minutes  to  hours,  depending  upon  the  pollutant  species  involved 
and  the  concentrations.  These  phase  changes  affect  not  only  the  subsequent 
dispersion  but  other  properties  that  may  be  of  concern. 

The  buoyancy  of  a  pollutant  plums  released  Into  the  atmosphere  can  be 
positive,  neutral,  or  negative.  Positively  buoyant  pollutant  plumas  rise  In 
the  atmosphere  until  they  reach  thermodynamic  equilibrium  with  the  surrounding 
air.  Thereafter  they  are  dispersed  by  ambient  turbulencm  and  wind.  Aircraft 
exhaust  plumes  are  positively  buoyant,  but  are  mixed  with  substantial  snglne- 
and  aircraft -wake-gene  rated  turbulence.  This  "Internal"  turbulence  dominates 
the  near -field  dispersion  until  It  is  quenched  by  entrainment  of  background- 
dominated  turbulent  sir.  Neutrally  buoyant  plumas  are  produced  by  sources 
that  emit  pollutants  near  ambient  temperature  and  density.  Negatively  buoyaot 
plumas  can  result  from  evaporation  of  cryogenic  liquids,  or  from  highly 
comceet rated  emissions  of  pollutants  having  molecular  weights  greater  than 
that  of  air,  or  from  apodal  conditions  Involving  gas-aerosol  phase  equilibria 
(as  in  tbs  csss  of  ammonia).  Negatively  buoyant  plumas  often  pose  special 
hasards  because  hlgher-than-noraal  concentrations  can  persist  near  ground 
level  for  longer  periods  because  of  the  suppression  of  the  dilation  that  would 
occur  In  neutral  or  positively  buoyant  situations.  Special  attention  to 
souroa  conditions  Is  required  to  determine  whether  or  not  negatively  buoyant 
conditions  apply. 
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Coaaerclal  aircraft  sources  ar«  unlikely  to  yield  negatively  buoyant 
pluaaa.  However,  hypargollc  fuala  are  aommtima a  carried  by  allltary  alrcreft 
and  are  often  atored  and  tranaported  In  allltary  fact lit lee  and  are  aubject  to 
accidental  releaeee.  Treatnent  of  negative  buoyancy  effecta  la  beyond  the 
a co pa  of  the  preaent  study. 


3.2.6  Air-Quality  Varlablea 

The  type  of  alr-quallty  nodal  and  lte  epeclally  required  feeturee 
depend  on  which  air-quality  variable  the  particular  application  calle  for. 
Soar  applications  require  concent ratlone  of  prloary  pollutante  and  othere  nay 
be  concerned  with  secondary  pollutants  or  derived  quantities.  Since 
secondary-pollutant  concent rat Iona  depend  upon  one  or  aore  conversion 
processes,  their  alaulatlon  can  be  quite  difficult.  Derived  quantities 
represent  a  third  class  of  alr-quallty  variables  of  concern  In  soae  applica¬ 
tions.  These  quantities  are  generally  properties  of  prlaary  or  secondary 
pollutants.  Bxaaplee  are  ataospherlc  extinction  or  related  quantities 
(visibility,  obscuration),  zones  of  flaaasbllity  or  detonablllty,  odor  levels, 
etc.  Since  these  properties  do  not  theaaelves  control  the  pollutant  concen¬ 
tration,  they  can  be  dealt  with  by  algorlthaa  appended  onto  pollutant 
concentration  aodels.  Hence,  applications  calling  for  derived  quantities 
require  a  pollutant  concentration  nodal  plus  an  appropriate  algorltha  to  yield 
the  required  property  of  Interest  In  the  application. 


3.2.7  Terrain  Characteristics 


Most  alr-quallty  aodels  work  beet  for  relatively  flat,  unifora 
terrain.  Plow  patterns  In  conplex  natural  terrains,  shorallne  envlronasnts, 
and  around  nan-aade  structures  generally  cannot  be  accurately  slaulated  with 
slaple  pollutant  concentration  aodels.  Special  flow  aodels  are  usually 
required  for  such  applications.  Fortunately,  however,  aoet  civilian  and 
allltary  applications  are  confined  to  relatively  slaple  terrain  conditions. 
Two  exceptions  are  worth  noting:  first,  on  a  saall  spatial  scale,  local  waka 
effecta  due  to  the  source  structures  theaaelves  or  to  nearby  structures  can 
lead  to  significant  changes  in  wind  flow  and  in  concentration  patterns.  For 
applications  requiring  pollutant  concentrations  that  are  likely  to  be 
Influenced  by  wakes  froa  structures,  field  aaaaureaente  or  peraaetrlsationa 
boned  on  physical  aodellng  are  aoet  desirable.  Pollutant  concentration  aodels 
should  include  peraaetrlsationa  to  handle  such  cases.  Second,  on  a  larger 
spatial  scale,  shoreline  environaente  can  exhibit  special  transport  and 
dispersion  conditions.  In  aany  cases  such  locations  are  associated  with 
periodic  recirculation  patterns  resulting  In  pollutant  accuaulatlon  over  a 
several-day  period,  ae  la  the  Los  Angeles  Basin.  In  addition,  aarlne 
Inversions  can  result  in  fairly  shallow  nixed  layers  and  hlgher-thaa-expected 
inland  pollutant  concentrations.  The  recirculation  problae  cannot  be  dealt 
with  using  slopllf led  aodels.  The  effects  of  shallow  alned  layers  can  be 


effectively  handled  using  simplified  models  provided  that  an  adequate 
description  of  the  particular  region  of  Interest  is  available  and  meteoro¬ 
logical  data  sufficient  to  characterize  typical  and  worst-case  conditions  can 
be  obtained*  Careful  accounting  of  the  diurnal  pattern  of  the  shoreline 
meteorology  is  usually  required*  In  some  cases  it  is  possible  to  simulate  the 
shoreline  meteorology  itself*  although  this  type  of  simulation  is  generally 
beyond  the  scope  of  the  usual  pollutant  concentration  models* 


3*3  METEOROLOGICAL  REQUIREMENTS 

Computation  of  atmospheric  concentrations  requires  meteorological  data 
either  on  a  short-term  (generally  hourly)  or  on  a  long-term  statistical 
basis*  Typical  short-term  data  required  by  source-oriented  models  Include 
wind  speed*  wind  direction,  depth  of  the  mixed  layer*  and  some  measure  of  the 
turbulent  properties  of  ^the  atmosphere*  such  as  the  Pasquill-Cif ford-Turner 
(PGT)  stability  class*  For  the  simpler  models*  one  representative 

meteorological  monitoring  station  is  required*  This  assumes  that  the  region 
of  model  application  is  fairly  uniform  with  respect  to  terrain  features  and 
surface  roughness* 

Several  points  are  worth  noting*  First*  even  in  cases  where  only 
moderate  variations  in  terrain  exist*  the  near-surface  wind  field  can  depart 
significantly  from  spatial  uniformity*  Second*  the  original  PGT  atmospheric- 
stability  classification  scheme  has  been  under  review  for  a  number  of  years* 
and  the  American  Meteorological  Society  (AMS)*2*  in  particular*  has 
recommended  some  alternative  schemes  that  require  additional  meteorological 
data*  In  cases  where  the  additional  data  are  not  available*  calculatlonal 
procedures  based  on  theoretical  considerations  and  meteorological  field 
experiments  have  been  proposed*  This  second  alternative  is*  unfortunately* 
not  guaranteed  to  produce  more  satisfactory  results  than  the  original 
scheme*  (This  latter  point  has  not  been  adequately  addressed  in  the 
literature*) 

Third*  estimates  of  the  depth  of  the  mixed  layer  are  often  not  readily 
available  to  model  users*  Provided  that  the  transport  distances*  or  spatial 
scale*  of  the  application  are  not  more  thao  a  few  kilometers,  the  depth  of  the 
mixed  layer  is  not  too  critical  unless  it  is  relatively  small*  say  a  few 
hundred  meters  or  less*  Unfortunately*  such  low  mixing  depths  can  occur  near 
shoreline  environments  and  can  often  occur  in  other  environments  for  periods 
of  ooe  or  more  hours  in  the  morning  and  sometimes  for  more  extended  periods 
during  the  evening  and  nighttime  hours  under  the  influence  of  urban  heat- 
island  effects.  Consequently*  for  worst-case  calculations*  it  is  important  to 
obtain  at  least  reasonably  good  representative  estimates  of  the  mixed-layer 
depth  as  a  function  of  time  of  day* 


4  OPERATIONAL  REQUIREMENTS  FOR  MODELS 


Operational  requirements  placed  on  computational  procedures  for 
estimating  emissions  and  atmospheric  concentrations  Include  usability v 
compatibility  with  available  computer  hardware,  and  documentation.  Of  these 
three,  usability  has  come  to  be  recognised  as  a  requirement  of  paramount 
importance*  Users,  of  course,  prefer  to  work  with  simple  models  that  do  not 
require  extensive  familiarity  with  technical  detail  or  manpower -intensive  data 
compilations.  At  the  same  time,  it  is  recognised  by  users  that  the  nature  and 
technical  characteristics  of  air-quality-related  problems  cover  a  fairly  broad 
spectrum.  Hence,  flexibility  is  also  an  important  requirement  of  a  model. 

The  rapid  growth  in  the  use  of  desk-top  terminals  together  with 
minicomputers  or  microcomputers  has  led  to  the  need  for  highly  efficient 
computer  software  that  requires  little  central  processing  unit  (CPU)  time  and 
core  storage.  This  requirement  Is  generally  Inconsistent  with  large,  complex, 
flexible  computer  codes  unless  those  codes  can  be  redesigned  in  a  modular 
fonsat. 


Another  aspect  of  usability  is  the  choice  of  Interactive-  or  batch-mode 
operation.  In  the  Interactive  mode,  the  user  sits  at  a  computer  terminal  and 
responds  to  a  series  of  questions  or  prompts  Issued  by  a  user-interactive 
computer  program.  The  user's  responses  provide  the  inputs  needed  by  the 
program  to  perform  the  computation.  In  the  batch  mode,  the  user  prepares  a 
deck  of  cards  or  a  computer  file  containing  all  the  necessary  Input  data. 
This  preparation  Is  generally  done  by  following  a  series  of  directions 
provided  by  a  user's  manual.  Once  the  input  file  is  completed,  the  user 
submits  a  job  that  requires  the  computer  program  to  read  In  the  input 
information  and  perform  the  calculation.  The  process  of  submitting  a  Job 
Involves  reading  the  deck  of  cards  into  the  computer  or  keying  in  a  brief  set 
of  instructions  via  a  computer  terminal. 

From  the  point  of  view  of  the  user,  the  Interactive  mode  is  often 
preferred.  This  Is  especially  true  for  users  having  little  or  no  familiarity 
with  the  computer  program.  Beginning  users  tend  to  be  intimidated  by  user's 
manuals  that  require  detailed  study  prior  to  using  the  computer  program.  In 
the  interactive  mode,  the  uninitiated  user  can  often  begin  "playing"  with  the 
computer  program  with  little  or  no  preliminary  study  of  documents.  However, 
for  serious  computations  the  differences  in  effort  required  nay  turn  out  to  be 
quite  esmll. 

The  disadvantages  of  the  Interactive  mode  are  that  the  user-interactive 
program  must,  of  necessity,  contain  substantially  more  programming  instruc¬ 
tions  than  a  program  written  for  batch-mode  operation.  This  naans  more  core 
storage  and  longer  running  time.  In  addition,  in  order  to  remain  within  the 
limitations  of  particular  computer  installations,  it  may  be  necessary  to  limit 
the  number  of  options  available  to  the  user  and  therefore  limit  the 
flexibility  and  degree  of  sophistication  of  the  main  computational  parts  of 
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the  computer  program.  Furthermore,  for  large  problems,  it  may  prove  tiresome 
to  remain  at  a  computer  terminal  for  long  sessions  using  the  Interactive  mode 
and  preparing  long  streams  of  input  Information. 

One  way  to  avoid  long  interactive  terminal  sessions  Is  to  use  Input- 
file  editors.  These  are  programs  that  permit  the  user  to  modify  existing 
large  input  data  sets  without  having  to  redevelop  the  data  sets  from  scratch. 

The  problems  of  Halted  core  storage  and  run  time*  as  well  as  the 
problem  of  limited  options  or  flexibility,  can  be  overcome  by  designing  mixed¬ 
mode  computer  programs  that  can  be  executed  In  parts.  For  example.  It  may  be 
desirable  to  divide  a  computer  program  Into  two  or  more  components.  The  first 
component,  which  takes  in  and  organizes  the  input  data,  could  be  operated  In 
the  interactive  mode,  and  the  second  part,  which  contains  the  main  computa¬ 
tional  algorithms,  could  be  designed  for  operation  in  the  batch  mode.  In 
addition,  for  Increased  convenience 
the  first  part  could  be  supplemented 
with  an  editor  program  that  could  be 
used  to  make  modifications  to  existing 
large  sets  of  Input  data.  The  overall 
structure  of  such  a  package  of 
computer  program  components  would  look 
something  like  the  macro  flowchart 
shown  in  Fig.  1.  Note  that  in  this 
flowchart  an  additional  component  for 
creating  various  output  displays  is 
also  Illustrated.  For  increased 
flexibility  or  to  accommodate  varying 
technical  requirements  of  different 
applications,  several  alternative 
computational  program  components  could 
be  substituted,  provided  thst  their 
input  and  output  structures  were 
compatible  with  the  other  program 
components. 

A  final  operational  requirement 
for  s  computational  system  Is  good 
documentation.  Generally,  the 

documentation  should  consist  of  the 
following  parts: 

1.  A  description  of  the 
Intended  uses  of  the 
program. 


Fig.  1  Macro  Flowchart  of 
a  Mixed -Mode  Computer 
Program 
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2.  A  technical  description  of  the  equations  end  algorithms 
used  and  a  tabulation  of  the  important  parameters  and 
constants* 

3.  A  description  of  the  cosputer  programs.  Including  a  macro 
flowchart  and  descriptions  of  all  subroutines;  possibly 
also  detailed  flowcharts  and/or  program  listings. 

A.  A  detailed  user's  manual  that  Includes  guidance  in 
selection  of  various  options,  if  there  are  any. 

S.  An  example  problem  that  can  serve  as  a  guide  to  the  user 
as  well  as  a  benchmark  for  comparison  with  results 
obtained  by  the  user  on  hie  own  computer  facilities. 
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b  SUMMARY  ANT>  FVA1  .HATTON  OF  THE  AVAP  COMpOTFR  MOOFL 


The  two  main  models  developed  sped fi cal lv  for  wllltarv  and  civilian 
air-quality  problems  involving  aircraft  are  the  AOAM  and  the  AVAP  model, 
reapecti vel v.  These  models  were  originally  developed  in  the  early  1970*  and 
have  subsequent  ly  undergone  various  modi  Mentions,  ref  inements ,  and  updates. 
Tn  Secs.  5  and  b,  the  AVAP  model  and  AOAM*  along  with  their  various  modified 
versions,  will  he  discussed  in  detail.  Their  intended  uses  and  strengths,  as 
well  as  their  limitations  and  weaknesses,  will  he  described.  Comparisons  will 
he  made  In  Sec.  7,  and  the  best  features  of  each  model  will  hr  indicated  for 
possible  future  Incorporation  into  a  proposed  lolnt  cl vi 11  an /mi  1 i tary  modeling 
package  (see  Sec.  9). 


S.1  THF,  AVAP  MOOFL 

The  AVAP  model  was  originally  intended  for  evaluating  the  air-qualitv 
Impacts  of  large  civilian  airports.  The  overall  structure  of  the  model  is 
illustrated  in  the  macro  flowchart  shown  in  Fig.  2.  All  of  the  components 
Illustrated  In  this  chart  are  contained  In  a  single  computer  program  composed 
of  a  number  of  subroutines  and  a  main  driver.  Both  the  emission-inventory  and 
dispersion  algorithms  are  contained  in  this  package  and  are  executed  together 
in  the  Hatch  mode.  The  outputs  can  he  generated  in  printed  formats  or  punched 
cards  or  stored  on  magnetic  tape  for  future  use.  The  user  can  choose  any 
conf iguratlon  of  airport  he  wishes,  including  several  runways  and  associated 
taxiwavs,  gate  areas,  etc.  In  addition  he  can  define  aircraft,  airport- 
nonaircraft,  and  environ  sources,  and  can  designate  anv  of  these  as  point, 
line,  or  area  sources.  The  line  sources  can  have  arbitrary  lengths  and 
orientations  in  three  dimensions.  The  model  is  designed  to  compute  emissions 
and  atmospheric  concentrations  of  relatively  inert  pollutants  on  an  bour-bv- 
hour  basis  in  24-hr  cvcles.  The  model  was  not  designed  to  compute  long-term 
averages.  However,  long-term  Averages  can  he  computed  hv  combining  the 
hour tv-average  values  either  in  sequence  or  in  some  statistical  fashion* 

A  uniform  receptor  grid  and/or  a  set  of  specially  located  receptors  can 
he  defined  hv  the  user.  To  allow  the  model  to  he  applied  to  an  arbitrary 
runway  configuration,  th«  user  la  asked  to  specify  runway  operations  in  terms 
of  four  90°  wind  quadrants.  These  quadrants  can  he  arbitrarily  oriented,  hut 
the  user  must  define  the  wind-quadrant  orientation  angle  that  measures  the 
orientation  relative  to  geographic  north.  (This  angle  should  he  between  0° 
and  90*.)  Then,  for  each  of  these  four  wind  quadrants,  the  user  must 
allocate,  for  each  class  of  aircraft,  the  runways  that  are  available  for 
arrivals  or  departures.  This  user-sunplted  Information  allows  the  computer 
model  to  determine  which  runways  are  being  used  during  each  hourly  set  of 
meteorological  conditions  for  each  aircraft  class. 


Overall  Structure  of  the  Airport  Vicinity  Air  Pollution  Model 
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The  simulation  of  both  airborne  and  ground-based  aircraft  operations  Is 
based  on  an  aircraft  classification  scheme  that  Is  described  In  the  user's 
manual.  This  scheme  can  be  altered  by  the  user,  provided  the  alterations  are 
done  in  a  self-consistent  manner*  The  classification  scheme  is  described  in 
the  next  several  paragraphs* 

Each  aircraft  type,  such  as  a  B727  or  DCIO,  Is  assigned  a  "Type  Index" 
value*  The  arrival  rate  Is  given  by  the  user  In  terms  of  the  number  of 
arrivals  of  each  aircraft  type  per  hour  to  the  airport*  Service-vehicle  types 
and  their  operating  times  are  also  given  in  terms  of  the  aircraft  "Type 
Index."  To  each  aircraft  type  Is  assigned  four  additional  Indices:  Aircraft 
Class,  Aircraft  Range,  Number  of  Engines,  and  Engine  Type* 

The  Aircraft  Class  Index  is  used  by  the  computer  model,  together  with 
the  hourly  wind  direction,  to  assign  takeoff  and  landing  runways  and  outbound 
and  Inbound  taxlways*  Runway  availability  Is  specified,  as  already  mentioned, 
by  the  user  input  data  in  terms  of  four  wind  quadrants  and  aircraft  class* 
(Taxiway  availability  Is  determined  by  runway  and  airline  or  gate  position.) 

The  Aircraft  Range  Index  Is  used  to  assign  certain  aircraft  operational 
parameters,  Including  taxi  speed  and  tlmes-ln-mode*  The  numbers  and  types  of 
engines  are  used  to  assign  pollutant  emission  factors.  An  example  of  the 
classification^  scheme,  which  was  used  in  the  example  problem  In  the  AVAP 
user's  manual.  Is  shown  In  Table  2* 

The  airborne  operations  simulated  In  the  model  Include  one  approach 
path  (defined  by  the  user)  and  one  or  two  departure  paths  (defined  by  the 
user)  for  each  runway*  These  airborne  optional  paths  are  not  Aircraft-Class 
dependent  except  to  the  extent  that  the  runways  themselves  are  Aircraft-Class 
dependent* 

One  diurnal  aircraft  arrival  pattern  Is  provided  by  the  user  for  each 
aircraft  type*  To  simplify  the  amount  of  input  data  required,  the  model 
applies  the  same  pattern  to  all  airlines  employing  a  given  aircraft  type* 


5.2  EXISTING  VERSIONS  OP  AVAP 

The  original  version  of  the  AVAP  model  la  from  circa  1973*  It  has  been 
published  by  National  Technical  Information  Service  (NT1S),  la  available  on 
magnetic  tape,  and  has  bae£  well  documented  wit^  both  technical  reports,6* 
validation  study  reports,  *  and  a  user's  manual.  The  user's  manual  contains 
benchmark  runs  and  an  example  application  to  Washington  National  Airport  (OCA) 
to  help  aid  the  user*  The  source-emission  factors  and  source-activity 
algorithms  were  developed  from  published  data  and  airport  observations  dating 
from  1969  to  1973*  The  aircraft  types  treated  In  the  model  have  not  changed 
appreciably  from  the  NTIS  publication  (which  appeared  In  197S-1976)  to  present 
times* 


Table  2  Aircraft  Clasaif 1 cat ton  Scheae  Uaed  in  tha 
Bxeaple  Problaa  In  tha  AVAP  Uaar'a  Guide 


Aircraft 

Aircraft 

Aircraft 

Aircraft 

Bnglnaa  par 

Englna 

Engine 

Typo 

Modal 

Clean 

Range 

Aircraft 

Type 

Modal* 

1 

Boa inf  727 

4 

3 

3 

3 

JT80 

2 

Douglaa  DC9 

4 

3 

3 

3 

JT80 

3 

Booing  737 

4 

3 

3 

4 

JT8D 

4 

Convalr  380 

5 

3 

2 

5 

A-501-D13 

5 

1AC  111 

7 

3 

2 

6 

Spey-511 

6 

NAHOO  YS1 1 

5 

3 

2 

5 

A-S01-D13 

7 

Beach  99 

6 

A 

2 

A 

TPE331 

8 

Fairchild 

rM-227 

3 

3 

2 

5 

A-S01-331 

9 

TWln  Otter 

6 

4 

2 

A 

TPR331 

10 

Piatoe 

Baglne 

9 

A 

2 

7 

d-320 

^Aircraft  angina  node  la  for  which  no  oaiaalona  data  are  available  are  rnplacad 
b f  nodala  having  ainilar  charactnriatlca. 


A  alnpllfied  vnralon  of  AVAP  that  will  bn  referred  to  na  the 
“abbreviated  vereion”  (circa  197})  waa  developed  aa  an  alternative  to  tha  full 
AVAP  aodel  for  acreenlng  applicattoae.  Mil la  tha  dieperaion  algorithea 
renal end  virtually  unchanged,  tha  Input  a true Cure  waa  coaaidarably  nod if lad, 
and  far  fewer  opt  Iona  were  nede  available  to  the  wear.  The  abbreviated 
vereion  la  Much  eaaler  to  uae,  niece  Moot  of  the  deciaioea  have  been 
“internalised. "  Only  one  runway-taxlway-teradaal  area  configuration  la  avail** 
able,  and  no  eourcaa  outalde  the  airport  (environ  eoureee)  ate  pornitted.  The 
aedel  caeputar  code  la  available  through  the  PAA  and  haa  bean  docunantod  with 
a  detailed  but  unpubllahod  uaer'a  annual.* 

An  updated  vereion  of  the  original,  or  VTlt-publiahed,  vereion  of  the 
AVAP  nadel  waa  prepared  and  applied  to  the  aoeeoanaat  of  air-duality  inpecta 
at  four  aajor  O.S.  airport a  around  19t0*  Thin  vereion,  which  will  be  referred 
to  aa  the  “epdatad  vereion,"  haa  net  been  fernally  pebllahed  or  estenalvely 
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documented.  However,  e  description  of  the  ndel,  eloog  with  the  results  of 
the  air-quality  assessment  mentioned  above,  iias  been  reported  in  three 
documents.  *  The  updated  version  contains  updates  to  both  the  emission  and 
dispersion  algorithms.  The  updates  are  based  on  new  emission-factor  data 
published  by  the  Environmental  Protection  Agency  (EPA),  observations  of 
aircraft  operations  at  several  major  airports,  and  field  programs  involving 
meteorological  and  air-quality  observations  at  Dulles  International  Airport 
and  Washington  National  Airport. 


5.3  COMPARISONS  AMONG  AVAP  VERSIONS 

Por  convenience,  the  three  versions  of  AVAP  mentioned  above  will  be 
designated  as  follows: 

Version  1:  Original  NT1S  version,  circa  1973 

Version  2:  Abbreviated  version,  circa  1975 

Version  3:  Updated  version,  circa  1980 

Comparisons  can  be  conveniently  made  in  the  categories  of  overall 
structure,  input  structure  and  user  options,  emission  parameters,  dispersion 
algorithms,  output,  and  documentation. 


5.3.1  Overall  Structure 


The  three  versions  all  contain  both  the  emission  and  dispersion 
algorithms  and  are  operated  In  the  batch  mode. 


5.3.2  Input  Structure  and  User  Options 

The  Input  structures  ot  the  three  versions  are  different.  Versions  1 
and  3  are  quite  similar  except  for  one  Important  difference.  Whereas  Version 
1  requires  only  the  arrival  rates  as  input,  Version  3  also  requires  the 
departure  rates  as  input.  Version  l  incorporates  an  empirical  algorithm  that 
takes  into  account  aircraft  turnaround  time  sod  overnight  parking  to  calculate 
departure  rates.  This  procedure  was  followed  because,  at  the  time  of  program 
development,  only  arrival  times  of  aircraft  were  readily  available  from  the 
commercial  airline  guide.  because  of  the  Influence  of  overnight  aircraft 
parking.  Version  1  performs  calculations  In  24-hr  cycles  (as  many  as  the  user 
chooses).  Version  3,  on  the  other  hand,  makes  use  of  both  user-supplied 
arrival  data  and  departure  data  available  from  more  recent  Issues  of  the 
commercial  airline  guide.  Hence,  Version  3  can  be  operated  on  an  hour-by-hour 
basis.  Also,  Version  3  contains  updated  emission  factors  and  updated  tlmae- 
1 a -mods  for  aircraft  sources.  These  updates  do  not,  however,  const Itwte 
changes  In  input  structure,  ns rely  changes  In  input  parameter  values. 
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Version  2  requires  only  arrival  rates  as  Inputs t  but  it  assumes  that 
the  number  of  arrivals  during  a  given  hour  equals  the  number  of  departures. 
Version  2  performs  only  one  one-hour  computation  at  a  time.  Version  2  also 
differs  from  Version  I  In  a  number  of  other  significant  ways.  Whereas  Version 
1  permits  the  user  to  specify  a  great  variety  of  aircraft*  airport- 
nonaircraft,  and  environ  sources  with  point,  line,  and  area  geometries,  the 
choices  permitted  by  Version  2  are  very  limited.  In  Version  2  the  user  can 
use  at  most  one  runway,  one  inbound  apron,  one  Inbound  taxiway,  one  outbound 
taxiway,  one  outbound  apron,  up  to  60  access-vehicle  straight-line  roadway 
segments,  one  aircraft  area  source,  and  one  nonaircraft  area  source. 

The  aircraft  area  source  includes  emissions  from  four  source  types: 

1.  Aircraft  ground -service  vehicles. 

2.  Auxiliary  power  units  (APU)  on  aircraft. 

3.  Aircraft  taxiing  within  the  area. 

4.  Aircraft  engine  idling  within  the  area. 

In  addition,  whereas  In  Version  l  the  user  must  supply  moat  airport  and 
aircraft  source  parameters,  this  burden  is  largely  removed  through  the  uee  of 
default  values  provided  in  Version  2.  The  ueer  has  the  option  of  overriding 
up  to  23  defaults  simply  by  entering  l's  In  the  appropriate  columns  of  a 
default  control  card  and  then  adding  the  necessary  input  cards  containing  hla 
own  values. 


5.3.3  Emission  Parameters 


The  emission  parameters  for  nonaircraft  airport  and  environ  source 
types  are  Identical  for  all  three  versions.  (Version  2  does  not  treat  environ 
sources.)  The  aircraft  emission  factors  are  the  same  In  Versions  1  and  2  but 
have  been  updated  in  Version  3  to  1977  values. 


5.3.4  Dispersion  Psrsmstsrs  and  Algorithms 

The  dispersion  psrsmstsrs  and  algorithms  In  Versions  1  and  2  are 
Identical.  There  are  several  differences  In  Version  3.  These  are: 

1.  An  empirical  plume -rise  equation  baaed  on  the  Dulles  3- 
Tower  field  Experiments  was  Introduced  into  the  taxi-idle 
aircraft  mode. 


2.  New  Initial  plums  dimensions  baaed  on  the  Dulles  and 
Washington  National  axperl manta  were  used. 
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3*  Initial  plume  dimensions  were  incorporated  into  the 
runway,  taxiway,  and  apron  line  sources  using  the  addition 
of  variance  in  lieu  of  the  virtue  1-point -source  method. 

All  other  source  types  used  the  vlrtual-polnt-source 
method. 

4.  Enhanced  vertical  dispersion  due  to  plume  rise  was 
Incorporated  into  the  treatment  of  all  sources. 

5.  Only  downwind-distance-dependent  vertical-dispersion  pa¬ 
rameters  were  used,  in  lieu  of  the  maximum  of  the  above 
parameters  and  the  time-dependent  vertical-dispersion 
parameters. 

5.3.5  Output 

The  model  outputs  were  changed  in  Version  3  to  provide  additional 
emission  and  concentration  information.  In  addition,  the  output  was  modified 
to  Insure  compatibility  with  a  contour-plotting  package. 


5.3.6  Documentation 

Version  1  has  technical  documentat ion,  model  validation-study 
documentation,  and  a  user's  manual,  all  published  by  FAA.  Version  2  has  only 
an  unpublished  user's  manual.  Version  3  has  supporting  technical  documenta¬ 
tion  but  no  formal  model  documentation  or  user's  manual. 


5.4  AVAP  MODEL  EVALUATION 


5.4.1  Intended  Uses  and  Strengths 

Version  1  (original  NTIS  version)  was  intended  to  be  used  to  assist  in 
(1)  preparation  of  EISs  for  new  construction,  (2)  modification  of  existing 
airports,  and  (3)  evaluation  of  alternative  strategies  for  emission  control. 
Examples  of  changes  in  airports  that  could  be  analysed  Include  fuel-use 
changes,  changes  In  the  physical  airport  layout  (additions  of  taxlways, 

runways,  access-vehicle  roadways,  etc.),  growth  In  operations,  and  environ 
source  changes.  In  Version  1,  the  user  la  provided  with  a  great  deal  of 

flexibility  regarding  his  choice  of  airport  operating  parameters,  emission 
factors,  and  treatment  of  source  types  in  terms  of  points,  lines,  or  areas. 
Listings  of  parameter  values  are  provided  In  the  user's  manual,  but  the  user 

can  use  alternative  values  If  he  feels  they  are  mors  suitable  for  his 

particular  application.  However,  regardless  of  whether  or  not  he  uses  the 
values  listed  in  the  user's  manual,  the  values  he  uses  must  be  entered  Into 
the  computer. 


tin 


I 
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In  Version  2  (abbreviated  version) ,  most  of  the  parameter  values  are 
provided  as  automatic  defaults  that  do  not  have  to  be  entered  by  the  user.  If 
the  user  wants  to  override  these  default  values,  then  he  must  provide 
additional  input  information.  Version  2,  however,  is  Intended  only  as  a 
screening  tool  to  obtain  a  rough  idea  of  the  air-quality  impacts  that  might 
occur  if  all  airport  operations  were  confined  to  a  simple  configuration 
consisting  of  only  1  runway  and  a  few  other  distinct  sources.  Hence,  Version 
2  does  not  provide  nearly  as  much  flexibility  to  the  user  as  Version  I  does. 

All  three  versions  of  AVAP  permit  the  user  to  readily  update  emission 
parameters.  Many  of  these  parameters  require  periodic  updating  or  adjustments 
to  account  for  faclllty-to-faclllty  variations.  All  three  versions  are 
designed  for  operation  in  batch  mode  and  combine  emission  and  dispersion 
calculations  in  a  single  computer  code.  Consequently,  while  they  are  not 
"ueer-f rlendly"  in  the  sense  of  an  Interactive  code,  they  are  more  efficient 
to  operate,  especially  on  large  problems,  and  do  not  require  the  long  terminal 
sessions  characteristic  of  large  interactive  codes. 


The  Version  1  dispersion  algorithms  are  designed  to  provide  a  somewhat 
conservative  estimate  of  relatively  inert  primary-pollutant  concentrations. 
Pollutants  such  as  S02 ,  CO,  NO^,  HC  or  MfHC,  and  total  suspended  particulate 
(TSP)  can  be  considered.  No  chemical  or  physical  transformations  are 
considered,  and  no  removal  processes,  such  as  wet  or  dry  deposition  or 
gravitational  settling,  are  accounted  for. 


The  dispersion  parameters  correspond  to  greater  dispersion  rates  than 
those  used  in  Turner's  Workbook  of  Atmospheric  Dispersion  Estimates *  because 
they  are  adjusted  for  one-hour  sampling  times  and  urban  dispersion 
conditions.  Comparisons  indicate  that  they  more  closely  represent  field 
observations.  The  absence  in  Versions  1  and  2  of  any  treatment  of  plume  rise 
or  enhanced  vertical  dispersion  from  aircraft  sources  leads  to  overestimates 
of  concentrations  near  such  sources.  This  problem  is  partially  corrected  in 
Version  3,  which  does  account  for  these  effects  for  slow-moving  aircraft.  The 


treatments  of  these  effects  in  Version  3  is  based  on  the  Dulles  3-Tower  field 
12 

Experiments. 


The  AVAP  model  Incorporates  simple  treatments  of  building  and  stack 
downwaeh  and  wake  effects  and  aircraft  wake  and  Jet  turbulence  effects. 
Briggs9  "rules -of -thumb"  are  used  to  indicate  when  the  building  and  stack 
effects  are  likely  to  occur,  and  such  effects  are  expressed  In  terms  of 
reduced  plume  rise  and/or  enhanced  initial  plume  dispersion,  finite  initial 
plume  sixes  are  handled  with  the  virtual-point -source  method.  Aircraft  jet- 
wake  and  turbulence  effects  are  similarly  treated  by  enhanced  initial  plume 
dispersion.  The  parameters  used  for  this  pumas  in  Version  3  are  baaed  on 
the  Dulles  and  DGA  Airport  field  Experiments.  Ho  attempt  is  made  to  treat 
these  complex  phenomena  in  detail  la  the  near  field.  The  model  applies  only 
after  the  plumes  become  passive  objects  Influenced  by  ehbient  atmospheric 
turbulence.  Hence,  the  model  should  not  be  used  to  estimate  concent  rations 
too  close  to  large  structures  or  aircraft. 
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Complex  terrain  effects  ere  Ignored  In  AVAP.  The  use  of  this  model  In 
valley  situations  or  near  shorelines  is  not  advised  unless  the  effects  of 
these  features  are  veil  known*  The  models  are  applicable  to  nearly  flat 
regions  on  the  order  of  10*20  ka  on  a  side* 

Since  all  versions  aasuae  steady-stats  conditions,  that  is,  constant 
emission  and  dispersion  rates  for  at  least  one  hour,  transient  phenomena 
cannot  be  accounted  for*  Repeated  transient  phenomena,  such  as  aircraft 
takeoffs,  are  treated  as  though  their  air-quality  lapacts  could  be  averaged 
over  one  hour.  This  treatment  Is  satisfactory  only  for  relatively  Inert 
pollutants  and  averaging  tlaea  on  the  order  of  one  hour.  Such  a  procedure 
must  be  aodlfled  for  situations  where  shorter  averaging  tlaes  are  required. 
Also,  fairly  reactive  ^  and/or  secondary  pollutants  or  highly  variable 
Meteorological  conditions  require  aore  detailed  treataents. 


5.4.2  Limitations  and  Weaknesses 


Llaltatlons  and  weaknesses  are  not  always  perceived  as  such  by  various 
users*  Soaetiaee  a  limitation  aay  be  regarded  as  a  blessing.  Nevertheless, 
an  attempt  has  been  made  to  list  what  way  be  regarded  by  some,  at  least,  as 
shortcomings  of  the  various  AVAP  versions.  These  are  discussed  under  the 
three  separate  headings  of  Usability  and  Availability ,  Documentation,  and 
Technical  Issues. 


Usability  and  Availability 

Version  i  tends  to  be  hard  to  use  because  Its  application  requires  a 
careful  review  of  the  user's  manual.  In  addition,  for  large  problems,  a 
considerable  amount  of  Input  data  must  be  prepared  by  the  user.  Version  2, 
however,  provides  many  default  values  and  alleviates  much  of  this  problem. 
However,  this  reduction  In  Input  preparation  time  Is  achieved  at  the  expense 
of  flexibility  and  accuracy  of  representation  of  actual  airports.  Of  the 
three  versions,  only  Version  1  Is  readily  available  through  NTIS. 

Because  the  emission  aod  dispersion  algorithms  are  combined  into  a 
single  package,  and  because  of  the  options  regarding  source  Inputs,  the  AVAP 
codes  all  require  rather  large  computer  core  storage*  In  addition,  the 
algorithms  are  not  optimally  written  and  therefore  can  consume  significant 
amounts  of  computer  time,  especially  when  large  numbers  of  repetitive 
calculations  are  involved*  Hence,  these  codes  require  fast,  large,  mainframe 
computers*  Also,  it  Is  Inefficient  to  use  them  on  small  problems  Involving 
only  a  few  sources* 


Documentation 

The  documentation  of  Version  I  was  quite  adequate  from  both  the  user 
and  technical  points  of  view  at  the  time  it  was  prepared*  However,  neap  of 
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the  emission  parameters  either  listed  In  the  user's  manual  or  Incorporated  as 
defaults  Into  the  emission  model  are  outdated.  In  addition,  technical  changes 
have  been  incorporated  Into  Version  3  that  have  not  been  adequately 
documented*  Further  improvements  have  also  been  planned  but  have  not  yet  been 
thoroughly  documented* 


Technical  Issues 

No  treatment  of  reactive  or  secondary  pollutants  (e.g*,  reactive 
hydrocarbons,  NO2.  O3,  or  secondary  fine  particulate  matter)  Is  available*  NO 
to  N02  conversion  Is  not  treated*  Aerosol  formation  is  not  treated*  There  is 
no  treatment  of  complex  terrain  features*  Applications  are  limited  to  an 
airport  and  its  near  vicinity.  Only  neutral  and  positive  buoyancy  are 
treated*  Treatment  of  individual  aircraft  plume  dynamics  is  not  adequate  tor 
research  studies  of  near-field  phenomena  (although  such  effects  may  not  be 
relevant  for  ElS-type  applications).  Building -wake  effects  are  not  adequately 
treated  to  provide  reliable  estimates  of  concentrations  near  terminal 
buildings*  Nonaircraft  emission  factors  are  out  of  date*  Plume  rise  In 
takeoff  mode  is  not  treated.  Only  Carson-Moses  and  Holland  plume-rise 
formulas  (see,  for  example,  Ref*  15)  are  currently  available  as  options  for 
nonaircraft  sources*  No  treatment  of  calm  conditions  other  than  persistence 
of  preceding  or  following  nonzero  wind  conditions  Is  Included*  The  present 
line-source  algorithm  Is  Inefficient  and  has  been  shown  to  be  Inaccurate  under 
certain  special  input  conditions*  The  dispersion  parameters  are  based  on 
earlier  dispersion  modeling  studies  and  the  PGT  stability  classification 
scheme*  The  recent  recommendations  of  the  AMS  should  be  evaluated  and 
considered  for  Incorporation  Into  the  AVAP  models*  Currently,  the  user  oust 
either  supply  a  value  for  the  depth  of  the  mixed  layer  or  choose  the  default 
procedure,  which  Is  based  on  the  Holzworth  climatological  value  for  the 
region*  Some  alternative,  more  adequate  default  procedure  should  be 

developed*  Specific  jet-engine  thrust  settings  are  assumed  to  correspond  to 
aircraft  operational  modes*  These  My  require  updating. 
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dates  Essential  for  Inert-Pollutant  and/or  Screening 
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Updated,  faster,  more  accurate  line-  and  point-source  dispersion 
algorithms  should  be  Incorporated*  Updated  documentation  to  Include  the 
latest  improvements  In  algorithms  and  various  parameter  values  should  be 
added*  Dispersion  parameters  and  the  turbulence  classification  scheme  should 
be  updated*  Updating  of  the  aircraft  plume-rise  algorithm  and  lnltlal- 
dlsperslon  parameters  for  all  aircraft  modes,  especially  takeoff.  Is  needed. 
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5*4.4  Additional  Updates  PiilttbU  for  lot rt -Pollutant 
and/or  Scrmlng  Applications 

If  nonaircraft  aourcaa  art  to  be  considered*  emission  factors  should  be 
updated.  If  cala  conditions  are  of  concern*  a  calm  algorithm  should  be 
incorporated.  The  original  OTIS  version  should  be  replaced  with  Version  3  and 
additional  updates  undertaken  as  indicated  above. 


5.4.5  Updates  Weeded  for  Reactive-Pollutant  Applications 

Development  and  testing  of  algorithms  to  treat  ilO/NOj  conversion  In 
aircrsft  plumes  should  be  undertaken  (see*  for  example*  kef.  17).  Development 
and  testing  of  algorithms  that  can  be  used  for  treating  more  general  reactive- 
pollutant  problems*  Including  photochemical -smog  formation  and  visibility 
degradation*  are  needed.  The  above  algorithms  should  be  incorporated  into  a 
general  computational  package  for  either  research  or  regulatory  applications* 
as  needed. 


5.4.6  Revisions  Needed  to  Enhance  Usability 

One  approach  to  making  AVAP  more  usable*  especially  to  a  broader 
spectrum  of  users*  would  be  to  redesign  the  computer  code  so  that  It  could  be 
installed  on  microcomputer  or  minicomputer  systems.  This  would  require  that 
the  emission  and  dispersion  parts  of  the  code  be  separated  and  that  each  part 
be  further  subdivided  Into  stand-alone  modules.  Bach  module  could  be  designed 
to  read  in  data  from  the  user  terminal  or  from  one  or  two  external  storage 
devices  and  to  write  data  out  onto  one  or  two  output  date-storage  devices  that 
could  be  accessed  by  subsequent  modules.  An  interactive  data-flle  editor 
could  be  designed  to  further  simplify  the  process  of  creating  new  data  files 
or  editing  existing  ones.  Using  the  interactive  mode  for  data-flle  operations 
and  the  batch  mode  for  running  co^utatlonal  codes  could  produce  the  optimum 
combination  of  ease  of  use  and  tlmt  saving. 

A  proposed  new  version  of  the  AVAP  model  designed  along  the  lines 
indicated  above  la  discussed  further  in  Sec.  9  of  this  report. 
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6  SUMMARY  AND  EVALUATION  OF  THE  AQAM  COMPUTER  MODELS 


6.1  THE  AQAM 

The  AQAM  it  Intended  for  use  In  evaluating  the  air-quality  Impacts  of 
large  military  aviation  facilities.  As  In  the  case  of  their  civilian  counter¬ 
part  (AVAP) ,  the  AQAM  computer  codes  are  user-oriented  codes  with  great 
flexibility.  The  user  can  elect  to  treat  virtually  any  conceivable  condona¬ 
tion  of  aircraft ,  air-base,  and  environ  sources  likely  to  be  encountered  in  or 
around  an  air  base.  Aircraft  sources  in  particular  can  be  treated  in 
considerable  detail  —  greater  detail,  in  fact,  than  In  the  AVAP  aodel.  In 
contrast  to  Its  civilian  counterpart,  AQAM  is  composed  of  several  stand-alone 
but  compatible  computer  codes: 

a  Source-Emisalon-lnventory  Model 

a  Short-Term  Emisslon/Dlsperslon  Model 

a  Long-Term  Dispersion  Model  (Research  Version  and  Appli¬ 
cations  Version) 

The  overall  structure  of  the  model  is  shown  In  Fig.  3.  The  dashed 
lines  are  used  to  distinguish  the  four  separate  component  parts,  which  Include 
the  three  separate  computer  codes  listed  above  and  a  "Meteoro log leal  Data 
Program**  that  la  operated  on  request  by  the  Air  Force  Weather  Service  (ETAC). 

The  Short-Term  Emisslon/Dlsperslon  Model  is  used  to  calculate  hourly 
average  source-emission  rates  and  pollutant  concentrations.  It  utilises 
essentially  the  same  point-,  area-,  and  line-source  dispersion  algorithms  as 
the  AVAP  model.  However,  various  updates  and  changes  over  the  years  have 
resulted  In  some  minor  differences  In  some  of  the  algorithms. 

The  Long-Term  Dispersion  Model  has  no  equivalent  In  the  AVAP  model.  It 
employs  a  statlstlcai-cllmatoioglcal-dlspersion  approach,  as  opposed  to  the 
hour-by-hour  approach,  to  compute  long-term  average  pollutant  concentrations 
on  a  monthly  or  annual  basis.  Such  averages  can  be  computed  for  several 
distinct  daily  time  intervals,  as  shown  in  Table  3.  As  Indicated  in  Fig.  3, 
the  Long-Term  Dispersion  Model  requires  msteorologlcsl  input  data  prepared 
specifically  for  this  purpose  by  STAC. 

There  ere  two  versions  of  the  Long-Term  Dispersion  Model:  the  Research 
Version  end  an  abridged  version  referred  to  as  tha  Applications  Varslon.  The 
abridged  version  is  somewhat  less  flexible  but  requires  substantially  less 
computsr  run  tlms  then  the  Research  Version. 

The  Source -Emission-Inventory  Model  is  also  e  physically  separate 
computer  code.  It  operates  on  user-input  source  date  end  produces  e  computer 
file,  containing  source  Information  end  annuel  average  emission  rates,  that 
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Table  3  Definition  of  Diurnal  Time 
Intervals  Used  in  the  AQAM  Long- 
Term  Dispersion  Model 


Period 

Tlae  Interval 
(LSI) 

Interpretation 

I 

0000  -  2400 

all  hours  of  day 

2 

0600  -  1800 

business  hours 

3 

0600  -  0900 

morning  rush  hours 

4 

0900  -  1300 

midday 

5 

1300  -  1800 

evening  rush  hours 

6 

1800  -  2100 

late  evening 

7 

2100  -  0600 

nighttime 

can  be  used  as  an  end  product  or  as  input  to  one  of  the  AQAM  dispersion 
models.  The  Source-Emission-Inventory  Model  organises  sources  into  aircraft, 
air-base  nonaircraft,  and  environ  source  types,  as  does  the  AVAP  model.  It 
also  defines  point,  line,  and  area  source  geometries.  However,  the  structures 
of  these  two  codes  differ  In  accordance  with  operational  differences  at 
military  and  civilian  facilities.  Different  aircraft  types,  operational 
modes,  and  engine  thrust  settings  are  defined  in  these  codes,  and  the  AQAM,  in 
particular,  places  greater  emphasis  on  aircraf t -type-dependent  operational 
parameters.  Whereas  runway,  and  therefore  taxiway,  assignments  are  based  on 
aircraft  range,  wind  direction,  and  airline  in  the  AVAP  model,  they  obey  a 
different  set  of  constraints  at  military  facilities.  The  dependence  of  runway 
choice  on  wind  direction,  for  example,  la  specified  by  the  user  in  the  form  of 
a  series  of  keys,  that  is,  a  sequence  of  ones  and  zeros,  indicating  whether 
the  runway  is  used,  for  each  of  16  wind-direct  ion  sectors.  This  is  a  simpler 
input  schema  than  the  wind-quadrant  procedure  used  in  the  AVAP  model.  Also, 
rather  than  specifying  s  taxlway-segment  aircraf t-act ivity  matrix  as  in  the 
AVAP  model,  AQAM  employs  taxiway  trajectories,  one  for  each  runway-end  and 
aircraf t -parking -area  combination.  Each  such  trajectory  may  contain  several 
straight-line  taxiway  segments.  This  scheme  Is  better  suited  to  military 
facilities  and  is  very  convenient  for  the  input-data  compiler  but  is  not 
computationally  optimal.  In  addition,  military  training  flights  and  espe- 
c  dally  **t  ouch-go**  operations  have  no  equivalent  in  commercial  airports  and 

therefore  require  apeclel  treatment  in  AQAM.  Additional  special  military 
sources  treated  in  AQAM  include  training  firee  end  engine  test-stand 
I  operations. 


'  * 


6.2  EXISTING  VERSIONS  OF  AQAM 
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The  original  AQA N,  as  11 lust rat ad  In  Fig.  3,  was  deacrlbed  In  detail  In 
a  technical  docunent1  dated  February  1975.  An  operator 'a  guide1  vaa 
publlehed  In  1976,  and  detailed  deacrlptlone  of  the  coaputer  codec  thenaelvea, 
Including  flowcharts  and  aubroutlne  deacrlptlone,  were  published  in  three 
separate  docuaants2  2  in  1977.  Between  the  tlae  that  the  original  technical 
docuaentatlon  was  drafted  and  the  coaputer-code  docuaantatlon  waa  drafted, 
soaa  changes  were  aade  In  the  codes.  Because  no  atteapt  waa  asde  to  update 
the  original  technical  docuaentatlon,  soaa  Inconsistencies  exlet  between  the 
coaputer-code  and  technical  docuasntatlons. 

After  1977,  the  Short-Tare  Salas Ion /Dispersion  and  Source-Eel as Ion- 
Inventory  codes  were  used  priaarily  at  Tyndall  Air  Force  Base  (AFB)  for 
environaental  lapact  studies,  at  Mlrsaar  Naval  Air  Base  for  studies  of  naval- 
aircraft  Inpact  and  for  validation  studies,  and  at  Argonna  National  Laboratory 
for  validation  and  sensitivity  studies  and  for  developnent  of  reflneaents  and 
iaproveasnts.  Each  of  these  groups  aade  aodlflcatlona  to  the  original 
coaputer  code  to  aset  their  nlsslons,  so  there  are  now  three  separate  versions 
of  this  aodel.  No  atteapt  will  be  aade  here  to  review  the  Navy  version  of  the 
AQAM  (see,  for  exaaple.  Ref.  23),  but  the  other  versions  will  be  discussed  In 
coaparlaon  with  the  original  version  In  eoas  detail  below. 


6.2.1  Current  Version  of  AQAM  In  Use  at  Tyndall  AFB  (circa  February  1982) 

This  version  of  the  Short -Tara  Ealseion/Dlspersion  Model  la  essentially 
the  saas  as  that  described  in  Ref.  20.  However,  llne-by-llne  coaparisons  have 
revealed  that  soaa  undocuasnted  changes  have  been  aade  In  the  Tyndall  version. 
These  Include  minor  changes  In  the  Inputs  and  outputs,  additions  of  soaa 
explanatory  consents ,  and  other  minor  coding  changes  of  no  consequence.  One 
change  that  will  affect  the  computations  la  the  Incorporation  of  a  terrain- 
correction  factor. 


6.2.2  Valldatloo-Btudy  Version 

Argonna  National  Laboratory  prepared  two  slightly  different  versions  of 
the  AQAM  coaputer  code  for  purposes  of  perforates  calculations  sod  coopering 
results  with  observations  asde  at  Mllllaas  ATI.2  *  These  versions,  AQAM  I 
and  AQAM  II,  are  briefly  described  In  Refs.  24  and  25.  AQAM  1  was  set  up  to 
opsrsts  la  precisely  the  saas  way  that  the  original  AQAM  aoraslly  would  bo 
operated  by  a  user.  Measly.  It  accepts  annual  average  eaisslea  data,  which 
are  subsequently  reduced  to  hourly  average  emission  rates  using  built-in 
algorlthas  that  require  as  input  various  tlaa-of-year,  day-of -weak,  sad 
t las -of -day  ealsslon-dlstrlbutlon  paraastars.  AQAM  11,  on  the  other  head,  woe 
set  up  to  utilise  obeerved  hourly  aircraft  activity  data  taken  at  Nil liana 
AFB.  Mth  regard  to  the  dispersion  algorlthas,  AQAM  1  and  AQAM  II  are 
identical.  However,  AQAM  1  and  II  are  not  identical  to  the  version  described 


reran* 
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In  the  original  code  documentation  mentioned  above*  The  differences  are 
listed  below: 

•  Different  treatment  of  calm  conditions* 

•  Different  output  structure* 

•  Use  of  a  modified  line-source  dispersion  algorithm, 
together  with  a  more  accurate  "error  function"  subroutine* 

•  Other  minor  differences  in  coding* 

The  different  treatment  of  calm  conditions  constitutes  a  major 
Improvement  In  the  Short-Term  Emission/Dispersion  Model*  Unfortunately, 
documentation  is  presently  limited  to  a  technical  description  of  the  calm 
algorithm  in  draft  form  and  a  computer  code  listing*  The  use  of  a  slightly 
modified  line-source  algorithm  and  a  more  accurate  ‘error  function"  resulted 
In  a  substantial  improvement  in  the  accuracy  of  the  code  under  certain  rather 
special  conditions  but  required  only  minor  changes  to  the  original  computer 
code. 


6.2.3  Other  Changes 

The  only  other  major  change  of  note  has  been  the  development  of  a 
completely  new  line-source  algorithm*  This  algorithm  has  not  been 
Incorporated  into  any  version  of  the  Short-Term  Emlsslon/Dlsperslon  Model  as 
of  yet,  and  it  has  not  been  formally  published,  although  a  technical 
description  of  the  algorithm  does  exist  in  draft  form  and  a  code  listing  is 
available*  This  development  results  in  a  substantial  Improvement  in  the 
computer  run  time  over  that  of  the  original  algorithm,  as  well  as  a 
significant  Increase  in  accuracy*  The  Increase  in  accuracy  is  only  really 
significant  for  certain  rather  special  combinations  of  input  parameters  that 
resulted  in  erroneous  output  from  the  original  code* 


6*3  COMPARISON  OF  THE  VARIOUS  VERSIONS  OF  THE  SHORT-TERM  AQAM 

For  purposes  of  this  comparison,  the  three  versions  of  the  Short-Term 
AQAM  shall  be  referred  to  as  follows: 

Version  1:  Original  documented  version  of  the  Short-Term 
Model 


l  * 


Version  2:  Version  used  in  the  Williams  APB  validation  9tudy 
Version  3:  Present  Tyndall  APB  version,  February  16,  1962 
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6*3.1  Treatment  of  Low  Wind  Speeds 

Versions  1  snd  3  both  sssu mm  that  If  the  triad  speed  is  less  then  1  mph, 
the  wind  speed  end  direction  ere  set  to  the  nearest  nonsero  hourly  value* 
Version  2  uses  the  calm  algorithm* 

6*3*2  Treatment  of  Lid  Height  or  Mixing  Depth 

Versions  1  and  3  both  require  hourly  values  as  input*  Version  2  uses 
either  the  Nosakl  equation  (see  Ref*  25,  p.  60)  or  hourly  input  values  froa 
acoustic  sounder  data  obtained  at  Williams  AFB.  If  the  trlnd  is  calav  Version 
2  assuaes  pollutants  are  unlforaly  mixed  in  the  vertical  direction  within  the 
alxed  layer* 


6.3.3  Error  Function 

All  three  versions  require  the  use  of  an  error-function  subroutine  that 
is  called  by  the  line-source  algorithm.  Versions  1  and  3  use  the  same  error- 
function  subroutine*  Version  3  uses  a  superior  subroutine*  In  addition, 
there  is  also  a  minor  change  in  the  line-source  algorithm  Itself  in  Version  2 
that  corrects  an  error  that  occurs  only  for  certain  combinations  of  input 
variables. 


6.3*4  Terrain  Correction  Factor 

Version  2  employs  a  terrain  correction  factor;  the  other  two  versions 

do  not* 

6*3*5  Plums  Rise 


Briggs fs  plume-rise  formula  has  been  added  as  an  option  to  Version  3 
but  not  to  the  other  two* 


6.3*6  Treatment  of  Large  Area  Sources 

All  three  versions  treat  large  area  sources  using  the  same  algorithm* 
However,  this  algorithm  is  not  documented  in  the  technical  report  dated 
February  1975*  It  is  an  improved  treatment  compared  with  that  described  in 
the  February  1975  document* 
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6.3.7  Input /Output  Structure 

There  are  iom  differences  in  the  Input/output  statements  appearing  In 
the  three  codes.  These  differences  do  not  affect  the  operation  of  the  codes 
themselves. 


6.4  AQAM  EVALUATION 


6.4.1  Intended  Uses  and  Strengths 

The  AQAM  package  was  originally  designed  for  use  on  large*  fast*  main¬ 
frame  computer  facilities  by  environmental  and  computer  staff  at  various  Air 
Force  Installations .  It  was  to  be  used  to  compute  emissions  and  air  quality 
on  either  a  short-term  or  long-term  basis.  The  models  are  very  flexible  and 
can  be  used  for  a  wide  variety  of  air-quality  problems.  Separate  emission  and 
dispersion  computations  can  be  performed.  Considerable  freedom  of  choice  Is 
available  to  the  user  regarding  types  of  sources  to  be  Included  in  a 
particular  calculation  and  the  level  of  detail  of  treatment.  At  the  time  of 
development*  the  computer  codes  incorporated  state-of-the-art  emissions  and 
dispersion  modeling  techniques  for  nonreactlve  pollutants  (over  time  Intervals 
of  one  hour)  and  for  transport  distances  of  the  order  of  several  kilometers 
away  from  the  military  installation.  Provisions  for  treating  a  wide  variety 
of  aircraft,  nonaircraft  air-base*  and  environ  sources  were  incorporated.  The 
codes  were  designed  to  accept  meteorological  data  routinely  collected  at  air 
bases  and  processed  by  ETAC.  The  Long-Term  Model*  In  particular*  operates  on 
a  climatological  Joint  stability  and  wind-rose  data  set  prepared  specifically 
by  ETAC  for  that  purpose.  The  original  version  of  the  AQAM  computer  codes  is 
well  documented*  and  the  Short-Term  Model*  especially*  has  undergone  extensive 
testing.  Tests  have  Incliuled  application  to  both  the  Washington  National 
Airport2  and  Williams  AFB2  ’  and  extensive  sensitivity  analysis.  The 

Long-Term  Model  has  been  tested  using  hypothetical  data  bases  only.  Both  the 
research  and  the  applications  versions  of  the  Long-Term  Model  have  been  tested 
using  hypothetical  data  bases  and  have  been  compared  against  each  other  and 
have  undergone  some  sensitivity  analyses • 


6.6.2  Limitations  and  Weaknesses 


Model  Usability 

The  generality  and  flexibility  built  Into  the  original  AQAM  computer 
codes*  along  with  the  ability  to  treat  aircraft  operations  at  Air  Force  bases 
in  great  detail*  has  generated  two  problems  that  have  limited  the  codes'  use 
at  various  Air  Force  Installations.  First*  the  potential  user  tends  to  be  put 
off  by  being  confronted  with  having  to  proceed  step-by-step  through  the  user's 
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guide,  which  requires  hla  to  asks  a  large  nuaber  of  technical  declalona  and  to 
gather.  If  the  problea  warrants  It,  a  large  body  of  inforaation;  this  can  be  a 
very  aanpower-lntenslve  process.  Of  course,  if  the  scope  of  the  problea  is 
saall,  the  input  inforaation  requireaants  are  correspondingly  saall. 
Nevertheless,  even  if  the  problea  is  a  staple  one,  the  novice  user  is  still 
faced  with  what  ha  aay  regard  as  an  awesoaa  task  of  becoalng  f  sal  liar  with  the 
general  input  structure  and  operation  of  the  coaputcr  codes.  The  other 
problea  is  related  to  coaputer  facility  requlreaants.  If  the  nuaber  of 
sources,  especially  line  sources,  used  in  the  problea  is  large,  then  the 
coaputer  run  tlae,  especially  on  older  facilities,  can  be  excessively  long. 
Furtheraore,  regardless  of  the  scope  of  the  problea,  the  coaputer  core-storage 
requlreaent  Is  large  and,  in  fact,  often  too  large  for  the  aany  coaputer 
installations  at  which  severe  constraints  are  placed  on  users  for  one  reason 
or  another. 


Docuasntatlon 

Although  the  original  AQAM  has  been  fully  docuaented,  a  nuaber  of  code 
changes  and  updates  have  been  iapleaanted  by  Air  Force  personnel  In  the 
Tyndall  version  and  by  Argon ne  staff  in  the  version  used  in  the  Vllllaas  AFB 
Validation  Study.  These  changes,  especially  those  lapleaented  at  Argonne, 
affect  model  performmnce  under  certain  conditions.  The  published  docuaen- 
tatlon  has  not  been  correspondingly  updated. 


Technical  Issues 


The  technical  Ualtations  and  weaknessas  of  the  AQAM  dispersion  codes 
are  essentially  the  saae  as  those  of  the  AVAP  coaputer  codes  (see  Sec.  5.4.2) 
with  the  following  exceptions: 


e  The  AQAM  codes  have  not  been  updated  to  include  the  results 
of  the  Dulles  and  Washington  National  airport  experiaants 
regarding  Jet-aircraft  pluas  rise  and  initial  dispersion. 

e  The  AQAM  Short-Tern  Model  (the  version  used  in  the  Will  leas 
AFB  Validation  Study)  contains  a  cal a  algorltha,  whereas 
the  AVAP  aodel  does  not.  AQAM  also  contains  a  superior 
"error  function"  subroutine  that  allainates  aoaa  of  the 
nuaerlcal  probleaa  arising  froa  certain  coabl net ions  of 
Inputs  to  the  AVAP  code. 


e  The  Wi lliaaa  APB  Validation  Study  version  also  contains  the 
default  option  of  using  the  Nosakl  equation  to  coaputa 
nixing  depth.  However,  this  equation  has  been  shown  to  ba 
inadequate.  A  better  default  option  la  required.  The 
default  algorltha  used  in  the  AVAP  aodel,  although  not 
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adequate  either,  la  probably  superior  to  the  Noaakl 
equation* 

e  The  level  of  treatment  of  aircraft  and  air-base  sources 
available  to  the  user  In  the  AQAM  is  probably  unnecessarily 
detailed*  For  moat  purposes,  the  number  of  source  types 
could  be  reduced  and  the  use  of  aircraft-dependent 
operational  parameters  could  be  replaced  with  the  use  of 
generic  operational  parameters* 


6*4.3  Updates  Essential  for  Inert-Pollutant  and/or  Screening  Applications 
See  Sec*  5.4.3* 

6.4.4  Additional  Updates  Desirable  for  Inert-Pollutants  and/or  Screening 
Applications 

See  Sec*  5*4.4.  In  addition,  the  documentation  should  be  updated,  and 
minor  differences  among  existing  versions  of  the  AQAM  should  be  eliminated* 


6.4.5  Updates  Weeded  for  Reactive-Pollutant  Modeling 
See  Sec •  5 .4 . 5 . 


6*4*6  Revisions  Needed  to  Enhance  Usability 

As  with  the  AVAP  computer  code,  the  best  way  to  improve  the  usability 
of  the  AQAM  codes,  especially  to  make  it  useful  to  a  broader  spectrum  of 
users,  Is  to  redesign  the  codes  so  that  they  can  be  operated  on  microcomputer 
or  minicomputer  systems*  Although  the  emission  and  dispersion  computer  codes 
for  the  Short-Term  Model  are  already  separated,  the  separation  is  Incomplete* 
tor  example,  emission  rates  dependent  on  the  specific  hour  of  simulation  or  on 
the  associated  meteorology  for  that  hour  are  still  computed  within  the  short¬ 
term  emisslon/dlsperslon  computer  program.  Furthermore,  and  even  more 
Important,  the  structure  of  the  AQAM  codes  and  the  else  of  the  arrays  used  to 
store  Input  data  preclude  the  use  of  computers  having  only  modest  core 
storage*  Hence,  a  course  of  action  similar  to  that  briefly  outlined  in  Sec* 
5*4*6  for  the  AVAP  model  could  be  taken  for  the  AQAM  as  well*  Proposed  new 
designs  for  both  the  AVAP  model  and  AQAM  are  discussed  In  Sec*  9* 


7  MODEL  COMPARISONS  AND  RECOMMENDAT I ONS 


Although  the  AVAP  Model  end  AQAM  started  off  having  virtually  the  same 
dispersion  algorithms*  these  algorithms  have  evolved  along  somewhat  separate 
paths  due  to  various  updates  and  refinements.  The  purpose  of  the  present 
comparison  is  to  point  out  the  similarities  and  differences  and  to  select* 
where  possible*  the  best  that  each  model  has  to  offer.  These  "best"  offerings 
could  then  be  used  to  guide  the  development  of  a  new  composite  or  joint 
mllltary/civlllan  dispersion  package.  For  some  features  of  the  dispersion 
problem*  neither  model  treatment  may  be  completely  satisfactory.  In  such 
cases,  other  alternatives  are  suggested. 


For  the  purposes  of  this  comparison*  the  version  of  the  AVAP  model  used 
In  the  Updated  Assessment  of  Air  Quality  Impacts  at  Major  U.S.  Airports^** ** 
and  the  version  of  the  AQAM  used  in  the  Williams  AFB  Validation  Study  (AQAM  I 
Is  technically  equivalent  to  AQAM  II)  will  be  used. 


Table  4  summarises  the  methods  used  by  each  of  these  two  models  to 
treat  each  feature  of  the  modeling  problem.  Recommendations  for  the  "best** 
method  of  treating  each  feature  are  also  given  In  the  table. 
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8  WHICH  DIRECTION  SHOULD  FUTURE  EFFORTS  TAKE?  —  A  DECISION  TREE 


A  combination  of  changing  operational  and  technical  requirements  for 
models  and  computer  codes;  advances  In  the  state  of  the  art  of  modeling; 
availability  of  modern,  low-cost  microcomputers  and  minicomputers;  and  numer¬ 
ous  Inconsistencies  between  various  currently  in-use  versions  of  AVAP  and  AQAM 
computer  codes  and  documentation  raises  a  number  of  Interrelated  questions 
regarding  where  future  efforts  should  be  headed* 

For  example,  should  research  and  development  (RAD)  tasks  be  undertaken 
to  advance  the  state  of  the  art  of  emission  and/or  dispersion  modeling  or 
should  future  efforts  be  directed  more  toward  applications  issues  and,  in 
particular,  toward  Improvement  of  model  usability? 

Should  the  Inconsistencies  that  have  been  pointed  out  between  various 
computer  codes  and  their  documentation  be  resolved  now  or  should  such  updating 
be  postponed  until  already-developed,  superior  algorithms  are  Incorporated 
Into  the  codes?  Alternatively,  should  the  AQAM  and  AVAP  packages  continue  to 
be  supported  as  separate  systems  or  should  a  joint  package  that  takes  advan¬ 
tage  of  the  best  elements  of  both  systems  be  developed?  Should  that  joint 
package  Incorporate  all  of  the  latest  algorithms?  Should  further  efforts  to 
Improve  or  update  the  current  large  models  be  suspended  in  favor  of  designing 
new  versions  specifically  for  modem  microcomputers  and/or  minicomputers? 

These  Issues  and  others  point  to  the  need  for  a  systematic  decision¬ 
making  procedure*  One  approach  to  such  a  systematic  procedure  la  to  formulate 
a  "decision  tree"  that  reveals  the  various  options  and  their  consequences* 
Figure  4  shows  a  first  attempt  at  a  decision  tree  developed  for  the  purpose  of 
sorting  out  sow  of  the  options  for  proceeding  with  the  AVAP  model  and  AQAM* 
No  attempt  Is  made  to  display  all  of  the  possible  options  or  to  even  fully 
characterise  the  options  that  are  given*  Rather,  this  decision  tree  should  be 
regarded  simply  as  a  guide  to  future  decision  making.  Alternative  branches 
can  be  readily  added  and  further  expanded. 

The  decision  tree  contains  question  marks  that  are  located  at  branch 
points  or  decision  points.  The  straight  lines  radiating  from  each  question 
mark  indicate  the  alternative  paths  that  various  decisions  will  lead  to.  For 
example,  on  the  first  diagram,  the  first  question  mark  requiring  a  decision 
refers  to  the  question  of  whether  to  support  applications  or  RAD-related 
work.  The  line  going  to  the  left  shows  the  "applications  branch."  Following 
along  the  "applications  branch,"  the  next  decision  concerns  the  choice  between 
"user"  and  "apeclallst"  models.  Following  the  "user  branch,"  the  next 
decision  concerns  the  choice  between  separate  AVAP  and  AQAM  packages  and  a 
Joint  Air  Quality  Modeling  Package  (JAQMF),  etc.  The  circled  letters,  which 
appear  at  the  lower  extremities  of  each  of  the  branches,  refer  to  continua¬ 
tions  that  are  shown  on  subsequent  pages. 
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9  OUTLINE  Of  A  PROPOSED  NEW  DESIGN  POR  THE  AVAP  AND  AQAM  EMISSION  MODELS 


A  new  coaputatlonal  ijitw  consisting  of  •  sat  of  stand-alone  coaputer 
progress  la  being  recoasended  for  use  In  caees  in  which  the  scope  of  the 
problea  does  not  warrant,  or  available  coaputer  facilities  do  not  perait,  the 
use  of  the  full  AVAP  aodel  or  AQAM.  To  date,  only  the  ealselon  portion,  or 
"front  end,"  of  the  eystea  has  been  considered  in  any  detail. 

A  descriptive  outline  of  the  ayetea  is  presented  in  this  eectlon, 
together  with  aoae  aacro  flowcharts  Illustrating  its  overall  structure  and 
aajor  coaponents.  This  eectlon  also  contains  aoae  discussion  of  the  rationale 
leading  to  the  design  choices  eabodled  in  the  new  eystea.  Detailed  flowcharts 
and  descriptions  of  data-file  structures  and  contents  are  contained  in 
Ref.  30. 


The  overall  structure  of  the  civilian  and  ailltary  versions  of  the  new 
eystea  is  sufficiently  siallar  that  only  a  single,  generic  systea  needs  to  be 
described.  When  necessary,  distinguishing  features  applicable  only  to  the 
ailltary  or  civilian  versions  will  be  noted. 

The  new  coaputatlonal  eystea  is  designed  in  a  aodular  fashion  to 
accoaaodata  saall  coaputara,  while  at  the  aaas  tlae  providing  the  user  with 
the  flexibility  to  treat  a  variety  of  source  types  and  configurations  and  to 
per fore  either  ealsslons  or  pollutant-concentration  calculations.  Each  aodule 
or  subprograa  is  designed  to  read,  at  aost,  one  or  two  input  data  files  and 
write,  at  aost,  one  or  two  output  data  files.  New  files  can  be  created  froa 
scratch  or  old  data  filaa  can  be  edited  as  required.  For  a  given  type  of 
problea,  only  those  aubprograas  and  corresponding  data  files  that  are  required 
ere  actually  used.  This  eaves  a  considerable  aaount  of  core  storage  as 
coapared  with  the  full  AVAP  or  AQAM  coaputer  codes. 

The  esdselon  and  dispersion  portions  of  the  systea  are  coapletely 
separated  froa  each  other,  ellalnatlng  the  need  to  provide  core  storage  for 
both  when  only  one  portion  aay  actually  be  required.  In  the  AVAP  aodel,  both 
portions  are  coablned  into  a  single  code.  In  AQAM,  the  source  ealselon 
Inventory  is  coaplled  and  annual  average  ealeslon  rates  are  coaputed  in  the 
Source-Baisslon-Inventory  Model  coaputer  code.  The  reaalnlng  short-tern 
ealselon  and  dispersion  aodel  calculations  are  perforaed  in  the  8hort-Tera 
Ealselon/  Dispersion  Model  coaputer  code.  Figure  3  illustrates  the  AQAM 
Short-Tera  Model  coaputer-code  structure.  All  of  the  subroutines  contained 
within  the  dotted  line  are  associated  with  the  coaputation  of  short-tera 
aulas Ion  rates.  These  subroutines  or  their  equivalents  would  be  placed  in 
separate  aodules  in  the  new  coaputatlonal  systea.  Purtheraore,  whereas  it  la 
now  necessary  to  place  the  entire  aalsslon  inventory  into  core  to  coapute  the 
pollutant-concentration  contribution  froa  each  source,  in  the  now  systea  only 
one  source  will  occupy  the  core  at  a  tlaaf  This  will  result  in  a  very 
substantial  reduction  in  core-storage  requlreasnts. 


For  user  convenience,  s  special  nodule  is  proposed  that  will  create  or 
edit  the  baste  Input  data  files.  It  will  be  operable  in  the  interactive  mode, 
while  the  sain  computational  routines  will  be  operable  in  the  batch  node. 
This  nixed -node  systen  will  simplify  the  user's  task  of  compiling  the  input 
data  files,  while  avoiding  the  necessity  of  long  sessions  at  the  computer 
ternlnal  during  the  longer  run  tines  required  by  the  computational  programs. 


9.1  OVERALL  STRUCTURE  OF  THE  NEW  SYSTEM 

An  overview  of  the  new  computational  system  is  shown  in  Fig.  6.  Note 
that  the  user  is  expected  to  have  access  to  a  microcomputer  or  minicomputer 
system  that  utilizes  at  least  two  permsnent  storage  devices  (e.g.,  disk 
readers/writers).  This  hardware  requirement  is  necessary  for  the  efficient 
interactive  processing  of  input  and  output  data  files.  Once  these  files  have 
been  prepared,  they  are  subsequently  used  by  the  computational  computer 
routines  to  generate  either  emission  rates  or  ambient  concentrations.  It  is 
desirable,  but  not  essential,  to  have  access  to  a  mainframe  computer  for  doing 
the  batch  mode  emission  and  dispersion  calculations.  (Microcomputers  or 
minicomputers  may  be  used  if  the  number  of  sources  to  be  analyzed  is  small.) 


Fig.  6  Interactive  and  Batch-Mode  Computer  Operations  — 
An  Overview  of  the  Computational  System 


Figure*  7a  and  7b  show  macro  flowchart*  of  the  ahort-term  (hourly) 
emission*  and  dispersion  portions  of  the  system,  respectively*  In  both  charts 
the  malor  processing  steps  are  shown  to  the  left  In  rectangular  boxes ,  with 
the  processing  sequence  proceeding  downward  as  Indicated  by  the  arrows.  The 
boxes  and  arrows  shown  on  a  particular  horizontal  line  represent  the  Input  of 
data  from  data  files  and  the  processing  of  data  to  create  or  edit  the  data 
files  needed  to  carry  out  the  malor  steps  shown  to  the  left.  Note  that  the 
various  data  files  are  identified  by  numbers  for  convenience. 

Each  of  the  data-processlng  steps  along  a  horizontal  line  or  the  malor 
steps  In  the  vertical  line  can  be  performed  as  separate.  Isolated  computer 
operations  or  clustered  together  as  dictated  by  the  optimal  use  of  the  user’s 
computational  facilities.  For  example,  beginning  on  the  first  line  of  the 
first  chart,  the  user  may  elect  to  input  data  via  the  Interactive  Data-File 
Program  (I.O.F.P.)  and  create  File  l  In  a  single  operation.  He  may  then 
create  File  2  in  another  atep.  Then  he  can  exercise  the  aircraft-source- 
emissions  code  that  requires  File  1  and  File  2  as  Inputs  and  that  generates 
File  7  as  an  output  In  a  single  batch-mode  operation. 

Note  that  at  the  end  of  the  second  chart,  the  user  is  given  several 
output  options.  In  general,  the  user  should  be  given  the  choice  of  outputting 
individual  or  accumulated  source  contribut ions  to  pollutant  concentrations  at 
each  receptor. 


9.2  LONG-TERM  EMISSION  ESTIMATES 

The  macro  flowcharts  shown  In  Figs.  7a  and  7b  Illustrate  only  the 
short-term  (hourly)  emission  and  dispersion  calculations.  In  addition,  it  may 
be  desirable,  for  emission-reporting  purposes,  to  also  compute  annual  average 
emissions.*  In  some  cases,  only  annual  average  emissions  may  be  required. 
For  the  sake  of  computet  tonal  efficiency.  It  is  worthwhile  to  have  separate 
computer  programs  for  estimating  hourly  and  annual  average  emissions.  Note 
that  simply  scaling  hourly  emissions  up  or  annual  emissions  down  may  not  be 
satisfactory,  since  the  hourly  average  emission  rates  of  Interest  may  not  be 
typical  for  the  entire  year. 

Figure  8  Illustrates  a  straightforward  procedure  for  structuring  an 
annual-average-emission  computation  routine.  Note  that  such  a  routine 
requires  Che  use  of  LTO-cycle-t ype  emission  factors. 


9.3  DATA  FILES  AND  THE  INTERACTIVE  DATA-FILE  PROGRAM 

The  new  computational  system  makes  extensive  use  of  dete  files.  These 
files  are  the  meent  by  which  one  module,  representing  one  step  in  the 


*Thls  has  bean  a  traditional  requiremant  for  military  purpoaaa 
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Fix*  7a  Macro  Flowchart  of  Short -Tara  Salsa Ion  Coaputatlona 
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READ  fLTO'S  FEE  A/C  TYPtf 


READ  LTO  EMISS.  FACTORS 


'fLTO'S  FIR  A/C  TTFE  FER  TEAR 


FILE  1  PERMANENT  A/C  DATA 


COMPUTE  A/C  4  SERVICE- 
VEHICLE  EMISSIONS 


READ  ACCESS- VEHICLE 
EMISSION  FACTORS 


FILE  4  A.V •  EMISS.  FACTOR  DATA 


READ  ACCESS-VEHICLE 
INPUT  DATA 


COMPUTE  ACCESS- 
VEHICLE  EMISSIONS 


OUTPUT 


'#  PASSENGERS  ARRIVING  AT 

AIRPORT  VIA  ACCESS  VEHICLES 
FER  LTO  OF  TYPE  I  A/C 

#  PASSENGERS  DEPARTING 

AIRPORT  VIA  ACCESS  VEHICLES 
PER  LTO  OF  TYPE  I  A/C 

X  PASSENGERS  USING 
CARS  (TAXIS) 

BUSES 

OTHER 


A/C  -  aircraft 


#  VEHICLE  MILES  (VM'S)  PER 
VEHICLE  TRIP  THROUGH 
AIRPORT 

#  MIN.  OF  IDLE  TIME  PER 
VEHICLE  TRIP  THROUGH 
AIRPORT  PER  VEHICLE  TYPE 

X  CARS  PARKING  FOR  >1  HR 

X  BUSES  PARKING  FOR  >1  HR 

AVG.  VEHICLE  SPEED 
(EXCLUDING  IDLING) 


Pig*  8  Annual-Avaraga-Eslaalon-Eatlmatlon  Coda 
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computational  eyatea,  paeeea  on  information  Co  «  eubeequent  stop.  It  la,  in 
face,  the  uaa  of  data  fllaa  Chat  anablaa  tha  othervlae  lengthy  eoaputar 
prograa  Co  ba  brokan  down  into  aaall  modules  that  can  bo  handlad  Independently 
of  each  other,  chua  greatly  reducing  core  atoraga  and  computer-run  cl  at 
requlreaenta. 

Two  general  cypaa  of  data  fllaa  are  envleloned  for  uaa  by  the 
coapuCatlonal  ayatea,  peraananc  data  fllaa  and  problea-epedflc  data  fllaa* 
The  flret  type  will  contain  data,  auch  aa  aalaaion  factora  and  other  paraaatar 
valuea,  that  are  not  expected  to  be  changed  vary  often.  The  eecond  type  will 
contain  data,  auch  aa  aource  deacrlptlone  and  activity  at  a  particular 
facility,  that  would  ba  expected  to  change  at  leant  partially  f row  computer 
run  to  computer  run. 

Strictly  aa  a  convenience  to  uaera.  It  la  worthwhile  conelderlng  a 
data-f ile-proceeslng  computer  prograa  that  could  be  uaed  either  to  create  new 
data  fllaa  from  acratch  or  to  modify  or  edit  exletlng  data  fllaa.  For  even 
greater  convenience,  the  computer  prograa  could  be  dealgned  to  operate  In  the 
interactive  node. 

An  Interactive  date-file  prograa  could  operate  aa  followa.  It  would 
coomunlcate  with  a  computer-terminal  uaer  by  ieauing  aeeeagee,  aeklng 
queatlona,  and  making  requeete  for  Input  data.  The  uaer  would  reepond  by 
keying  In  anawera,  coaaanda,  or  data.  The  prograa  ahould  be  able  to  read  in 
data  keyed  In  from  the  ueer'a  terminal  or  from  data  fllaa  on  atoraga 
devlcea.  It  ahould  aleo  be  able  to  Hat  the  contenta  of  data  fllaa  at  the 
computer  terminal,  generate  hard  copy,  or  write  data  fllaa  onto  atoraga 
devlcea.  When  modifying  an  exlatlng  data  file,  the  prograa  ahould  provide  the 
option  of  either  writing  over  the  old  data  file  or  leaving  it  intact  and 
writing  the  modified  data  file  onto  an  entirely  new  apace. 

For  convenience.  It  may  be  worthwhile  for  the  data-file  prograa  to  be 
capable  of  working  with  a  number  of  different  typea  of  data  fllaa.  For 
example.  It  may  be  ueeful  to  have  one  prograa  that  could  creata  or  modify 
aeveral  of  the  data  fllee  Indicated  in  Flga.  7a  and  7b.  To  uaa  a  prograa  with 
thla  capability,  the  uaer  would  elaply  aupply  the  file  type  or  number  upon 
requeet  by  the  data  file  program. 

A  macro  flowchart  of  an  example  of  an  Interactive  data  file  prograa  la 
llluetrated  In  Fig.  9.  A  detailed-example  flowchart  la  given  in  Ref*  30. 


9.4  SOURCE  TYFES  TO  BE  OOMSIOCRED 

Except  in  epeclal  caaee.  It  la  generally  adequate  to  treat  only  the 
major  emit tore  at  aa  aircraft  facility.  Theme  include}  aircraft  operatione 
(eepeclally  ground  operatione),  aircraft  aervlce  vehlclee,  and  civilian  and 
military  acceea  vehlclee. 
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OR  END  SESSION 


WRITE  EDITED  DATA 
FILE  ONTO  NEW  FILE 
OR  OVER  OLD  FILE  SPACE 


EDITED  DATA  FILE 


Fig.  9  Macro  Flowchart  of  Interactive  Data -File  Prograa  (IDFP) 


Ealsslons  f  row  airborne  aircraft  operations  are  Initially  much  wore 
dispersed  than  ground-based  ealsslons  and  therefore  contribute  negllgtbly  to 
ground-level  pollutant  concentrations  except  In  rather  special  cases  (see  Sec. 
10.3).  For  this  reason.  It  Is  not  necessary  to  treat  ealsslons  froa  airborne 
portions  of  aircraft  LTO  cycles  In  as  great  detail  as  ground-based 
ealsslons.  In  addition,  ealsslons  froa  touch-go  operations  at  allttary 
training  bases  will  asks  anch  saaller  contributions  to  ground-level  pollutant 
concentrations  than  standard  LTO  operations.  However,  for  certain  purposes, 
such  as  eat lasting  lapacts  of  aircraft  ealsslons  on  phot ocheal  ca 1-saog 
foraatlon  and  visual-range  reduction,  etc*,  It  aay  be  desirable  to  Incorporate 
ealsslons  estlaates  for  these  operations  and  also  to  provide  seas  crude 
estlaate  of  source  location. 
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In  general,  the  greatest  attention  to  detail  ahouid  be  given  to 
locating  aircraft  ground  operations*  This  includes  aircraft  parking  and 
servicing  areas,  taxiway  segments,  aprons  adjacent  to  the  ends  of  runways, 
where  queuing  occurs,  and  runways* 

In  general,  ground  operations  that  deviate  from  a  staple,  nonstop 
aircraft  aoveaent  should  be  accounted  for,  since  such  operations  determine 
Initial  pollutant  concentration  patterns*  Particular  attention  should  be 
given  to  queuing  at  runways  or  service  areas  and  congestion  that  leads  to 
slower-than-normal  taxiing*  Provisions  for  both  queuing  and  reduced  taxiing 
speeds  should  be  Incorporated  Into  the  model  design* 

Most  other  sources,  either  because  of  their  low  emission  rates  or 
Intermittent  operations,  do  not  make  significant  contributions  to  overall  air 
quality*  However,  in  certain  special  cases  it  asy  be  necessary  to  treat  one 
or  more  of  the  following: 

s  Hydrocarbon  evaporative  losses  from  fuel  storage  and 
handling*  Parked  vehicles  are  a  major  source  of 
evaporative  losses* 

•  Emissions  from  power-generating  or  space -heating/ coo ling 
facilities* 

e  Incinerators* 

s  Engine  test  stands* 

e  Training  fires* 

In  addition,  if  the  problem  involves  toxic  substances,  such  sources 
should  be  given  special  treatment*  Only  emission  factors  for  the  "criteria** 
pollutants  should  be  stored  in  the  computer  codes  and  data  files*  Special 
provision  should  be  made  for  inputting  emissions  from  nonroutine  sources 
directly  rather  than  storing  a  great  variety  of  special  emission  factors* 

It  should  be  borne  in  mind  that,  since  there  is  no  health-related 
standard  for  hydrocarbons  (HCs),  it  is  not  necessary  to  produce  detailed 
dispersion-model  estimates  for  this  pollutant  category*  Of  course,  HC 
emissions  from  combustion  sources  can  be  routinely  produced,  If  needed, 
without  any  significant  additional  effort,  provided  the  necessary  emission 
factors  are  stored  in  the  data  files*  Evaporative  HC  emissions  are  another 
matter*  Fortunately,  studies  at  commercial  airports  indicate  that  HC 
emissions  due  to  evaporation  tend  to  make  up  a  fairly  smell  fraction  of  total 
HC  emissions  due  to  both  combustion  end  evaporation*  Evaporative  loesee 
associated  with  storage  tanks  io  particular  constitute  a  smell  fraction  of 
total  evaporative  losses  at  both  commercial  and  military  facilities* 
Evaporation  from  parked  vehicles  exceeds  that  associated  with  storage  tanks  at 
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both  types  of  facilities.  Also,  depending  on  how  spillage  during  aircraft 
filling  and  fuel  line  venting  Is  handled,  these  nay  constitute  Important 
sources  of  evaporative  HCs.  See  Sec.  10.5  for  further  details. 


► 


9.5  AIRCRAFT  SOURCES  AND  OPERATIONAL  MODES 

The  aircraft  operational  Modes  that  have  been  used  for  Modeling 
Military  and  civilian  aircraft  facilities  are  listed  in  Table  5.  Note  that 
there  are  three  special  Modes  unique  to  air  bases,  namely  idle  during  anting, 
Idle  during  disarming,  and  touch-go  training  flights.  These  latter  operations 
do  not  occur  at  all  military  facilities  and,  when  they  do  occur,  they  are 
generally  limited  to  certain  squadrons  only. 

It  is  important  to  observe  In  Table  5  that,  basically,  only  six  engine 
modes  are  defined  for  both  the  military  and  civilian  aircraft  operations. 
Hence,  the  computer  programs  used  to  compute  aircraft -source  emissions  can  be 
structured  in  essentially  the  same  way  for  both  types  of  facilities. 

For  simplicity,  it  is  recommended  that  emission  factors  be  stored  by 
engine  node  on  a  per  aircraft-type  basis  rather  than  on  a  per  engine-type 
basis. 


It  is  further  recommended  that,  for  purposes  of  alr-quallty  computa¬ 
tion,  only  the  ground-based  aircraft  operations  be  Included.  (See  Sec.  10.3 
for  a  possible  exception  to  this  rule.)  The  airborne  operations  are  known  to 
make  only  very  small  contributions  to  ground-level  pollutant  concentrations. 
However,  because  the  user  may  wish  to  test  this  point,  or  because  he  may  wish 
to  compute  the  emissions  due  to  airborne  operations.  It  Is  advisable  to 
provide  an  option  to  Include  the  airborne  sources  and  computation  of  emissions 
from  airborne  operations.  If  such  an  option  is  provided,  the  airborne  sources 
should  be  restricted  to  one  approach  path  and  one  departure  path  per  runway. 
Each  such  path  could  have  up  to  two  components,  one  extending  from  the  ground 
to  some  nominal  height  of  say  100  m  or  200  m  and  the  other  from  there  to 
perhaps  1000  m.  This  Is  the  way  the  present  AVAP  model  is  designed.  The  AQAM, 
however,  presently  uses  aircraft-type-dependent  approach  and  departure 
paths.  This  latter  degree  of  detail  disproportionately  Increases  the 
computational  requirements  without  attendant  significant  changes  in  accuracy 
of  emissions  or  pollutant -concent rat Ion  calculations. 


With  respect  to  ground-based  aircraft  sources,  It  is  recommended  that 
two  additional  sources  be  considered  for  Inclusion:  first,  one  or  two  line 
sources  per  runway  for  queued  aircraft  and  second,  a  short  taxi way  segment  to 
connect  the  end  of  the  taxiway  with  the  start  end  of  the  runway  takeoff 
roll.  Inclusion  of  these  two  source  types  will  greatly  improve  the  accuracy 
of  the  spatial  distribution  of  the  taxiway  emissions  and  consequently  will 
improve  the  spatial  resolution  of  the  pollutant-concentration  calculations* 
During  hlgh-aircraft-actlvity  periods,  up  to  SOX  or  more  of  the  ground-based 
aircraft  emissions  can  be  attributed  to  such  sources  (primarily  the  line 
sources  representing  the  aircraft  queues)* 


*r  >«Jl 
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Table  5a  Military  Aircraft  Mode  Dafinltlona 


Aircraft  Operational  Mods 

Eng tea  Mods  or 

Thrust  Setting 

Engine  Mods 
Nuwber 

1 

Idle  at  startup 

Idle 

3 

2 

Taxi  before  takeoff 

Idle 

3 

3 

logins  check  at  runway  end 

Military 

4 

A 

tuaway  roll 

Military  or  afterburner 

5 

5 

Cllaboot  Step  1 

Military  or  afterburner 

5 

6 

Cllaboot  Step  2 

Military  or  afterburner 

5 

7 

Approach  Step  1 

lateraadlate 

1 

S 

Approach  Step  2 

Into rand late 

1 

9 

Landing  on  runway 

Mixed 

2 

10 

Taxi  after  landing 

Idle 

3 

11 

Idle  at  shutdown 

Idle 

3 

Special  Modes 

12 

tdle  during  a ruing 

Idle 

3 

13 

Idle  during  dlsarula| 

Idle 

3 

U 

Touch-go  runway  roll" 

Idle  (?) 

3 

13 

Auxiliary  power  unit  operation  — 

6 

*Tooch-go  approach  path  ■  coolant tonal  approach 

Touch-go  cllaboot  path  -  conventional  cllaboot  but  spatially  displaced* 


Table  5b  Civilian  Aircraft  Mode  Def Initione** ^ 


Engine  Mods  or  tnflna  Mods 

Aircraft  Operational  Nods  Threat  Sat  tin#  amber 


1 

Approach  step  1 

Approach 

2 

Approach  atop  2 

Apircoach 

3 

Landing 

Mixed 

4 

Inbound  taxi 

Idle 

3 

Idle  at  shutdown 

Idle 

6 

Idle  at  startup 

Idle 

7 

Outbound  taxi 

Idle 

S 

Manway  roll 

Takeoff 

9 

Cl 1 about  1 

Cllaboot 

10 

Cl 1 about  2 

Cllwbout 

11  Auxiliary  poeer  unit  operatlaa 


<hdsalM  fact  ora  are  stored  for  each  diet  tact  aircraft  type, 
engine  woda  amber  (for  niUtaryi  S  aircraft  ♦  1  ami  11  ary 
power  eatt)t  for  civilian!  5  aircraft  ♦  1  auxiliary  power 
unit)*  and  pollutant  type. 

Nlaar  daaa  ant  have  to  eater  the  amber  of  engines  par  aircraft 

type- 
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The  short  tins  source  connecting  the  taxiway  with  the  start  of  the 
runway  takeoff  roil  is  Intended  to  be  assigned  the  emissions  associated  with 
the  final  aircraft  msneuvers,  generally  including  a  90°  or  180*  turn  plus  a 
pause  prior  to  takeoff  clearance* 

The  line  sources  associated  with  the  runway  queue  would  be  assigned  all 
extra  emissions  associated  with  delays  in  takeoff  due  to  queuing*  A  recom¬ 
mended  procedure  for  dealing  with  runway  queues  is  described  in  Sec*  9*6, 

One  additional  refinement  is  recommended  for  inclusion  in  treating 
taxiway  segments*  A  minimum  time  spent  on  any  segment  should  be  defined 
(nominally  30  sec)*  The  time  spent  on  a  particular  segment  should  be  taken  as 
the  maximum  of  this  minimum  time  and  the  time  computed  using  taxiing  speed  and 
segment  length*  This  procedure  will  avoid  any  unrealistically  short  times  and 
allow  for  aircraft  turning* 

Generally,  the  user  should  be  advised  to  keep  the  number  of  Individual 
taxiway  line  segments  to  a  minimum  to  reduce  computational  and  input 
requirements* 


9.6  AIRCRAFT  QUEUING 

In  view  of  the  Importance  of  aircraft  runway  queuing,  or  aircraft 
queuing  in  other  locations  for  that  matter,  it  is  worthwhile  devoting  some 
detailed  attention  to  it*  To  cover  the  widest  possible  number  of  situations, 
it  is  recommended  the  user  be  given  three  options: 

1*  Ignore  queuing  altogether* 

2*  Include  queuing  as  treated  by  queuing  algorithms 
Incorporated  in  the  computer  program* 

3*  Include  queuing,  but  user  supplies  his  own  delay  times* 

If  option  2  or  3  above  is  selected,  the  user  should  then  be  required  to 
provide  the  following  information: 

#  For  option  2  or  3,  glee  the  Identification  number  of  one  or 
two  taxiway  segments  that  will  contain  the  aircraft  queue 
for  each  runway*  Taxiway  segments  must  always  be  straight- 
line  segments*  Because  of  taxiway  configuration  and 
anticipated  queue  lengths,  it  may  be  necessary  to  use  two 
taxiway  segments  to  contain  the  entire  queue*  (Queuing  of  i 

greater  complexity  could  be  handled,  but  the  queuing 
algorithm  used  by  the  computer  code  becomes  cumbersome.) 


e  For  option  3,  give  the  average  delay  time  for  each  aircraft 
entering  the  queue*  If  the  user  selects  this  option,  it 
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will  be  utuaid  that  the  utrt  Million*  due  to  queuing  ere 
uniformly  distributed  over  the  taxiway  aegaents  designated 
by  the  user  for  the  queue* 

If  the  user  selects  opr  ton  2,  the  computer  algorithm  can  compute  the 
queue  length  and  compare  it  with  the  length  of  the  taxiway  segments  designated 
by  the  user.  If  the  actual  queue  length  differs  significantly  from  the  length 
of  the  taxiway  segment,  then  the  code  can  define  new  line  sources  collnear 
with  the  taxlwsy  segments  to  represent  the  queue* 

The  queuing  algorithm  Itself  can  be  based  upon  the  following  hypo¬ 
thesis,  which  Is  itself  based  on  dasslcel  queuing  theory*  The  propoeed 
hypothesis  says  that  all  aircraft  must  stop  and  wait  to  be  served  (by  the 
runway)  and  that  the  average  number  of  aircraft  waiting  to  be  served  is  given 
by: 


where: 

N  -  the  number  of  departing  aircraft  (N^)  waiting  to  be  served 
plus  the  number  of  arriving  aircraft  (N#)  waiting  to  be 
served, 

V  •  the  aircraft  traffic  volume,  l*e*,  the  number  of  aircraft 
using  the  runway  for  both  arrivals  (nfl)  and  departures 
(n^)  per  hour,  and 

C  -  the  runway  capacity  (the  maximum  posaible  number  of 
arriving  and  departing  aircraft  that  can  be  served  per 
hour)* 


Note  that  this  is  a  steady-state  hypothesis*  It  Is  aseumsd  that  the  aircraft 
traffic  volume,  V,  persists  for  more  than  one  hour.  The  actual  number  of 
aircraft  in  the  queue  will,  of  course,  fluctuate  from  time  to  time*  N  Is  to 
be  regarded  aa  an  average  value  over  the  hour* 


Since  the  queue  contains  both  arriving  (M#)  and  departing  (H^)  aircraft 
in  ganeral,  and  since  only  the  departure  queue  Is  of  concern,  a  way  must  be 
found  to  determine  the  number  of  aircraft  In  the  daparture  queue  (H^) 
separately.  Here  It  Is  assumed  that  Nj  is  proportional  to  the  total  number  of 
departing  aircraft  (n^)*  Hence: 


or 


(2) 
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n. 

It  follows  that  »  —  H  or,  using  Eq.  1: 


(3) 


If  all  of  tha  aircraft  In  tha  quaua  wars  of  tha  sane  typa ,  tha  length 
of  tha  departure  quaua  would  ha  given  by: 

L  -  Mtd  (4) 

where  D  ■  distance  occupied  by  one  aircraft.  If  several  aircraft  types  are 
Involved,  than  tha  nunber  of  aircraft  of  each  typa  (Hj)  nust  be  estlnatad. 
Then  tha  length  of  the  departure  queue  Is  given  by: 

L  -  ID^  (5) 


Tha  distance  occupied  by  each  aircraft  type  (Dj)  can  be  roughly  estlnated  to 
be  2  to  2.5  aircraft  lengtha.  The  nunber  of  aircraft  of  each  type  In  the 
queue  (Hj)  can  be  roughly  eatlnated  fron  the  ratio  of  the  nunber  of  that  type 
ualng  the  runway  for  departure  (n^)  to  the  total  nunber  of  aircraft  using  that 
runway  for  departure  (n^).  That  la: 


(6) 


Finally,  the  average  extra  tine  In  ninutes  per  departing  aircraft  due 
to  queuing,  not  Including  the  taxi  tine  in  the  absence  of  queuing.  Is  the  saae 
for  all  aircraft  types  and  Is  given  by: 


T 


(7) 

I 


The  one  paraaeter  left  to  be  estlnated  Is  the  runway  capacity,  C.  A 
detailed  evaluation  of  C  based  upon  all  the  available  data  la  beyond  the  scope 
of  the  present  effort  to  outline  the  new  eoaputatlonal  aye ten.  However,  a 
prelinlnary  review  of  sons  data  suggests  that  C  is  in  tha  neighborhood  of 
46.  For  exaaple,  If  the  hourly  rate  of  departures  (nj)  and  arrivals  (na)  on  a 
runway  ware  30  and  12,  respectively,  then  fron  Eq.  1,  the  average  nunber  of 
aircraft  awaiting  runway  service  would  be  seven.  Of  these,  five  would  be  In 
the  departure  queue  according  to  Eq.  3.  Consequently,  each  aircraft  that 
eaters  the  departure  queue  would  have  to  wait  6.25  aln  according  to  8q.  7. 


It  Is  worth  noting  that  if  a  eonaarclal  runway  warn  used  exclusively 
for  departures.  It  could  accoanodats  about  60  aircraft  per  hour,  tf  it  ware 
aaad  exclusively  for  arrivals.  It  could  accoanodats  about  40  aircraft  per 
hour,  allowing  for  a  separation  tins  of  1.5  aln.  Hence,  the  value  of  C  is 
expected  to  lie  between  the  values  of  40  and  60  aircraft  per  hour.  Generally, 
runways  are  not  used  exclusively  for  arrlvele  or  departures. 
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It  Ls  worth  stressing  here  that  N  and  T  are  average  values  only. 
Furthermore,  the  proposed  queuing  algorithm  ls  only  a  hypothesis,  and  the 
suggested  value  of  C  Is  only  an  estimate.  It  ls  strongly  recommended  that  the 
hypothesis  be  tested  and  that  the  value  of  C  (which  would  be  stored  In  the 
permanent  aircraft  data  file)  be  evaluated  using  data  for  various  aircraft 
facilities. 

It  should  be  noted  that  the  AQAM  did  not  use  any  queuing  algorithm  and 
the  AVAP  model  used  a  queuing  algorithm  based  on  O'Hare  Airport  data  only. 
Unfortunately,  very  little  information  exists  on  the  basis  for  the  form  of  the 
AVAP  queuing  algorithm.  Essentially,  the  latter  algorithm  has  the  form: 

T  -  (j*  ♦  1)  (8) 


where: 


T  -  delay  time  In  minutes  per  aircraft,  and 

nd  *  number  of  departures. 

For  the  example  given  above,  if  nd  *  30  aircraft  per  hour,  T  •  11  min.  The 
AVAP  queuing  algorithm  has  a  simple  linear  dependence  on  the  aircraft  traffic 
volume,  whereas  In  the  algorithm  proposed  here  (see  Eq.  1),  V  appears  In  both 
the  numerator  and  denominator  In  such  a  way  that  as  V  approaches  the  runway 
capacity,  C,  the  number  of  aircraft  In  the  queue,  and  therefore  the  delay 
time,  rapidly  Increases.  This  la  as  expected,  since  once  the  traffic  flow 
exceeds  the  runway  capacity,  the  queue  length  simply  Increases  with  time,  and 
there  ls  no  longer  any  meaning  to  the  Idea  of  an  average  (steady-state)  length 
of  the  queue.  When  V  ls  relatively  small  compared  with  C,  the  delay  time  Is 
smaller  for  the  proposed  algorithm  than  for  the  original  AVAP  algorithm. 

The  proposed  queuing  algorithm  can  be  coded  In  a  subroutine  called 
"QUEUE,"  which  would  function  as  follows.  QUEUE  is  called  from  the  Aircraft- 
Source-Emission-Inventory  code  whenever  the  queue  flag  1QUEPL  *1.  If  1QUEFL 
*  0,  queuing  is  Ignored.  If  1QUEFL  *  2,  the  user  must  input  the  queuing  times 
QUET  (IE)  for  each  runway  (IE).  All  aircraft  types  suffer  the  same  time  delay 
caused  by  the  queue  at  a  runway* 

QUEUE  performs  the  following  functions  for  each  runway: 

l*  It  computes  the  number  of  aircraft  In  the  queue* 

2.  It  apportions  this  number  amongst  the  different  types  of 
aircraft*  The  number  of  aircraft  of  a  given  type  in  the 
queue  ls  likely  to  be  a  noninteger* 

3*  It  estimates  the  physical  length  of  the  queue  (QL). 
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4.  It  compares  QL  with  th«  length  of  the  first  taxiway 
aagaant  (TSL)  to  which  the  queue  la  assigned  hy  the 
user.  If  (JL  -  TSL  *  20%,  e  variable  celled  "FLAG"  Is  set 
■  1,  end  the  extra  ealsslons  due  to  queuing  ere  uniformly 
distributed  over  the  first  texlwey  segment  used  for 
queuing.  If  QL  <  TSL  -  20%,  FLAG  is  set  •  2  end  a  new 
line  source,  collnear  with  the  texlwey  segment.  Is 
defined.  The  extra  emlsslone  due  to  queuing  are  then 
assigned  to  this  new  line  source.  If  QL  >  TSL  +  20%,  FLAG 
is  set  ■  3,  end  the  queuing  emissions  are  apportioned  to 
the  two  taxiway  segments  designated  by  the  user  for  the 
queue.  The  first  texiway  segment  is  assigned  emissions 
corresponding  to  the  fraction  T8L/QL  of  the  total  delay 
time.  The  remainder,  namely  that  sesocleted  with  the 
fraction  (1  -  TSL/QL),  Is  assigned  to  s  new  line  source  of 
length  QL  -  TSL  that  is  collnear  with  the  second  texlwey 
segment  used  for  queuing. 

3.  Finally,  the  subroutine  QUBUS  returns  the  value  of  FLAG 
along  wltH  the  value  of  the  extra  queuing  time  per 
aircraft  "TQUE"  (for  FLAG  ■  1  or  2)  and  the  value  of 
"TQUEP"  if  FLAG  -  3.  "TQUE"  and  “TQUEP"  ere  the  names  of 
variables  used  to  store  the  delay  times  assigned  to  the 
first  and.  If  used,  the  secood  texiway  segments  (or  line 
sources)  designated  by  the  user  to  contain  the  queue. 

If  a  new  line  source  Is  defined,  Its  second  end-point  coordinates  ere 
Identical  to  those  of  the  texiway  segment  with  Vhleh  it  Is  collnear.  Its 
first  end-point  coordinates  ere  given  by: 

XQl  *  xl  +  <X2  '  Xl>»  <9> 

»«i  -  ’i  ♦  I5Ss2‘  (,2  -  v 

where  (XpTj)  end  (XjtTj)  refer  to  the  first  end  secood  end  points  of  tbs 
texlwey  segment  collnear  with  the  now  line  source,  respectively. 

A  detailed  flowchart  of  subroutine  QOIUS  Is  given  In  Ref.  30. 


9.7  AIRCRAFT  MRVICB-VBN1CLR  EMISSIONS 

To  simplify  the  computation  of  emissions  from  aircraft  service 
vehicles,  a  single  emission  fnctor  should  bn  defined  that  represents  the 
emissions  from  ell  service  vehicles  need  to  service  a  single  aircraft  of  a 
given  type.  For  a  civilian  airport,  the  user  la  asked  to  enter  the  number  of 


aircraft  of  aach  typa  that  *111  bt  serviced  during  tha  hour  of  Intareat  at 
each  aircraft  gate  area.  For  a  military  air  base,  tha  uaer  la  aaked  to  enter 
the  number  of  aircraft  of  each  typa  that  *111  require  abut  down  servicing 
during  tha  hour  of  interest  at  aach  aircraft  parking  area  and  the  number  that 
will  require  startup  service  during  tha  same  hour  (not  necessarily  the  same 
aircraft)*  The  difference  between  civilian  and  military  operations  is  that 
commercial  aircraft  are  routinely  turned  around  with  a  single ,  more  or  leas 
continuous  servicing  operation  lasting  30-60  min#  whereas  mtlltary  aircraft 
may  be  serviced  In  two  separate  operations,  depending  on  whether  the  aircraft 
are  transient  or  permanently  assigned  to  a  base*  Some  of  the  military 
aircraft  may  undergo  only  one  turnaround,  while  others  may  be  turned  around 
several  times  in  a  given  day*  Hence,  whereas  one  service-vehicle  emission 
factor  for  each  aircraft  type  will  suffice  at  a  commercial  airport,  separate 
service-vehicle  emission  factors  for  shutdown  and  startup  operations  are 
required  for  each  aircraft  type  at  a  military  air  base* 

The  military  service-vehicle  emission  factors  originally  used  in  AQAM 
have  recently  been  updated  by  the  staff  at  Tyndall  AFB*  These  updated  factors 
should  be  implemented  in  the  new  computational  system  in  the  permanent 
aircraft  data  file* 

The  civilian  service-vehicle  emission  factors  originally  used  in  AVAP 
were  based  on  a  study  at  O'Hare  Airport  (circa  1972)  and  have  not  been 
updated*  The  types  of  service  vehicles  used  at  O'Hare  Airport  are  listed  In 
Table  6  along  with  the  total  service  times  (sum  of  service  times  for  all 
service  vehicles  of  the  same  type)  per  aircraft  type.  This  information  was 
compiled  from  quest lonnalres  sent  to  the  various  airlines*  The  emission  rates 
for  each  of  these  vehicle  types  were  determined  by  comparing  vehicle 
characteristics  with  those  of  vehicles  for  which  EPA  emission-factor  data  had 
been  published*  They  should  be  updated  before  being  incorporated  in  the 
permanent  aircraft  data  file  of  the  new  computational  system* 

A  procedure  for  computing  aircraft  service-vehicle  emissions  in  airport 
gate  areas  is  outlined  below*  (A  similar  procedure  could  be  used  for  startup 
and  shutdown  servicing  operations  at  a  military  facility*)  The  objective,  of 
course,  is  to  compute  total  emissions  per  pollutant  type  per  aircraft 
serviced*  Since  the  total  time  required  to  perform  all  services  ranges  from 
30  min  to  60  min,  service-vehicle  emissions  must  be  allocated  to  the  proper 
model  hour**  That  Is,  all  aircraft  arriving  at  the  airport  during  a  given 
hour  may  not  all  be  serviced  during  that  hour*  However,  other  aircraft  that 
arrived  during  a  previous  hour  nay  be  serviced  during  the  model  hour*  The 
uaer  should  have  the  option  to  specify  any  number  that  he  wlehee  for  aircraft 
bslng  serviced  during  the  model  hour*  That  number  dose  not  have  to  beer  any 
relationship  to  the  numbar  of  aircraft  arriving  or  dapartlng  during  tha  model 
hour* 


*Hour  for  which  emission  calculation  la  bslng  parformad 
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Table  6  Total  Minutes  of  Service-Vehicle  Operation 
Allocated  to  Servicing  Each  Aircraft  Type 


Aircraft  Type 

1* 

2* 

3* 

4* 

5 

6 

7 

8 

9 

10 

Vehicle  Type 

727 

DC  9 

737 

C5 

BAC 

YS 

B9 

PH 

TO 

GA 

1  Tractor 

66 

48 

83 

55 

50 

50 

0 

0 

0 

0 

2  Belt  Loader 

28 

15 

30 

0 

25 

25 

0 

0 

0 

0 

3  Container  Loader 

6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4  Cabin  Service  Truck 

12 

0 

15 

0 

0 

0 

0 

0 

0 

0 

5  Lavatory  Truck 

15 

15 

15 

10 

10 

10 

5 

5 

5 

0 

6  Water  Truck 

0 

10 

0 

10 

10 

10 

5 

5 

5 

0 

7  Food  Truck 

17 

17 

20 

10 

10 

10 

0 

0 

0 

0 

8  Fuel  Truck 

20 

15 

15 

10 

20 

20 

10 

10 

10 

0 

9  Tow  Tractor 

10 

5 

5 

5 

5 

5 

0 

0 

0 

0 

10  Conditioner 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

11  Aircraft  Starting 

Unit,  Transporting 
and  Diesel  Engines 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

12  Ground  Power  Unit, 

Transporting  and 
Gasoline  Engines 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

13  Ground  Power  Unit, 

Diesel  Engine 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

14  Transporter 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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ALGORITHM:* 

1.  USER  INPUTS  TOTAL  TIMES  REQUIRED  TO  SERVICE  1  TYPE  IAC 
AIRCRAFT  WITH  TYPE  ISV  SERVICE  VEHICLES 

-  SVTIME  (ISV, IAC)  (min) 

2.  USER  INPUTS  EMISSION  FACTOR  FOR  SERVICE  VEHICLE  TYPE  ISV 
AND  EACH  POLLUTANT  IP 

-  SVEMFT  (IP, ISV)  (g/min) 

3.  COMPUTE  TOTAL  EMISSIONS  OF  POLLUTANT  IP  FROM  ALL  SERVICE 
VEHICLES  FOR  ONE  AIRCRAFT  OF  TYPE  IAC 

SVEF  (IP, IAC)  -  Y.  SVEMFT(  IP,ISV)*SVTIME(  ISV, IAC) 

ISV 

4.  PRINT  OUT  SVEF  (IP, IAC)  OR  STORE  ON  FILE.  (This  data  must 
be  entered  Into  Pile  1,  Permanent  Aircraft  Data*) 


9.8  ACCESS  VEHICLES 

Access  vehicles  Include  all  ground  vehicles  transporting  personnel  and 
equipment  Into,  out  of,  and  around  an  aircraft  facility.  At  a  commercial 
facility  this  traffic  will  consist  mostly  of  passenger  vehicles.  At  a 
military  base  a  larger  fraction  of  the  vehicles  Is  likely  to  be  Involved  In 
transporting  personnel  Into,  out  of,  and  around  the  base. 

In  order  to  compute  the  emission  rates  from  each  roadway  segment  or 
parking  lot,  It  la  first  necessary  to  define  a  scenario  that  describes  the 
nature  of  the  vehicle  operations  on  that  segment  or  lot  and  the  vehicle  mix 
Involved.  Next,  the  emission  factors  for  each  scenario  must  be  computed. 
Once  the  applicable  scenarios  for  each  source  are  defined  and  the  emission 
factors  computed,  lt  is  straightforward  to  compute  the  emission  rates  for  each 
source  by  multiplying  the  emission  factors  by  the  levels  of  vehicle  activity. 
Unfortunately,  the  process  of  co^suting  the  emission  factors  Is  extremely 
tedious.  Portunately,  existing  moblle-source-emlsslon  computation  routines 


♦User  may  accept  default  values  for  both  SVTIME  (ISV, IAC)  and  SVBMPT  (IP,ISV) 
or  substitute  his  own  Input  values*  Values  of  SVTIME  (ISV, I AC)  must  be  based 
on  observations  at  airports  or  air  bases.  Values  of  SVEMPT  (IP, ISV)  can  be 
estimated  based  on  emission  factors  for  slses  of  similarly  fueled  vehicles 
that  are  published  In  Ref.  32.  Por  military  air  bases,  it  is  appropriate  to 
define  two  service-vehicle  emission  factors  per  aircraft  types  a  factor  for 
inbound  aircraft  of  type  IAC  and  pollutant  IP,  SVEPI  (IP,  I  AC),  and  a 
corresponding  factor  for  outbound  aircraft,  SVEPO  (IP, IAC). 
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are  now  readily  available*  The  procedure  recommended  here  Incorporates  the 
use  of  M0BILE2 J  or.  equivalently,  output  f roe  that  or  similar  computer  codes* 

A  set  of  emission  factors  must  be  computed  by  MOBILE233  for  each 
emission  scenario*  An  emission  scenario  must  be  assigned  to  each  physical 
source  (roadway  segment  or  parking  lot)  so  that  the  appropriate  set  of 
emission  factors  can  be  applied  to  each  source*  Depending  upon  the  variety  of 
combinations  of  conditions  encountered,  it  may  be  necessary  to  define  several 
scenarios  to  accommodate  all  of  the  physical  sources  at  an  aircraft  facility* 

Generally,  three  types  of  trips  will  be  encountered  at  a  civilian 
airport  during  the  course  of  a  one-hour  period* 

1*  A  vehicle  enters  the  airport  complex9  travels  to  a 
terminal  building  along  several  roadway  segments9  parks 
temporarily  near  the  building  (often  with  the  engine  left 
Idling) 9  then  departs  from  the  airport*  The  vehicle  is  in 
the  hot-stabilised  running  condition  for  the  entire  trip* 

2*  A  vehicle  enters  the  airport  complex  and  travels  to  a 
parking  lot  where  it  remains  for  at  least  one  hour* 
Emissions  result  from  driving  plus  hot-soak  evaporative  HC 
losses  from  the  carburetor.  (A  small  additional  source  of 
evaporative  HC  losses  is  due  to  diurnal  ambient  tempera¬ 
ture  changes  that  result  in  fuel  tank  evaporation  losses*) 

During  the  one-hour  period  the  hot-soak  loss  predominates* 

3*  A  vehicle  leaves  the  parking  lot  and  then  the  airport 
complex*  (All  operations  are  assumed  to  take  place  in  the 
cold-start  running  condition.) 

In  addition,  during  the  hour  of  interest 9  a  number  of  vehicles  are  expected  to 
remain  parked  in  the  lots*  Only  fuel-tank  evaporative  HC  losses  are 
associated  with  these  vehicles. 

Depending  on  the  design  of  the  access-vehicle  roadway  system,  the  same 
roadway  segment  may  serve  both  the  parking  lots  and  the  through  traffic* 
Hence,  vehicles  in  both  the  cold-start  and  hot-stabilised  running  conditions 
will  occupy  this  segment  simultaneously*  Such  a  mix  can  be  easily  handled 
through  the  scenario  parameters  required  as  input  to  M0BILE2* 

The  four  vehicle  operating  modes  of  Interest  are: 

l*  Cold-start  condition  —  vehicles  leaving  parking  lots* 

2*  Hot-start  condition  —  following  short  shutdown  interval* 


3*  Hot-stabilised  condition  —  all  vehicles  entering  airport* 


4.  Idle  at  hot -stabilized  condition  —  all  vehicle*  entering 
airport. 

The  following  types  of  emi salon  factors  must  be  determined  from  M0BILE2 
runs  for  each  scenario  type: 

1.  Driving  emissions  (for  each  pollutant  type,  for  each 
vehicle  type)  (g/VMT)*  CO,NMHC,NOX. 

2.  Idle  emissions  (for  each  pollutant  type,  for  each  vehicle 
type)  (g/mln)  CO,NMHC,NOX. 

3.  Hot -soak  evaporative  losses  from  carburetor  for  each  hot- 
soak  period  (g/hot  soak)  NMHC. 

4.  Fuel-tank-breathing  evaporative  losses  (g/hour)  WWC.** 


9.8.1  Scenario  Type  Definition 

Each  scenario  type  la  defined  using  the  following  quantities: 

1.  Region  of  country. 

2.  Calendar  year. 

3.  Average  vehicle  speed  over  source  (roadway  segment  or 
parking-lot  lanes);  either  a  single  value  for  all  vehicle 
types  or  up  to  eight  values  for  separate  vehicle  types. 

4.  Ambient  temperature. 

5.  Percentage  of  VMT  on  source  in  cold-start  running 
condition. 

6.  Percentage  of  VMT  on  source  in  hot-start  running  condition 
(M0BILE2  assumes  X  cold  start  +  t  hot  start  +  X  hot 
stabilized  -  100X). 

7.  Vehicle  mix  (percentage  of  VMT  on  source  by  vehicle  type). 

8.  Correction  factors  to  basic  emission  factors. 


♦VMT  -  vehicle  miles  traveled  (number  of  vehicles  *  length  of  roadway  segment 
In  miles). 

♦♦Items  1  and  4  are  required  for  each  vehicle  type  and  can  be  determined  from 
M0BILE2  output  using  a  special  procedure  described  In  ftec*  9.8,2. 
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In  addition  to  the  above  quantities,  it  is  also  necessary  to  specify 
the  number  of  trips  per  day  (actually  equivalent  to  the  number  of  hot-soak 
periods  per  day)  and  the  average  dally  mileage  (see  Sec*  9.8.2). 


9.8.2  Procedure  for  Determining  Evaporative-Loss  Emission  Factors 

H0BILE2  requires  the  input  of  tiro  special  quantities  for  each 
scenario:  the  number  of  trips  per  day  and  the  average  dally  mileage.  These 

two  quantities  are  not  used  to  compute  the  driving  emission  factors.  They  are 
only  used  by  H0BILE2  to  convert  the  units  of  the  fuel-tank  HC  breathing  loss 
from  grams  per  day  to  grams  per  mile  and  the  units  of  the  carburetor  hot-soak 
HC  loss  from  grams  per  hot  soak  to  grams  per  mile.  Unfortunately,  the  latter 
units  are  not  the  units  needed  by  the  new  computational  system  to  compute 
evaporative  losses  from  access  vehicles.  The  required  units  are: 

For  carburetor  hot -soak  loss:  grams  per  hot  soak 


For  fuel  tank  evaporative  loss:  grams  per  hour 

In  other  words,  M0BILE2  computes  the  quantity: 

.  a  Mb 
E  - 

where: 

E  -  evaporative  emissions  (g/ml), 

a  -  diurnal  average  fuel-tank  breathing  loss  (g/24  hr), 

N  -  number  of  trips  or  hot  soaks  per  day, 
b  -  hot-soak  carburetor  evaporative  loss  (g/hot  soak),  and 
L  -  total  number  of  miles  traveled  per  day. 


(10) 


What  Is  needed  for  access-vehicle  emission  computations  are  the 
separate  quantities  a  and  b.  These  two  quantities  can  be  obtained  from  the 
results  of  two  separate  H0BILK2  runs  as  follows.  For  the  first  run,  set  the 
number  of  trips,  H  -  0,  and  the  number  of  miles  traveled  per  day,  L  *  1.  Then 
the  output  of  H0BILE2  will  be  Ej  •  a.  For  the  second  run,  set  N  «  l  and  L  » 
1.  Then  the  output  of  H0BILR2  will  be  E2  •  a  ♦  b.  Given  these  two  outputs, 
E|  and  Ej*  the  separate  quantities  a  and  b  are  computed  as: 

a  -  E|  and  b  •  Ej  -  Ej  (11) 

where  a  is  fuel-tank  breathing  loss  in  grams  per  day,  and  b  is  the  hot-soak 
loss  in  grams  per  hot  soak.  To  get  the  breathing  loss  in  units  of  grams  per 
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hour,  iltply  divide  a  by  24.  A  flowchart  illuatrating  a  computational 
procedure  for  determining  a  and  b  la  given  in  Pig*  10. 


9.8.3  Access-Vehicle  Physical-Source  Types 

Two  typea  of  physical  sourcea  are  used  for  access  vehicles:  (1) 

straight-line  roadway  segments  (may  contain  several  parallel  lanes)  and  (2) 
square-area  parking  lota  (several  may  be  required  to  cover  irregularly  shaped 
parking  lots).  Two  types  of  information  are  required  for  each  source  type: 
(1)  physical  description  of  source  and  (2)  traffic  information. 


Physical  Description 

1.  Source  type  (1.  Roadway,  2.  Parking  lot) 


2.  If  type  1:  elevation  >  0,  above  grade 

■  0,  at  grade 
<  0,  below  grade 
overall  roadway  width 
number  of  lanes 
lane  width 

end-point  coordinates  (center  of  roadway 

cross-section  at  each 
end  of  segment) 


3.  If  type  2:  number  of  levels 

overall  height  (•  0,  if  only  ground  level) 
length  of  side  of  square  area 
coordinates  of  center 


Traffic  Information  (for  a  one-hour  computation) 

1.  Scenario  type  (a  number  identifying  the  applicable 
scenario  emission  factors  from  M0BILE2) 

2.  Por  source  type  1  (roadway  segment): 

-  Number  of  vehicles /hour  of  each  type 

-  Number  of  minutes  of  extra  idle  tlme/vehlcle/vehlcle 
type 
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For  source  type  2  (perking  lot) 


-  Nuaber  of  vehicles  of  each  type  entering  lot  during  hour 
(assuae  all  are  In  hot-stablllsed  running  condition)* 
Each  such  vehicle  will  have  driving  ealsslons.  Idle 
ealsslons,  hot -soak  carburetor  evaporative  HC  ealsslons. 

-  Nuaber  of  vehicles  of  each  type  leaving  lot  during  hour 
(assuae  all  are  In  cold-start  running  condition).  Each 
such  vehicle  will  have  cold-start  driving  ealsslons  and 
cold-start  Idle  ealsslons.  Unfortunately,  M0BILE2  does 
not  produce  cold-start  Idle-ealsslon  factors.  Hence, 
extra  driving  esdsslons  aay  be  Introduced  to  coapensate. 

-  Nuaber  of  vehicles  reaslnlng  parked  during  hour  (fuel- 
tank  evaporative  HC  ealsslons  only). 

-  Average  Inbound  alleage  (aaasured  along  path  followed  by 
average  vehicle). 

-  Average  outbound  alleage  (augmented  by  additional 
alleage  to  account  for  extra  Idle  tlae  if  any). 


9.8.4  Outline  of  the  Coaputatlon  of  Access-Vehicle  Ealaalon  Hates 

Figure  11  outlines  the  overall  procedure  for  coaputlng  ealsslon  factors 
using  H0BILE2  and  source  ealsslon  rates.  Guidelines  for  selecting  Input 
paraaeter  values  for  running  M0BILE2  are  given  In  Tables  7a  and  7b.  The 
ealselon-factor  coaputatlonal  procedure  Is  outlined  In  eoaewhat  greater  detail 
In  Fig.  10.  A  detailed  flowchart  showing  the  coaputatlon  of  access-vehicle 
source  ealsslon  rates  Is  given  In  kef.  30. 


9.9  HYDROCARBON  EVAPORATIVE  LOSSES  DUE  TO  FUEL  STORAGE  AND  HANDLING 

As  Indicated  eerller,  the  evaporative  HC  esdaslons  tend  to  be  a  saall 
fraction  of  the  total  HC  ealsslons  due  to  coabustlon  and  evaporation. 
Nevertheless,  In  soae  cases,  the  user  aay  wlah  to  include  an  astlaata  of 
ealsslons  froa  this  source  class,  which  Includes  fuel  storage  and  handling. 
The  original  AQAH  Included  a  detailed  treataent  of  evaporative  HC  ealsslons 
using  essentially  tha  procedures  outlined  below.  The  AVAP  aodel  Ignored  all 
evaporative  losses. 
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n*.  11  Macro  Flowchart  of  Access-Vehicle 
Source-Balsa  lon-Rate  Coaput at  Ion 


Thera  are  four  types  of  evaporative  losses  at  a  typical  aircraft 
facility: 

1 .  Breathing  loases  due  to  diurnal  taaperature  changes  that 
result  In  emissions  from  fixed-roof  storage  tanks  as  well 
as  tank  trucks  and  vehicle  fuel  tanks  that  do  not  have 
vapor-control  systems. 

2.  (forking  losses  due  to  displacement  of  vapors  during 
filling  operations  where  vapor-control  systems  are  not 
used. 


Table  7a  Input  Data  for  M0B1LE2  Computational 
Control-Section  Inputs* 


Recommended 
Variable  Input  Value 


Meaning 


SPDFLG 

VMFLAC 

ALTPLC 

IMFLG 

ALHFLC 

l  FORM 

PRTFLG 

ICKVFG 

IDLFLC 


0  Indicates  user  mill  use  a  single  average  speed  for 

all  vehlclea  on  roadway  section 

1  Indicates  user  will  enter  own  VMT  mix  In  PARAMETER 

SECTION*  This  will  enable  user  to  define  aa  many 
distinct  vehicle  mixes  aa  are  needed  to  describe  all 
access-vehicle  roadways* 

0  MO BILE 2  emission  factors  will  be  used* 

0  No  Inspection/maintenance  credit* 

l  User  will  utilise  additional  correction  factors  (see 

MOBILE 2  Manual). 

0  or  2  Indicates  user  wants  numerical  output  for  subsequent 

use  by  other  computer  codes* 

0  Indicates  user  wants  output  of  all  available 

pollutant  (THC,CO#NOX)  emission  factors* 

3  User  supplies  number  of  trips  per  day  and  mileage  per 

day  to  compute  evaporative  HC  emission  factor*  User 
wants  separate  output  of  these  emission  factors* 

1  User  wants  output  of  idle  emission  factor. 

User  wants  output  of  NMHC  in  lieu  of  total 
hydrocarbon  (THC)  emission  factor* 


NMHFLG 


l 
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Table  7b  Input  Data  for  M0BILE2  Computational 
Parameter-Section  Inputs  and  Flags* 


Record  Variable 


Meaning 


SCENARIO 


One  record  for  each  scenario* 


IREJN 

CY 

SPO 


TAMB 


Flag  Identifying  geographic  region* 

Last  two  digits  of  calendar  year. 

Single  eve rage  route  (roadway-segment)  speed.  If 
the  fraction  of  idle  time  to  total  time  on  roadway 
segment  la  substantially  different  from  that  used 
in  the  M0BILE2  assumptions,  then  the  excessive  Idle 
time  should  be  treated  separately*  This  extra  Idle 
time  may  be  combined  with  the  Idle  eadsslon  factors 
to  compute  the  additional  idle  emissions  on  the 
roadway  segment. 

Ambient  temperature. 


PCCN 

PCHC 

PCCC 


VMT  MIX 


Percentage  of  VMT  In  cold-start  mode  by  noncatalyt- 
lcally  equipped  vehicles. 

Percentage  of  VMT  In  hot-start  node  by  catalyt- 
lcally  equipped  vehicles* 

Percentage  of  VMT  In  cold-start  node  by  catalyt- 
lcally  equipped  vehicles. 

If  VMFLAG  ■  1,  user  must  supply  fraction  of  total 
VMT  traveled  by  each  vehicle  type.  This  feature 
allows  user  to  simulate  a  broad  range  of  mixtures 
of  vehicles  on  public  as  wall  as  llnltsd-access 
roadways  used  by  cargo  and  service  vehicles. 


ADDITIONAL  LIGHT-DUTY  If  ALHFLG  -  1  or  2,  us or  must  supply  approprlats 

GASOLINE-POWBRED-VEHICLB  corrsctlon  factors  (ses  M0BILB2  Manual). 
CORRECTION  FACTORS 
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3.  Standing  storage  loaaaa  due  to  poor  eeala  on  floating-roof 
storage  tanka* 

4.  Spillage* 

The  firat  three  types  of  loaaaa  can  be  estlaated  with  the  uae  of  well- 
known  equations  that  require  Inforaatlon  about  the  tanka  aa  well  aa  sows 
aeteorologlcal  Inforaatlon*  The  required  foruulae  are: 


Standing  Storage  Losses 

SL  •  155  “i"*'5  (t5X=t)°','!'’ciSc3 


where: 


(12) 


SL  •  standing  loss  (lb/day), 

W  -  liquid  density  (lb/gal), 

Kj  ■  tank  construction  factor, 

■  0.045  for  welded  tanks, 

■  0.13  for  riveted  tanks, 

D  -  tank  disaster  (ft),* 

?  •  true  vapor  pressure  of  the  bulk  liquid  at  its  average 
storage  teaperature  (pala), 

Vtf  •  average  wind  speed  (aph), 

Cj  “  tank  seal  factor  (for  alapllclty  adopt  only  one  value  for 
each  tank  fara), 

*  1.0  for  tight-fitting,  aodern  seals, 

•  1.33  for  loose-fitting  eeala  (typical  of  those  built 
before  1942), 

Cj  ■  fuel  factor, 

-  1.00  for  gasoline, 

■  0.96  for  naphtha  (JM), 

■  0*83  for  karoeeae. 


•If  D  >  150  ft,  aaa  D<1»0)0*5  in  lieu  of  D1*3. 


t 


19 


f 
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•  0.79  for  distillate  oil,  and 

•  0.75  for  crude  oil, 

C3  •  tank-color  factor  (for  simplicity  use  0.95  for  all 
tanks), 

*1.0  for  light  gray  or  aluminum  finish, 

-  0.9  for  white, 


Vapor  Pressure 

P(T)  -  exp  (o  -  0/T)  (13) 

where: 

P(T)  -  vapor  pressure  (psia)  as  a  function  of  ambient  tem¬ 
perature  T, 

T  -  ambient  temperature  (*A),  and 
a  and  0  -  parameters  that  depend  on  the  fuel  type. 


Breathing  Losses 


BL 


wk  d1*73 

365  2U 


( _ L_\o- 

\  14.7  -  P/ 


68  0.51 
H 


at°*5c4c5 


whe  re : 


(14) 


BL  -  breathing  loss  (lb/day), 

*2  *  liquid-dependent  factor, 

*  0.014  for  crude  oil, 

*  0.019  for  distillate  oil, 

*  0.020  for  kerosene, 

-  0.023  for  JP4, 

-  0.024  for  gasoline, 

H  *  average  vapor  space  height  Ut)  (select  one  value  for 
entire  baee). 


AT  *  average  diurnal  temperature  variation  (*F), 
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C4  •  finish  factor,  which  varies  f  row  1  to  1.S8  (for 
slapllclty  choose  one  representative  value  for  entire 
base),  and 

C5  ■  adjustaant  factor  for  tanka  <20  ft  in  dieaeter. 

For  tank  trucks,  use: 

4  (1  -  C  )TC 

- - y-5 -  (15) 

*D 


where: 

“  ratio  of  aaount  of  fuel  left  in  tank  to  tank  capacity, 
and 

TC  ■  tank  capacity. 


Working  Losses 

VL  -  K3WPVK4/365  (16) 

where: 

WL  -  working  loss  (lb/day), 

K3  -  liquid-dependent  factor, 

•  2.25  a  10^  for  crude  oil, 

-  2.76  *  10-4  for  distillate  oil, 

•  2.95  *  10”*  for  kerosene, 

■  3.00  a  10”*  for  gasoline, 

-  3.24  a  10"4  for  JP-4, 

V  •  annual  gallons  of  liquid  puaped  into  tank,  end 
I4  •  turnover  factor  based  on  number  of  tank  turnovers/yr. 


Turnover  Factor 


*4  -  (30.34/TO)  ♦  0.137  (17) 

where: 

TO  ■  number  of  turnovers  per  year. 
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-  V/CT,  and 
Cj  -  tank  capacity, 
or  alternatively, 

K4  -  (30.34  Cj/V)  ♦  0.157  (18) 


9.9.1  Poaalble  Sources  of  Evaporative  Loaoee 

1.  Fixed-roof  tanka  —  breathing  losses,  working  losses,  and 
spillage. 

2.  Floating-roof  tanks  —  standing  storage  losses  and 

spillage. 

3.  Tank  trucks  —  breathing  losses,  working  losses,  and 

spillage. 

4.  Jet  aircraft  fuel  tanks  —  breathing  losses,  working 
losses,  and  spillage  (including  drainage  of  fuel  lines). 

5.  Piston-engine  aircraft  —  carburetor  and  fuel  tank  losses 
and  spillage. 

6.  Access  and  aircraft  service  vehicles  —  carburetor  and 

fuel  tank  losses  and  spillage. 

7.  Other  sources. 

M0B1LE2  aay  be  used  to  generate  carburetor  and  fuel-tank  losses  for 

various  vehicle  types  (sea  Sec.  9.8). 

In  addition  to  the  above  sources,  there  are  also  HC  evaporative  losses 

associated  with  a  variety  of  nonalrcraft-related  operations,  Including  dry 

cleaning,  paving,  and  spraying  and  finishing  of  surfaces.  One  additional 

aircraft-related  source  Is  deicing.  The  contribution  of  these  sources  Is 

swell,  but  the  user  way  wish  to  Include  such  sources  under  the  "OTHER  SOURCES'* 
category  for  cowpleteness. 


9.9.2  User  Input  Raoul rewsnts 

In  principle,  to  cowpute  evaporative-loss  rates,  the  uaer  aist  aupply 
weteorological  data,  Including  temperature  (T)v  wind  speed  (Vw),  and  diurnal 
temperature  variation  (AT),  and  complete  descriptions  of  every  source*  In 
practice,  there  are  far  too  many  sources  to  treat  each  one  on  an  Individual 


basis*  Furthermore,  since  dispersion  calculations  are  not  necessary  for 
evaporative  HCs  and  since  they  only  account  for  on  the  order  of  10Z  of  the 
total  HC  emissions  anyway,  the  locations  of  Individual  sources  are  not 
important*  Hence,  it  Is  recommended  that  the  Individual  sources  be  aggregated 
into  the  following  areas: 

i*  Storage-tank  farms. 

2.  Ground-vehicle  filling  and  service  stations,  including 
areas  where  tank  trucks  are  filled* 

3.  Ground-vehicle  parking  areas  (do  not  double-count 
evaporative  losses  from  access  vehicles). 

4*  Aircraft  parking  and  service  areas  (or  gate  areas)* 


For  each  of  these  four  areas,  the  following  Information  is  required: 


1  •  Storage-Tank  Farms  (working,  breathing,  standing,  and  spillage 
losses) : 


NL0C1 
X( I) , Y( I) 
NDIA(I) 

ma(i,j) 

FINFAC(I,J) 
ANNGALO ,  J) 

CAP( I  ,J) 
CONFAC(I,J) 
SEAL( I , J) 
NFXTKS( I ,J ,K,L) 


HFLTKS(I,J,K) 
SPILLI(I , J) 


no.  of  locations 

coordinates  of  center  of  farm  I 

no*  of  different  diameters  at  farm  I 

value  of  diameter  1  at  farm  J 

fixed-roof-finish  factor  of  diameter  I 

at  farm  J 

I  gal  pumped/yr  into  fixed-roof  tanks  of 
diameter  I  at  farm  J 

capacity  of  tanks  of  diameter  I  at  farm 

J 

construction  factor  of  floating-roof 
tanks  of  diameter  I  at  farm  J 
seal  factor  for  floating-roof  tanks  of 
diameter  I  at  farm  J 

no*  of  fixed-roof  tanks  of  construction 
type  I,  storing  liquid  J,  having 
diameter  K,  at  farm  t 

no*  of  floating-roof  tanks  storing 
liquid  I,  having  diameter  J,  at  farm  X 
#  of  gal  of  liquid  T  spilled  at  farm  J 
per  day* 


2*  Ground-Vehicle  Filling  and  Service  Stations  (working  and 
spill;  ge  losses) i 


HLOC2 

X(t)#Y(I) 


•  no*  of  locations 

•  coordinates  of  center  of  station  1 


A0'*l4'/  MM 


I MP  At.  T  (II  A  I  RtRAI  I  EMISSIONS  ON  AIR  QUAl  [I  Y  IN  INI 
V  U.  I N 1  M  01  AIRPORTS  V  IUI  AROONNI  NAIIONAl  L  AH  II 
INIHiiY  AND  I  NVIRONMI  NIAl  StSUMS  DIV..  .D  M  WOK  ,JAN  Ml 
l*i>  l  /.  .Ml  II  |i  I  AA  I  I  HA  1(1*11  AO  I  R  3  A  H)SSf>  l/i.u/y 


thruput(i,j) 


*  rf  £!#.  ■:  w  r 
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*  #  of  gal  of  liquid  of  typa  1  puaped  at 
station  J  par  day 
VTKFIL(ItJ)  -  #  of  lb  of  vapor  diaplacad  while  filling 
vehicle  tankaf  per  103  gal  of  liquid  of 
type  I  puaped  at  atatlon  J 

SPILL2(I,J)  -  #  of  lb  of  liquid  spilled,  per  103  gel 
of  liquid  of  type  I  puaped  at  atatlon  J. 


3.  Ground-Vehicle  Parking  Areas 


NLOG3 
X(I),Y(I) 
MVEH( I ,J) 

NHOTSKS(I,J) 


•  no.  of  parking  arena 

-  coordinates  of  parking  area  I 

-  no.  of  vehicles  of  type  I  parked  In  lot 
J  per  day 

•  no.  of  hot  soaks  of  vehicle  type  I  In 
lot  J  per  day 


4.  Aircraft  Parking  and  Service  Areas  or  Gate  Areas 


NLOC4 
X(I),Y(I) 
NACRP(I ,.!) 

AVRFRT(I,J,K) 

ACFIL(t,J) 

SPILL3(I,J) 


-  no.  of  aircraft  parking  areas 

■  coordinates  of  area  1 

•  no.  of  aircraft  of  type  1  refueled  In 
area  J 

■  average  #  of  gal  of  liquid  I  puaped  Into 
aircraft  type  J  In  area  K  per  day 

•  #  of  lb  of  vapor  displaced,  per  tO3  gal 
of  liquid  of  type  1  puaped  Into  an 
aircraft  of  type  J 

-  #  lb  of  liquid  aplllad,  per  103  gal  of 
liquid  of  type  I  puaped  Into  an  aircraft 
of  type  J 


9.9.3  Outline  of  Hydrocarbon  Evaporat lve-Loaa  Calculations 

Figure  12  shows  a  aacro  flowchart  of  a  procedure  that  can  be  used  to 
coapute  DC  losses  due  to  evaporation*  Care  should  be  taken  not  to  duplicate 
evaporative  losses  froa  vehicles  already  Included  under  aceesa  vehicles.  No 
detailed  flowcharts  have  bean  constructed  for  the  evaporative  loss  calcula¬ 
tions  because  of  the  uncertainty  of  the  actual  need  to  perfora  such  calcula¬ 
tions.  However,  If  the  need  arises,  such  a  flowchart  could  be  constructed  In 
a  staple,  straightforward  wanner  using  the  equations  and  user-supplied  Inputs 
defined  earlier  In  this  section. 


9.10  OTHER  SOURCES 

It  goes  alaoet  without  saying  that,  la  addition  to  the  aajor  source 
categories  discussed  in  Sec*  9.9,  there  are  nuastaus  other  sources  that  aay 


retrieve  hc  xv ap.  loss  data 


PILE  10 

BC  XV AP.  LOU  DATA 


■BAD  STORACE-TAIK  DATA 


STORAGE-TANK  INPUT  DATA 


COMPUTE  XMXUXOU  FROM 
STORAGE-TANK  FARM 


TANK-FARM  LOOP 


COMPUTE  FXXBD-ROOP- 
BREATHXN6  A  WORKING  1 


COMPUTE  FLOATINC-ROOP- 
TANK  STANDING  LOSSES 


ESTIMATE  SPILLAGE  IP  ANT 


I  CO  TO  NEXT  TANK  FARM 


COMPUTE  EMISSIONS  IKON 
VEHICLE  FILLING  STATIONS 


FILLING-STATION  LOOP 


COMPUTE  DISPLACED  VAPORS  BOEING  VEHICLE-TANK  FILLING 


A/C  »  aircraft 


ESTIMATE  SPILLAGE 


00  TO  NEXT  STATION 


COMPUTE  BOSSIOU  IRON  A/C  PARKING  A  EEPQELINB 


Flf .  12  INcW  Flowchart  of  NpAraoarboa 
Evapotat  1w-Uh  CaleaLatiooa 
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Fig*  12  (Coat'd) 


operate  In  or  near  an  aviation  facility.  Rather  than  burden  the  aodel  and  the 
coaputer  with  a  ayrlad  of  detalle  to  try  to  cover  all  of  theae  aource  types. 
It  aeena  aore  prudent  to  provide  the  near  with  the  option  of  being  able  to  odd 
additional  aoureee  for  which  ha  aaat  perfora  hie  own  Independent  eolesloo-rate 
eatlaatea  ualng  ealaalon  factors  froa  IH2.J  These  ealasioae,  together  with 
the  other  necessary  source-related  Input  data,  could  then  be  Incorporated  Into 
the  eaiseloa  aodele  to  be  accnanlated  with  other  source  ealeelona  or  to  be 
placed  in  the  special-source  ealaalon  Inventory  file  for  use  by  the  dispersion 
algorithae. 

The  structure  of  the  special-source  eaiseloa  inventory  file  Is 
described  In  Ref.  30,  along  with  these  of  all  the  other  source  files  required 
by  the  new  eoapotatlonal  procedare. 


10  SENSITIVITY  OP  EMISSIONS  CALCULATIONS 
TO  INCLUSION  OP  VARIOUS  SOURCE  TYPES 


10. 1  INTRODUCTION 

Because  Che  compilation  of  an  emission  inventory  at  a  military  or 
civilian  airport  la  a  manpower-intensive  undertaking,  it  is  Important  to 
recognize  which  emission  sources  are  most  important.  Unfortunately,  there  Is 
no  simple,  unique  answer  to  this  question.  The  relative  Importance  of  a 
particular  source  type  depends  on  the  pollutant  of  interest,  mode  of  operation 
of  the  source,  and  a  number  of  facility-dependent  factors,  including  facility 
geometry,  source  mix,  and  operational  procedures.  The  question  is  further 
complicated  by  the  purpose  to  be  served  by  an  emissions  calculation.  Por 
example,  whereas  the  overall  contribution  to  such  regional  Impacts  as 
visibility  degradation  and  photochemical  smog  formation  may  be  relatively 
significantly  Influenced  by  aircraft  airborne  operations,  the  ground-level 
concentrations  of  primary  pollutants  (CO,  NO^,  etc.)  are  more  Influenced  by 
nearby  ground-level  emissions  from  aircraft  (during  taxiing,  queuing,  etc*)* 
Hence,  in  seeking  guidance  on  how  best  to  expend  efforts  on  the  compilation  of 
an  emission  inventory,  the  purpose  for  the  calculation  should  be  given  proper 
consideration* 

With  these  considerations  in  mind,  it  Is  possible  to  develop  some 
general  guidelines  regarding  the  importance  of  including  certain  source  types 
or  modes  of  operation  on  the  net  emission  inventory*  In  the  sensitivity 

analysis  reported  here,  emphasis  Is  placed  on  the  percent  contribution  of 

individual  source  categories  and  modes  of  operation  to  the  overall  emissions 
from  an  airport  or  air  base*  It  is  also  possible  to  examine  the  Impact  of 
various  emission  source  terms  on  the  air  quality  (l*a*,  pollutant  concentra¬ 
tions),  although  considerably  more  effort  is  required  since  one  must  then 
consider  the  type  of  dispersion  algorithm  to  be  used,  the  meteorological  con¬ 
ditions,  the  positions  of  receptors  relative  to  sources,  the  mix  of  sources, 

and  the  levels  of  source  activity*  Although  no  air-quality  calculations  are 
used  in  the  present  sensitivity  analysis,  the  iaportance  of  emissions  from 
selected  source  types  to  air  quality  is  examined  with  the  help  of  a  combina¬ 
tion  of  linear  emission  density  calculations  (emissions  per  unit  length  of  a 
line  source)  and  other  straightforward  source-receptor  considerations* 

The  approach  used  here  primarily  involves  the  use  of  previously 
reported  results  (although  not  necessarily  reported  la  the  opea  literature)  of 
emission  calculations  for  a  number  of  civilian  and  military  facilities  using 
AVAR  and  AQAM,  respectively,  and  alee  the  use  of  simple  emission  estimates 
baaed  on  the  use  of  standard  LTD  cycles*  While  the  latter  is  much  simpler, 
the  former  provides  mare  realistic  estimates  of  the  importance  of  eomrcee  at 
actual  fad lit tea  because  measured  tlmes-ln-moda  are  used  for  each  facility* 


10.2  OVERALL  EMISSION  SOURCE  CONTRIBUTIONS  AT  CIVILIAN  AND  MILITARY 
FACILITIES 

Table  8  (how  eh*  result*  In  metric  toon  per  jrwr  (t/yr)  of  reported 
emission  calculations  at  four  aajor  civilian  airports  and  one  military  air 
baaa  using  the  AVAP  and  AQAM  amission  models,  respectively.  Tour  source 
categories  are  indicated:  aircraft,  aircraft  service  vehicles,  base  or  access 
vehicles,  and  stationary  sources*  This  table  clearly  demonstrates  the 
dependence  of  the  percent  contribution  of  a  given  source  clees  to  the  overall 
emleslons  on  pollutant  and  facility.  For  example,  baaa  end  access  vehicular 
traffic  make  up  21-59X  of  the  N0X  amissions.  It  should  be  noted  that,  with 
the  exception  of  the  Williams  AFB  data,  all  vehicle  amission  data  are  pre- 
1975.  Hence,  the  current  vehicular  emission  contributions  are  expected  to  be 
somewhat  smaller.  This  is  confirmed  by  a  study34  that  estimated  the  vehicular 
emlselona  to  be  45X  of  the  total  00  and  9X  of  the  total  N0X  produced  at  Dulles 
International  Airport  in  1976  and  36X  of  the  CO  and  7X  of  the  N0X  in  I960. 
Table  8  also  shows  that  whan  compared  with  THC  amissions  from  evaporation  and 
combustion,  evaporative  losses  vary  from  7X  to  28X  of  the  total  at  civilian 
airports  and  62X  of  the  total  at  Williams  AFB.  Of  this  62X  contribution  at 
Williams,  76X  was  due  to  spilling  during  aircraft  refueling  and  venting  of 
aircraft  fuel  lines. 

For  civilian  facilities,  the  aircraft  contribution  ranged  from  22X  to 
S8X  of  the  total  00  amissions  and  from  59X  to  78X  of  the  total  N0X 
amissions.  Hence,  aircraft  are  the  major  source  of  N0X  amissions  and  a 
significant  source  of  00  at  airports.  At  Williams  AFB,  which  is  a  training 
baaa  for  pilots  of  fighter  aircraft,  aircraft  contribute  72X  of  the  00  and  39X 
of  the  N0X*  Aircraft  also  contribute  the  major  portion  of  HCs.  Hence,  from 
the  point  of  view  of  photochemical  smog  precursors,  aircraft  are  the  major 
contributors  at  airports.  From  the  point  of  view  of  local  direct  pollutant 
impacts  (a.g.,  00  concentrations) ,  aircraft  are  significant,  hut  given  the 
proximity  of  the  public  to  vehicular  traffic,  the  latter  la  likely  to  be  a 
more  significant  contributor  to  local  ambient  00  levels. 


10.3  EMISSIONS  FROM  AIRBORNE  VS.  GROUND-BASED  AIRCRAFT  MOMS 

Am  ami salon-mo dal-lnde pendent  ami  facility-independent  way  of  enaminlng 
the  relative  laportanoa  of  Airborne  vs.  grouad-beeed  aircraft  emissions  la  to 
congers  the  percent  amissions  from  various  nodes  of  a  standard  LTD  cycle. 
Standard  LTD  eye  lee  ate  defined  for  military  aircraft  in  Refs.  35  and  36  and 
for  civilian  aircraft  in  Ref.  32.  A  standard  LTO  cycle,  which  nay  or  nay  not 
be  aircraft-type  specific,  consists  of  a  specification  of  the  engine  thrust 
settings  (or  nodes  or  foal  flow  rates)  and  the  times  agent  in  each  of  a 
ooRoones  of  several  aircraft  operational  nodes.  Table  5  in  Sac.  9.5  lists  the 
typical  military  and  civilian  aircraft  operational  modes.  Raf erases  35 
contains  tabulations  of  tines- la-nodsa  for  a  member  of  military  aircraft. 
Those  latter  data  are  based  an  data  gathered  at  air  bases.  Similar  data  are 
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Table  8  Pollutant  Ralaelone  by  Major  Source  Category 
at  Several  Major  Facllltlaa  (t/yr) 


Coabuetlon 

Evaporative 

Source 

CO 

BC 

HC 

N0X 

TSP 

SO 

Wlllleae  APB 

Aircraft 

1966 

253 

414 

50.3 

2.2 

24.4 

Service  Vehldaa 

60 

4 

— 

1.2 

1.0 

0.3 

Traffic 

568 

69 

67 

49.0 

24.0 

21.0 

Faellltiae 

150 

10 

62 

30.0 

1.0 

82.0 

Total 

2744 

336 

543 

130.5 

28.2 

127.7 

O'Hara  Airport 

Aircraft 

8279 

4888 

75 

2913 

493 

Service  Vehiclee 

3292 

735 

— 

211 

8 

Traffic 

3650 

624 

282 

468 

29 

Pacilltiaa 

IS 

15 

145 

139 

209 

Total 

15200 

6260 

500 

3730 

739 

Maahlngton 

Rational  Airport 

Aircraft 

1638 

388 

1096 

187 

Service  Vehielea 

1396 

193 

63 

5 

Traffic 

4443 

629 

412 

23 

PacilitiM 

- 

- 

- 

- 

Total 

7477 

1210 

475* 

1571 

215 

Atlanta  Airport 

Aircraft 

4959 

2415 

2072 

Service  Vahlclaa 

1626 

224 

57 

Traffic 

1870 

211 

212 

Facllltlaa 

136 

53 

313 

Total 

8600 

2900 

375 

2654 

*Fron  Mat*  26* 

Sour cm >  iMlllaaa  APB,  laf.  37;  O' Bare  Airport,  Baf.  6;  Maahlngtoo  Rational 
Airport,  uafabllabad  raaulta  obtained  aa  part  of  work  on  updated 
aaaaaanant  of  air  pollution  lapacta  at  aajor  alrporta  (ana  Baf* 
13);  and  Atlanta  Airport,  Baf*  31* 
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available  by  aircraft  claaa  for  civilian  aircraft  and  average  values  are 
reported  In  Kef*  32*  For  aany  purposes,  estlaates  of  aircraft  enlsslons  can 
be  based  on  the  use  of  such  LTO-cycle  data* 


10*3*1  Military  Aircraft  Operations 

For  purposes  of  cou paring  airborne  end  ground-based  enlsslons.  Table  9 
contains  a  tabulation  of  enlsslons  by  pollutant  and  selected  node  for  several 
military  aircraft  types.  Airborne  00  and  HC  enlsslons  constitute  only  2-7% 
and  1-9.4%  of  the  total  LTO-cycle  enlsslons,  respectively*  On  the  other  hand, 
airborne  N0X,  TSP,  end  S0%  enlsslons  constitute  29-5 IX,  28-54X,  and  23-36%  of 
the  total  LTO-cycle  emissions,  respectively*  Hence,  In  terns  of  mass  of 
enlsslons  only,  the  pollutants  fall  Into  two  groups  —  those  for  which 
airborne  enlsslons  are  relatively  unlnportant,  namely  00  and  HC,  and  those 
whose  airborne  enlsslons  are  relatively  lnportent,  nasal y  N0X,  TSP,  and  S0X* 

Somewhat  better  Insight  Into  the  question  of  which  aircraft  nodes  are 
most  lnportent  can  be  gained  by  examining  the  linear  emission  densities  asso¬ 
ciated  with  the  line  sources  used  to  represent  the  various  aircraft  nodes. 
The  linear  enlsslon  density  (p)  Is  defined  as  the  pollutant  enlsslons  per  unit 
length  of  the  line  source  (g/kn).  In  reality,  p  Is  raraly  independent  of 
poeltion  along  a  line  source.  This  Is  especially  true  for  runway  line  sources 
used  to  represent  the  takeoff  and  landing  sodas*  However,  for  present  pur¬ 
poses  It  is  satisfactory  to  consider  the  linear  enlsslon  densities  of  taxiway 
segments,  queuing  Unas,  and  segments  of  approach  and  departure  paths  to  be 
approxlnately  constant*  Given  this  assunptlon,  comparisons  can  readily  be 
made  between  the  linear  enlsslon  densities  of  these  line  sources*  It  Is  also 
worth  noting  that,  when  compering  total  enlsslons  per  node.  It  Is  necessary  to 
specify  the  total  times  spent  in  each  nods*  In  the  case  of  coopering  linear 
enlsslon  densities,  neither  the  tines  la  the  nodes  nor  the  lengths  of  the  Una 
sources  are  required*  Furthermore,  whereas  total  enlsslons  cannot  be  directly 
related  to  air  quality,  because  the  spatial  distribution  of  the  enlsslons  Is 
not  specified,  linear  enlsslon  densities  can  be  directly  related*  That  Is, 
lima  sources  having  the  sane  orientation  relative  to  a  receptor  will  produce 
roughly  the  sane  pollutant  concentrations  If  their  linear  enlsslon  densities 
are  equal,  asswuiag  plena  dynaaics  are  net  significantly  different.  Selow, 
cooperlsoas  will  first  be  nade  of  linear  emission  densities  without  regard  to 
lino  source  orientation*  Later,  the  affect  of  orientation  will  also  bo  con¬ 
sidered.  Before  proceeding  with  these  comparisons,  the  following  quantities 
should  be  defined.  Ift 

•  emission  rote  In  taxiing  mode  (g/s). 

If  •  taxiing  speed  (km/e),  and 


Pf  •  linear  end  salon  density  (g/km), 
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Table  9  Percent  of  Total  LTO-Cyde  Ealaelone  due  to 
Selected  Military  Aircraft  Modea  of  Operation 


Attack  Boaber  Fighter  Trainer 

Mode  A7  A37  B52D/F  C5A  F5  FI  5  T38 


CO  Kales lone 


Idle  at  Startup 

54 

26 

22 

28 

18 

31 

18 

and  Shutdown 

Taxiing 

37 

54 

53 

66 

59 

49 

65 

Cliabout 

0.3 

4 

0.6 

3 

4 

7 

4 

Approach 

2 

6 

3 

6 

7 

3 

6 

Total  Airborne 

2.3 

10 

3.6 

9 

11 

10 

10 

HC  Swiss Ions 

Idle  at  Startup 

57 

29 

23 

26 

21 

33 

20 

and  Shutdown 

Taxiing 

40 

61 

57 

63 

69 

53 

71 

C 11 about 

0.03 

0.7 

0.04 

0.2 

0.3 

0.4 

0.3 

Approach 

0.7 

3.0 

0.6 

5.6 

4.1 

9.0 

3.1 

Total  Airborne 

0.75 

3.7 

0.64 

5.8 

4.4 

9.4 

3.4 

HO  Emissions 

idle  at  Startup 

4 

11 

13 

4 

8 

11 

9 

and  Shutdown 

Taxiing 

3 

22 

27 

10 

26 

18 

31 

Cliabout 

23 

19 

23 

40 

19 

40 

21 

Approach 

9 

16 

15 

2 

10 

11 

20 

Total  Airborne 

32 

35 

38 

42 

29 

51 

41 

TSF  Emissions 

Idle  at  Startup 

10 

5 

12 

13 

5 

13 

5 

and  Shutdown 

Taxiing 

7 

10 

24 

32 

16 

21 

19 

Cliabout 

18 

26 

25 

23 

26 

33 

28 

Approach 

12 

16 

16 

3 

10 

14 

26 

Total  Air boras 

30 

42 

41 

28 

36 

47 

54 

SO.  Bales loos 

idle  at  Startup 

26 

16 

17 

16 

11 

17 

11 

and  Shutdown 

Taxiing 

18 

33 

38 

40 

34 

27 

52 

Cliabout 

11 

14 

11 

17 

15 

28 

18 

Approach 

13 

14 

13 

6 

10 

8 

8 

Total  Airborne 

24 

28 

24 

23 

25 

36 

26 
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then  t>j  -  Et/St.  Similar  definitions  follow  for  the  line  source  segments 
representing  other  aircraft  modes  of  operation.  To  compare  linear  emission 
densities  of,  for  example,  cllmbout  paths  to  taxlwaya,  one  can  simply  compute 
the  ratios  of  the  corresponding  linear  emissions  densities.  That  la: 

PC  f'C  /SC 
_  C1  Ct 

ci’«r  "“WT 

where: 


C|  -  1th  leg  of  the  cllebout  pathf 


pr  •  linear  emission  density, 

M 

Kr  -  emission  rate*  and 
C1 

-  average  aircraft  speed  on  the  1th  leg* 
ci 


Table  10  shows  the  values  of  the  ratios  of  N0X  linear  emission  densities  for 
several  airborne  line  sources  for  several  aircraft  types.  The  subscripts  C|, 
C2,  A| ,  and  A2  refer  to  the  first  and  second  legs  of  the  cllmbout  path  and  the 
first  and  second  legs  of  the  approach  path9  respectively*  For  Military 
aircraft  these  legs  or  line  source  segments  are  defined  aa  follows  (all  angles 
are  aircraft  dependent): 


1.  Approach  leg  #1  is  fro*  3000  ft  to  1000  ft,  and  leg  #2  la 
*rnm  1000  ft  to  ground  level* 

2.  Cliwbout  leg  #1  Is  fro*  ground  level  until  the  afterburner 
Is  shut  off,  and  leg  #2  is  frosi  afterburner  cutoff  to  3000 
ft*  Military  thrust  setting  is  used* 


Table  10  Ratios  of  W0X  Linear 
Emission  Densities 


Aircraft 

Pc/^T 

pC2,pT 

pA,,pT 

pA2/pT 

rS  and  T36 

1.123 

0.46 

0.123 

0.123 

T37 

0.71 

0.37 

0.14 

0.14 

fk 

2.36 

1.66 

0.31 

0.46 

Cl  306 

0.64 

0.44 

0.10 

0.06 

CU1 

2.47 

2.06 

0.67 

0.71 

s 


As  can  readily  be  seen  from  Che  results  in  Table  10,  the  N0x  linear 
emission  densities  for  airborne  aircraft  line  source*  range  froa  a  snail 
fraction  of  to  two  and  one-half  tlaes  as  large  as  the  taxlway  linear  emission 
density.  The  highest  and  second  highest  values  of  the  ratios  occur  for  the 
first  and  second  legs  of  the  cllmbout  paths,  respectively,  regardless  of 
aircraft  type.  Hence,  all  other  things  being  equal,  N0x  emissions  from 
cllmbout  path  line  source  segments  would  be  expected  to  have  comparable  air- 
quality  Impacts  to  NO^  emissions  from  taxlways.  Of  course,  all  other  things 
are  not  equal.  These  line  sources  are  generally  not  oriented  In  the  same  way 
relative  to  receptors.  Firstly,  and  most  Importantly,  cllmbout  paths  are,  of 
course,  inclined  at  an  angle  to  the  ground,  whereas  taxi way  emissions  are  at 
ground  level.  Secondly,  plume  dynamics  are  expected  to  be  somewhat  different 
for  low-speed  ground-level  and  hlgh-apeed  airborne  aircraft  plumes.  A  detailed 
dispersion  model  calculation  would  be  required  to  properly  determine  the 
effects  of  these  differences.  However,  one  can  easily  appreciate  the  fact 
that  since  airborne  plumes  must  grow  in  the  vertical  direction  in  order  to 
impact  ground-level  receptors,  their  Impact  at  ground-level  will  be  smaller 
than  for  plumes  emitted  at  ground  level*  Hence,  on  this  basis,  taxiway  NOx 
emissions  will  have  greater  Impacts  than  cllmbout  N0X  emissions  on  ground- 
level  receptors  equidistant  from  the  sources.  The  term  equidistant  is 
critical  here  because  cllmbout  paths  (at  least  their  lowest  portions)  could, 
in  principle,  pass  closer  to  the  public  than  taxlways.  Hence,  if  the  public 
resided  immediately  adjacent  to  the  lowest  legs-  of  the  cllmbout  paths,  the 
air-quality  Impacts  could  be  greater  for  those  paths  than  for  taxlways. 
Combined  with  runway  emissions,  the  cllmbout  leg  #i  emissions  may  not  be 
insignificant  in  such  situations.  How  Important  the  corresponding  air-quality 
impacts  would  be  would  require  calculations  with  a  dispersion  model  or 
measurements. 


10.3.2  Civilian  Aircraft  Operations 

The  following  discussion  la  baaed  primarily  on  emissions  computed  in 
connection  with  a  study  to  update  the  assessment  of  air-quality  Impacts  at 
several  major  airports.*  That  study  included  emission  and  alr-quallty 
calculations  for  a  one-hour  period  at  each  of  four  commercial  airports.  Table 
11  gives  the  times  and  runway  activities  used.  The  fractional  numbers  given 
for  numbers  of  arrivals  and  departures  arise  from  the  way  the  AVAP  model 
distributes  aircraft  to  the  runways  that  are  used  for  a  particular  wind 
direction*  The  particular  hour  used  for  each  airport  was  selected  on  the 
basis  of  high  (but  not  necessarily  peak)  aircraft  activity  and  the  probability 
that  "worst  case"  (l.a. ,  low  wind  spaed,  particular  wind  direction,  and  poor 
vertical  mixing)  meteorological  conditions  were  likely  to  occur.  Average 
times  in  modss  observed  separately  at  sack  airport  wars  used  in  these 
calculations. 

Aircraft  emissions  of  00  sod  M0X  for  several  airborne  and  ground -level 
modes  of  apodal  intaraat  ate  listed  ie  Tables  12  end  13  for  the  four  cl vi  lien 
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Table  11  Aircraft  Activity  Used  In  the  One-Hour  Ealaalon 
Ca leu let Iona  at  Cooaerclal  Airports* 


Airport 

Tine 
of  Day 

Wind 

Direction 

Quadrant 

Runway 

No. 

No.  of 
Arriving 
Aircraft 

No.  of 
Departing 
Aircraft 

DCA 

9-10aa 

1 

1 

20.1 

19.7 

2 

2.9 

3.3 

ORD 

S-9aa 

2 

2 

15 

* 

5 

15 

- 

6 

15 

- 

7 

- 

34.5 

8 

- 

34.5 

Jft 

7-8pn 

3 

1 

2.9 

8.4 

2 

5.8 

4.2 

7 

11.6 

12.6 

8 

8.7 

16.8 

LAX 

8-9an 

3 

l 

2.2 

0.8 

2 

7.3 

20.4 

3 

13.4 

24.5 

4 

5.2 

4.3 

*OCA  “  Waahington  National  Airport;  ORD  -  O'Hare  Inter¬ 
national  Airport;  JFK  ■  John  F.  Kennedy  International 
Airport;  and  LAX  ■  Los  Angeles  International  Airport. 


Table  12  Percent  of  Total  Aircraft  00  Ealsslons 
from  Selected  Modes  at  Coaaerdal  Airports 


Airport 

Queuing 

Taxiing 

Cltabout 

Approach 

DCA 

36 

23 

2 

9 

ORD 

48 

38 

1 

4 

JFK 

20 

10 

1 

4 

LAX 

39 

44 

1 

3 

100 


Table  13  Percent  of  Total  Aircraft  N0X  Emissions 
from  Selected  Modes  at  Commercial  Airports 


Alrpor 

Queuing 

Taxiing 

Cllmbout 

Approach 

DCA 

6 

4 

44 

13 

0RD 

8 

6 

45 

7 

JFK 

3 

9 

48 

6 

LAX 

5 

6 

51 

5 

► 


1 

4 

a 


airports.  For  convenience,  only  the  percentages  of  total  aircraft  emissions 
are  shorn.  As  already  Indicated  for  military  aircraft,  civilian  airborne 
aircraft  sources  of  CO  contribute  relatively  little  to  the  overall  aircraft  CO 
emissions  according  to  Table  12.  Queuing  for  takeoff  is  seen  to  be  a  major 
source  of  00  emissions.  This  Is  particularly  significant  because  queuing 
line  sources  are  substantially  shorter  than  the  combined  taxiway  line 
sources.  Hence,  the  linear  emission  densities  for  queuing  line  sources  are 
greater  than  for  taxlways,  and  It  follows  that  the  ground-level  alr-quallty 
Impacts  would  be  correspondingly  greater* 

Table  13  shows  that  N0x  emissions  are  greater  for  the  airborne  than  for 
the  ground-level  modes  with  cllmbout  contributing  the  major  fractions  (nearly 
50Z),  Tables  12  and  13  clearly  show  that,  from  the  point  of  view  of  emissions 
only,  the  airborne  sources  are  unimportant  for  00  emissions  but  Important  for 
NOx  emissions,  regardless  of  the  airport.  These  tables  also  show  that  the 
contributions  from  various  ground-level  aircraft  modes  vary  considerably  from 
airport  to  airport.  This  Is  largely  due  to  differences  In  airport  configura¬ 
tion  and  to  the  aircraft  activity  levels  used  for  the  emission  calculations 
(see  Table  11).  O'Hare  Airport  (0RD)  in  particular  has  a  relatively  large  CO 
contribution  from  queuing  due  to  the  large  number  of  departures  from  two 
runways.  More  will  be  said  about  queuing  In  Sec.  10.4  below. 

Although  there  are  differences  between  civilian  and  military  aircraft 
modes,*  the  same  arguments  regarding  linear  emission  density  and  airborne  vs. 
ground-level  sources  apply  and  will  not  be  repeated  here. 


•Only  one  approach  leg  (500  ft  to  ground  level)  Is  used  In  the  AVAP  model. 
The  first  cllmbout  leg  extends  from  ground  level  to  500  ft  and  the  second 
from  500  ft  to  2500  ft. 


10. 4  DEPENDENCE  OP  QUEUING  EMISSIONS  ON  THE  QUEUING  ALGORITHM 


The  queuing  algorithm  currently  uaed  in  the  AVAP  code  has  the  following 


form: 


q  "  180  *  Sff 


where  Tq  la  the  average  time  (hr)  spent  per  aircraft  In  a  departure  queue  and 
n^  la  the  number  of  departures  per  hour.  The  term  1/60  accounts  for  the  time 
spent  on  the  queuing  apron  when  no  aircraft  are  queued  up*  An  alternative 
expression  for  Tq  based  on  queuing  theory  was  discussed  In  Sec*  9*6*  Equation 
7  expressed  in  hours  rather  than  minutes  is; 


(20) 


where  Tq  1  is  the  time  (hr)  spent  per  aircraft  In  the  queuet  C  Is  the  runway 
capacity  (maximum  number  of  arrivals  plus  departures  that  could  theoretically 
be  serviced  by  the  runway  in  an  hour)t  V  is  the  aircraft  volume  (arrivals  plus 
departures)  during  the  hour,  and  the  term  1/60  has  been  added  to  account  for 
time  spent  on  the  queuing  apron  when  no  aircraft  are  queued  up*  The  diffl* 
culty  with  using  Eq.  20  Is  that  the  runway  capacity  must  be  specified*  As 
suggested  in  Sec*  9*6,  a  value  of  C  -  48  seems  to  be  reasonable,  given  the 
limited  data*  Table  14  gives  the  values  of  T^  and  Tq*  computed  from  Eqs*  19 
and  20,  respectively,  for  the  aircraft  activity  listed  In  Table  11*  It  can  be 
seen  from  these  results  that  the  queuing  times  per  aircraft  predicted  by  the 
AVAP  algorithm  (Eq*  19)  are  two  to  three  times  larger  than  for  Bq.  20  for  the 
aircraft  activities  given*  This  Is  expected  to  be  the  case  except  tdien  the 
traffic  volume  approaches  the  runway  capacity  (e.g*,  during  peak  traffic 
periods  or  periods  of  reduced  runway  capacity)*  In  fact,  as  can  be  seen  from 
Eqs*  19  and  20,  Tqv  will  be  >  Tq  when  the  following  condition  holds: 

C(f"Ty  -  75o  *  or  c<c  ”  v>  -  180  <2I) 


Por  C  *  48,  It  follows  that  Tq»  £  T-  when  the  total  number  of  arrivals  plus 
departures  per  hour  (V)  Is  greater  tnan  or  equal  to  44.25*  Two  hypothetical 
examples  are  Included  as  the  last  two  entries  In  Table  14*  In  both  of  these 
•uipUi t  th*  tot*i  traffic  voluat  la  *5.  Hanes,  '  >  Tq.  Furthermore,  with 
V  •  45  and  C  •  48,  It  la  aaan  that  tha  largar  th*  proportion  of  daparturaa  to 
arrival* ,  th*  largar  1*  th*  dlffarane*  between  T  '  and  T^.  It  la  not 
poaalbla,  given  th*  praaant  data  baa*»  to  aaka  a  definitive  determination  of 
which  algor 1 the  beat  represents  th*  real  world.  Th*  praaant  calculation*  only 
ahow  that  calculated  queuing  tlaa  1*  a  aanaltlv*  function  of  th*  algorithm 
uaed  and  that  emission*  da*  to  queuing  (which  are  proportional  to  queuing 
tlaa)  constitute  an  important  fraction  of  th*  total  ground-level  aircraft 
•missions. 
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Table  U  Queuing  Tinas  Computed  with  Two 
Alternative  Algorlthne  (C  “  48) 


Airport 

Runway 

No.  of 
Arrivals 
(n.) 

No*  of 
Doporturoo 
<nd) 

(hr) 

V 

(hr) 

DCA 

1 

20.1 

19.7 

0.126 

0.067 

2 

2.9 

3.3 

0.035 

0.018 

ORD 

7 

0 

34. 5 

0.208 

0.070 

JFK 

7 

11.6 

12.6 

0.087 

0.028 

8 

8.7 

16.8 

0.110 

0.032 

LAX 

3 

13.4 

24.3 

0.153 

0.067 

Hypothetical 

Ex.  1 

— 

20 

25 

0.156 

0.174 

Cxo  2 

10 

35 

0.211 

0.243 

10. 5  EVAPORATIVE  HYDROCARBON  EMISSIONS 

Evaporative  HC  anlaalona  are  anong  the  aoat  difficult  and  tedious 
emissions  to  estimate  and  compile,  and  are  generally  subject  to  large 
uncertainty.  Much  of  the  uncertainty  arises  from  the  fact  that  the  least  well 
quantified  sources  make  the  largest  contributions.  This  point  la  well 
Illustrated  by  soma  results  reported  for  the  Williams  APB  Source  Emission 
Inventory.  Table  IS  lists  several  sources  of  HC  emissions  Included  In  that 
inventory.  Of  the  543  t  of  HC  emissions  attributed  to  evaporative  loeaea,  1SX 
are  due  to  aircraft  fuel  tank  venting  and  38X  are  doe  to  spillage  during 
aircraft  fuel  tank  filling.  This  laaves  only  24Z  due  to  the  myriad  other  HC 
storage  and  handling  operations.  The  evaporative  loeaea  -attributed  to 
aircraft  are  quite  easy  to  compile  because  only  one  emission  factor  la 
required  to  specify  each  loss  per  operation.  In  contrast,  the  compilation  of 
the  looses  due  to  fuel  storage  and  handling  la  very  complex  (see  Sec.  9.9). 
On  the  other  hand,  the  one  number  needed  to  characterise  emissions  due  to 
aircraft  fuel  tank  venting  (B4  L/flll)  are  quite  uncertain.  Even  if  the 
actual  volumes  of  fuel  loot  were  representative,  much  of  this  fluid  may  and  up 
flowing  Into  drains  rather  than  being  evaporated  lata  the  atmosphere. 

la  terms  of  the  total  HC  emissions  duo  to  both  oaMbmetloa  aad  ovmpoio- 
tion,  the  category  representing  storage  ami  handling  eonstltetos  only  1SZ  of 
the  total.  Table  IS  gives  a  breakdown  of  this  latter  source  category.  It  com 
bo  seen  from  this  breakdown  that  tbs  working  loeaea  constitute  rnsnt  of  the 
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Table  15  Hydrocarbon  Balaalona  Inventory  Coaputtd  with 
AQAH  for  Wllllaae  Air  Force  Baee  (t/yr) 


Coabuetloa 

Evaporative 

Source 

Balsslons 

Loeeea 

Aircraft  Operation 

253 

Aircraft  Fuel  Venting* 

96.3 

Aircraft  Spillage* 

317.5 

Aircraft  Service  Vehicles 

3.7 

Base  Vehicular  Traffic 

6B.9 

Base  Facilities 

10.4 

Fuel  Storage  and  Handling 

129 

(including  parked  vehicles) 

Total 

336 

543 

Grand  Total  ■  879 

*Assuaes  2  l  vented  per  arrival  and  per  departure  for 
•oat  aircraft. 


*Assuaes  4-t  apt 1 lag*  per  fillup  for  aoet  aircraft. 
Source:  Kef.  37. 


Table  16  Breakdown  of  Hydrocarbon  Evaporative  Loeeea  frow  Fuel  Storage 
and  Handling  Other  than  Aircraft  Venting  and  Spillage  Loeaee  (t/yr) 


Fixed -Koof 

FixedHtoof 

Floating* 

Tanka, 

Tanks, 

Koof  Tanka, 

Source 

Working 

Breathing 

Standing 

Spillage 

Other 

Storage  Tanka* 

33.4 

0.12 

2.6 

Filling  and  Handling 
Fetrolaua  Storage 

20.3 

0.002 

3.3 

Farhad  Tank  Trucks 

1.7 

Farked.  Vehicles 
Others* 

67.1 

1.0 

Total 

53.7 

66.9 

2.6 

3.3 

1.0 

Grand  Total  •  129 


*2B  tanks. 

Varts  cleaning,  paint,  and  thinner. 
Seureet  Kef.  37. 
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storage  tank  loaaes  (33.4  t/yr  for  storage  ♦  20.3  t/yr  for  filling  opera** 
tlons).  Breathing  losses,  by  contrast,  are  relatively  saall.  Evaporation 
froa  parked  vehicles  and  tank  trucks  (68.8  t/yr)  constitutes  the  largest 
source  category. 

How  do  the  various  source  contributions  at  Wllllaaa  AFB  coapare  with 
those  at  coaaerclal  airports?  Tables  17  and  18  show  the  results  of  evapora¬ 
tive  loss  calculations  for  O'Hara  International  Airport.  These  tables  show 
that  evaporative  lossea  constitute  only  7X  of  the  total  HC  ealasiona  due  to 
coabustlon  plus  evaporation.  Host  of  the  evaporative  losees  are  due  to  fuel 
storage  and  handling  (Including  evaporation  froa  parked  vehicles)  at  O'Hare  In 
contrast  to  Vllllaas  AFB,  where  76X  were  due  to  aircraft  fuel  venting  and 
spilling.  At  O'Hare,  any  spilled  fuel  Is  assuaed  to  run  down  the  dralna  In 
the  paveaent.  The  storage  tank  braathlng  lossea  are  cpaparable  to  the  working 
losses  at  O’Hare  In  contrast  to  Vllllaas  where  the  working  losees  are 
predoainant.  However,  In  agreeaent  with  Vllllaas,  evaporation  froa  parked 
vehicles  constitutes  the  largest  source  of  ealasiona  In  the  fuel  etorage  and 
handling  class. 

Hydrocarbon  ealasiona  data  for  other  airports  are  very  Halted. 
Evaporative  losses  froa  fuel  etorage  and  handling  at  Vashiogton  National 
Airport  (DCA)  is  reported  to  constitute  HZ  of  the  total  HC  cessions.  °  No 
breakdown  la  available. 

At  Atlanta  Airport  (LAX)  fuel  storage  ealasiona  were  reported  to  aake 
up  1 IX  of  the  total  HC  eaissions.^  Again  no  breakdown  was  reported. 


Table  17  Hydrocarboo  Ealasiona  Calculated  for 
O'Hare  International  Airport  (t/yr) 


Source 

Coabustlon 

Ealsslons 

Evaporative 

Losses 

Aircraft  Operation 

4909 

Aircraft  Pilling* 

75 

Aircraft  Service  Vehicles 

735 

Access  Vehicles 

624 

Airport  Facilities 

15 

Foal  Storage  and  Handling 
(Including  parked  vehicles) 

428 

Total 

Grand  Total  -  6788 

6283 

503 

*Any  spillage  is  sssaasd 


to  flow  down  drains. 
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Table  18  Hydrocarbon  Evaporative  Loaaaa  at 
O'Hara  Intarnatlonal  Airport  duo  to  Pual 
Storago  and  Handling  ( t  /  yr  ) 


Plsed-Koof 

Fixed-Roof 

Floating- 

Tank a. 

Tanka, 

ftoof  Tanka, 

Working 

Breathing 

Standing 

Storage  Tanks 

78 

34.3 

22.1 

Service  Vehicle 
Pilling 

Parked  Vehicles 

11.4 

282 

(carburetors) 

Total 

89.4 

318.3 

22.1 

Grand  Total 

-  428 

Sourca:  tef.  6* 


It  aaana  claar  fro*  tha  above  r« ported  raaulta  that  NC  evaporative 
loaaaa  f roa  foal  atoraga  and  handling  conatltutaa  a  relatively  snail  fraction 
of  tha  total  HC  aalaalona  froa  coaaarclal  airports  and  that  a  substantial 
fraction  of  thaaa  evaporation  loaaaa  at  both  nllltary  and  coaaarclal  facili¬ 
ties  la  froa  parked  vehicles.  Only  In  tha  case  of  HI  1  Haas  API  ware  aircraft 
venting  and  a pi 11 age  aajor  aoureaa  of  HC  asdsalona. 
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II  MNUIT  AND  RCOOHHEMOATIONS 


The  first  part  of  this  report  addressed  the  statue  of  the  AVAP  model 
and  IQW  frea  the  pa re pec tire  of  the  modeling  requlreaente  of  users  concerned 
with  alr-gasllty  problems  in  aviation.  The  spproech  used  wee  to  begin  with  a 
brlof  daocrlptloo  of  the  typos  of  probleoo  likely  to  be  encountered  In  both 
the  ollltary  ood  del 11 an  sectors,  followed  by  s  detailed  discussion  of  those 
characteristics  of  the  prohlsas  that  determined  the  technical  requirements  for 
the  applicable  computational  procedures  or  models.  This  was  followed  by  a 
discussion  of  the  operational  or  user  requirements  of  the  models.  A  review 
ami  evs  last  lam  of  the  AVAP  model  end  AQAM  was  then  presented.  In  which  the 
Intended  usee,  strengths,  sad  weaknesses  were  described.  Finally,  the  methods 
mead  by  the  AVAP  model  sad  AQAM  to  treat  various  aspects  of  the  emission  or 
diapers  loo  calculations  sore  compared,  and  the  bast  methods  were  selected,  or 
altermmtiuee  warm  roc ones ad ad  where  appropriate. 

The  latter  pert lea  of  the  report  addressed  the  future  of  the  AVAP  modal 
ami  AQMk  Prso  the  evaluation  given  In  the  first  pert  of  the  report.  It  was 
clear  that  future  efforts  were  required  In  a  number  of  areas.  Including: 

a  Incorporation  of  technical  improvements. 

e  Updating  of  documentation  to  remove  inconsistencies  with 
uereleme  of  the  computer  codes  currently  in  use. 

e  Bectiflcatieo  of  the  problem  of  usability  of  the  computer 


because  of  the  number  of  interrelated  problems  sad  decisions  required, 
a  systematic  approach  to  ths  problem  was  developed.  The  “decision  tree”  that 
resulted  la  regarded  os  a  first  stop  towards  systematically  laying  out  which 
alternatives  ere  evailable  sad  whet  decisions  ere  required.  This  device 
should  at  least  simplify  the  decision-making  process  by  clarifying  the  types 
of  tasks  that  naturally  felloe  from  various  alternative  paths. 

Iho  final  section  of  the  tent  was  devoted  to  on  outline  of  a  proposed 
asm  computational  system  that  should  alleviate  at  least  some  of  ths  problems 
Identified  in  earlier  sections.  Iho  objectives  were  paramuot  In  the  new 
designs  to  mobs  the  model  easier  to  use  end  to  bo  able  to  implement  ths  model 
am  modern,  smell  computers.  In  the  process  of  designing  the  mew  eye  tea,  it 
one  found  that  net  only  could  these  objectives  be  met,  but  that  in  earns 
respects  the  model  could  even  be  Improved  technically.  Only  the  emission  or 
f  reef  end  pert  lea  of  the  now  system  has  been  add  reseed  in  this  report, 
bof arenas .  30  contains  detailed  flowcharts  and  other  Information  needed  to 
guide  the  development  of  the  actual  computer  codes  tor  ths  amloelea  portion  of 
the  model,  the  design  of  ths  dispersion  portion  of  the  new  system  remains  to 
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Th*  rnalndtr  of  this  ooctloa  concern*  rtcoinJitlow  for  future  | 

work*  Thooo  recoenendatlons  or*  divided  Into  two  group*:  those  Chat  portal*  1 

to  application-type  taaka,  ouch  as  th*  davalopeant  and  tasting  of  coapatar 

coda*  for  th*  n*w  coaputatlonal  ayatea*  and  tho*a  that  coocorn  futura  KM 
•ffort*.  Th*  overriding  racouuandat Ion  1*  that  both  types  of  effort*  b* 

pursued  in  parallel  to  avoid  future  ahortcoulnga  in  either  usability  or 
technical  quality* 

Vlth  respect  to  the  appllcat lone-related  Issues*  It  is  strongly 
recooaeuded  that  the  new  coaputatlonal  package  be  adopted  in  a  fora  aiailar  to  . 

that  outlined  In  See*  9.  In  particular*  it  la  laportant  to  coapletely  ' 

separate  th*  aaleelona  coaputatlons  froa  th*  dispersion  coaputatioa*.  It  Is  j 

also  laportant  to  sake  use  of  th*  nodular  nature  of  th*  cod*  and  th*  data-f 11a  1 

structure*  Thee*  structures  will  greatly  reduce  th*  coaputer-cora-atoraga  and  I 

run-tlaa  requlraasnts.  Further,  every  effort  should  be  aed*  to  aelntain  th* 
sane  overall  structure  In  both  the  allitary  aod  civilian  versions  of  the  ^ 

ealsslon*  portion  of  th*  systea*  It  would  be  desirable  to  use  the  seas  ■ 

prograaalng  language  for  both  versions*  If  possible*  | 

Before  the  ealsslon  portion  of  th*  new  systea  can  be  fully  laplaaented,  j 

several  saall  tasks  should  be  undertaken*  These  Include  th*  following:  j 

1*  Service-vehicle  ealsslon  factors  should  be  updated* 

2.  The  propoeed  queuing  algor itha  should  be  evaluated  using 
available  or  new  data  froa  allitary  and  coaaerclal 
facilities*  Queuing  data  has  been  taken  at  VI 1  Haas  APB 
and  at  eaveral  eoaaercial  airports*  1 

j 

3*  Aircraft  eatssion  factor*  should  be  updated,  and  new 
aircraft  types  should  be  added* 

i 

A.  Approach-  and  departure-path  peraaeter*  should  be  selected 
as  defaults*  Soae  sensitivity  tests  regarding  th*  lapect 
of  airborne  aircraft  operations  on  ground-level  concen¬ 
trations  should  be  conducted* 

J 

Frelialnary  desiga  of  th*  dispersion  portion  of  the  new  coaputatlonal 
systea  should  bs  undertaken  as  early  as  possible,  so  that  it  can  be  ready  for 
lnpleuentatlon  soon  after  the  ealsslon  portion.  The  dispersion  package  will 
require  a  driver  cods  that  can  read  in  the  aateorologlcal  data,  read  In  the 
source  data  (preferably  one  source  at  a  tiae),  select  the  appropriate 
dispersion  algetic  has,  end  output  the  results*  It  is  recoaasoded  that  the 
dispersion  package  be  aa  analyse  of  th*  boot  features  of  bach  the  AFAF  aodel 
and  AQ4H  aad  that  it  incorporate  roflaeasots  already  developed  bet  not 
previously  lapleasated. 


_ L 
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••for*  completing  the  design  of  the  now  dispersion  package.  It  la 
sided  that  the  following  Ml)  taaka  ha  eoaplstsd: 


5.  Compare  conpu tat Iona  ualng  Brubaker' a  advanced  single- 
aircraft  -plums  reaearch  model1'  with  treatments  of  lint* 
aource  ewlaalone  uaed  by  the  AVAP  model  and  AQAN. 
Determine  under  what  condition*  the  more  aophlatlcated 
modeling  approach  la  required. 


8.  Utilise  the  high  time-resolution  data  obtained  at  O' Ha re 
Airport,  together  with  raflnamanta  of  the  data  analysis 
techniques  developed  In  connection  with  the  O'Hare 
program,1  to  obtain  optimum  values  of  parameter*  needed 
to  define  aircraft  takeoff -plume  behavior. 


7.  Critically  review  the  surface  observations  at  Dulles  and 
Washington  national,  together  with  the  tower  measurements 
at  Dulles,  to  determine  to  what  extent.  If  any,  better 
definitions  of  aircraft  pluaw  behavior  can  be  obtained  for 
the  following  aircraft  modes:  taxiing,  queutng,  and  take¬ 
off.  Include  an  evaluation  of  che  types  of  analyses 
performed  and  the  degree  of  completeness  of  the  analysis. 
(Hot  all  data  were  analysed  in  detail.  Sow  data  received 
very  little  attention.)  Also  determine  If  the  data 
themselves  contain  the  necessary  characteristic*  to 
warrant  the  use  of  more-sophisticated  analysis  techniques. 


8.  On  the  basis  of  the  result*  of  taaka  5-7  above,  evaluate 
the  need.  If  any,  for  additional  field  axperiawats 
designed  to  further  elucidate  aircraft  plums  behavior. 


9.  Critically  review  all  data  related  to  plume  dynamics. 
Including  effects  of  plums  rise,  initial  plums  d lepers Ion, 
enhanced  vertical  dispersion,  wind  direction  relative  to 
aircraft-exhaust  velocity,  aircraft  mode  of  operation, 
etc.  Prepare  recommendations  for  paramtterltlng  these 
effect*  and  for  additional  anelyaee  of  existing  data  or 
for  acquisition  of  new  data.  Incorporate  parameter ise- 
tlona  Into  the  Joint  Air  Quality  Hod* ling  Peckagn. 

10,  Critically  review  the  advantages.  If  any,  of  incorporating 
recent  AM  recommendation*  regarding  alternative 
atablllty-claaelflcatioo  schemas  and  dispersion  parame¬ 
ters.  Consider  additional  na  neurons  wt  burdens  for 
research  and  regulatory  applications. 
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11.  Develop  en  alternative  to  the  Nosakl  equation  aa  a 
default  for  estimating  the  hourly  average  mixing  depth. 

Acquire  appropriate  data  aeta  from  representative  altea 
around  the  country  for  testing  purposes. 

After  completion,  each  portion  of  the  new  systea  should  be  carefully 
tested  and  then  the  portions  should  be  coablned  and  tested  together.  It  is 
recoaasnded  that  the  new  systea  be  coapared  with  either  the  original  AVAP 
aodel  and  AQAM  or  with  data  obtained  froa  the  field  program*  at  Wllllaas  ATI. 
Washington  National  Airport,  or  O'Hare  International  Airport. 

In  addition  to  designing  and  lap lamenting  the  new  coaputatlonal  systea, 
it  is  slso  laportant  to  continue  BAD  efforts  to  develop  the  state  of  the  art 
of  aircraft  air-quality  modeling  and  to  investigate  phenomena  that  were  beyond 
the  scope  and  capabilities  of  the  original  AVAP  aodel  and  AQAM.  Therefore,  it 
is  recoaasnded  that  the  following  tasks  be  considered  for  future  work: 

12.  Conduct  a  field  prograa  Involving  aeasureaents  of  NO  and 
NOj  concentrations  at  various  tlaes  of  the  year  (espec¬ 
ially  during  the  high-Oj  season),  using  two  or  aore 
monitoring  sites  to  obtain  a  aore  coaplete  characterisa¬ 
tion  of  the  NO /NO 2  problea.  Include  examination  of 
effects  due  to  aultlevent  pluae  interactions  and  to 
taxiing  and  queuing  aircraft.  Emphasis  should  be  placed 
on  deteralnlng  potential  peak-to-aean  ratios,  since  SPA 
will  probably  proaulgate  annual-average  standards  that 
are  presumed  to  protect  the  public  against  short-tera 
peak  concentrations. 

13.  Conduct  further  theoretical  studies  of  NO  and  NO2 
concentrations  to  deteraine  if  peak  NOj  froa  aircraft 
could  poae  potential  short-tera  health  hasards.  Examine 
the  applicability  of  the  existing  NO/NOj  conversion  data 
to  other  civilian  and  to  military  facilities.  Evaluate 
differences  between  civilian  and  military  aircraft,  if 
any,  expected  as  a  result  of  differences  in  NO2  to  NO 
eaission  ratios. 

14.  Conduct  a  field  and  laboratory  prograa  to  continue  pre¬ 
liminary  efforts  at  collection  and  characterisation  of 
organic-pollutant  saaploa  at  various  locations  around 
both  civilian  and  ailltary  facilities.  Baaaios  both  in- 
pluas  and  aablent  saaplas  froa  various  source  types. 


19*  further  elucidate  photochemical  reaction 
Involving  aircraft  hydrocarbon  ealssioos. 
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16.  Conduct  theoretical  reactive-pollutant  atudlee  In  support 
of  both  the  field  and  laboratory  atudlee  of  organic 
components  of  aircraft  exhaust  to  determine  their  Impact 
on  photochemical-smog  formation. 

17.  Conduct  experimental  studies  to  determine  the  relation¬ 
ship  between  aircraft  emissions  end  odors. 

18.  Investigate  the  applicability  of  certain  components  of 
the  new  dispersion  package  to  specialised  military 
applications.  Examine  ways  in  which  that  utility  could 
be  enhanced  to  fully  exploit  the  flexibility  of  the 
modeling  package. 
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